НОВЫЕ СООТНОШЕНИЯ ДЛЯ АНАЛИТИЧЕСКИХ РЕШЕНИЙ ЛОКАЛЬНО-НЕРАВНОВЕСНОГО ТЕПЛООБМЕНА Текст научной статьи по специальности «Математика»

УДК 539.3

DOI: 10.18698/1812-3368-2023-6-4-24

НОВЫЕ СООТНОШЕНИЯ ДЛЯ АНАЛИТИЧЕСКИХ РЕШЕНИИ ЛОКАЛЬНО-НЕРАВНОВЕСНОГО ТЕПЛООБМЕНА

Э.М. Карташов

professor.kartashov@gmail.com

РТУ МИРЭА, Москва, Российская Федерация

Аннотация

Развиты обобщенные модельные представления локально-неравновесного теплопереноса в терминах теории нестационарной теплопроводности для уравнений гиперболического типа. Модель включает в себя одновременно три системы координат: 1) декартовы координаты-массивное тело, ограниченное плоской поверхностью (одномерный случай); 2) сферические координаты-массивное тело с внутренней сферической полостью (центральная симметрия); 3) цилиндрические коорди-наты-массивное тело с внутренней цилиндрической полостью (радиальный поток теплоты). Рассмотрены три вида интенсивного нагрева (охлаждения): 1) температурный; 2) тепловой; 3) нагрев средой. Приведены примеры локально-неравновесного теплообмена, имеющего волновой характер, с учетом конечной скорости распространения теплоты. Волновой характер выражается наличием ступенчатой функции Хэвисайда в аналитическом решении базовых задач для частично ограниченных областей. Построены изохроны для температурных функций. Показано, что на поверхности фронта бегущей волны температурный профиль имеет разрыв. Это приводит к задерживанию оттока теплоты за границу разрыва — характерная особенность аналитических решений волновых уравнений. Приведены функциональные конструкции аналитических решений базовых задач локально-неравновесного теплообмена в области с движущейся во времени границей — практически новый случай в аналитической теплофизике. Представленные результаты имеют практическое применение

Ключевые слова

Обобщенные математические модели, локально-неравновесный теплообмен, области с подвижной границей

Поступила 25.01.2023 Принята 13.06.2023 © Автор(ы), 2023

Введение. В последние годы усилился интерес к изучению различных локально-неравновесных процессов, связанных с переносом энергии (теплоты, массы, импульса или их аналогов). С одной стороны, это связано с естественным направлением развития науки — от равновесных в целом систем к локально-равновесным и далее к локально-неравновесным. С другой стороны, интенсификация различного рода технологических процессов, использование материалов со сложной структурой (полимеров, жидких кристаллов, капиллярно-пористых и других дисперсных систем), возможность получения сверхвысоких (или сверхнизких) температур и давлений объясняют с практической точки зрения необходимость математического описания систем и физических процессов в экстремальных, локально-неравновесных условиях. Изучение таких систем может быть основано на различных версиях локально-неравновесной термодинамики: кинетических, молекулярно-динамических, феноменологических. Наибольшее развитие получила феноменологическая теория Максвелла — Каттанео — Лыкова — Вернотта, учитывающая инерцию теплового потока, связанную с конечной скоростью распространения теплоты в твердых телах. С точки зрения модельных представлений эта теория приводит к краевым задачам для уравнений гиперболического типа. Математические модели в виде краевых задач для уравнений гиперболического типа называют краевыми задачами обобщенного типа, которые значительно отличаются от классических на основе феноменологии Фурье по сложности их решения. Их специфика заключается в относительной простоте исходных математических моделей и трудности решения в аналитически замкнутом виде. Отсюда весьма незначительные успехи в нахождении точных аналитических решений таких задач и в основном для частично ограниченной области в декартовой системе координат при граничных условиях первого рода (температурный нагрев). Наиболее приемлемый метод их решения — операционный, но он приводит к аналитическим решениям в пространстве изображений по Лапласу в виде сложных функциональных конструкций, оригиналы которых не содержатся в справочниках по операционному исчислению. На этом пути возникают серьезные трудности вычислительного характера.

Рассмотрена открытая проблема аналитической теплофизики в терминах обобщенной модели локально-неравновесного теплообмена. Исследована проблема тепловой реакции на нагрев (охлаждение) массивного тела, ограниченного изнутри плоской поверхностью (упругое полупространство в декартовой системе координат), цилиндрической поверхностью (упругое пространство в цилиндрической системе координат с внутренней цилин-

дрической полостью) или сферической поверхностью (упругое пространство в сферической системе координат с внутренней сферической полостью). Рассмотрены такие случи интенсивного нагрева (охлаждения) границы тела (области), как температурный, тепловой, средой. Следует отметить, что развитие обобщенных гиперболических моделей, затрагивающих многочисленные практические приложения, — одно из направлений аналитической теории теплопроводности, которое не получило развития до настоящего времени. Рассмотрены принципиально новые случаи локально-неравновесного теплообмена в области с движущейся во времени границей, которые не были описаны в литературе.

Цель работы — описать функциональные конструкции в качестве аналитических решений обобщенных краевых задач нестационарного теплопереноса одновременно в трех системах координат и в условиях нагрева (охлаждения).

Обобщенные модельные представления для локально-неравновесного теплообмена. Классические модели аналитической теории переноса теплоты обязаны появлением линейно-градиентного соотношения

которое Ж.-Б. Фурье предложил в докладе «О распространении тепла в твердом теле» (Париж, 1807). Окончательно Фурье представил теорию в 1822 г. в работе «Аналитическая теория тепла», которую У. Кельвин назвал «великой математической моделью» [1]. Совместно с уравнением энергии для изотропных твердых тел

сР дТ (М, *) = - div д(М, г)+¥(М, г) дг

закон Фурье приводит к уравнению параболического типа для нестационарного теплопереноса вида

и соответствующим для (1) краевым задачам с начальным и граничными условиями

q (M, t) = -?igradT(M, t),

ST(M, t) dt

aAT(M, t) + — F(M, t), M e D, t > 0,

cp

(1)

T(M, t) t = 0 = 00(M), M eD,

(2)

ßx ÖT(M, t) + ß2T(M, t) = рэф(M, t), M e S, t > 0.

(3)

Здесь D — конечная или частично ограниченная выпуклая область изменения М(х, у, г); 5 — кусочно-гладкая поверхность, ограничивающая область D; п — внешняя нормаль к 5 (вектор, непрерывный в точках 5). Входящие в (1)-(3) параметры — теплофизические характеристики среды — постоянные величины в интервале значений температуры, не выходящих за точки перехода [2]. В литературе неоднократно отмечались некоторые парадоксы при использовании модельных представлений (1)-(3):

- отсутствие инерционности процесса теплопроводности в законе Фурье ввиду неучета механизма переноса теплоты элементарными частицами вещества (электронами, молекулами, ионными решетками) и времени релаксации, связанного со временем свободного пробега микрочастиц и, как следствие, вытекающего из аналитических решений моделей (1)-(3) вывода о бесконечной скорости распространения теплоты;

- сингулярный характер теплового потока и скорости движения изотерм в области х > 0, t > 0 при х ^ 0 (£ ^ 0).

Тем не менее перечисленные парадоксы не ограничивают область применения математических моделей (1)-(3), охватывая все новые содержательные объекты и все большее число самых разнообразных приложений [3-8]. В последние годы усилился интерес к изучению процессов, протекающих в локально-неравновесных условиях, что обусловлено широкими возможностями их практического применения [9-12]: создание новых технологий получения наноматериалов и покрытий с уникальными физико-химическими свойствами (бинарные многокомпонентные металлические сплавы, полимерные материалы, металлические полупроводниковые стекла, наножидкости, коллоидные био- и криосистемы); оптимизация режимов лазерной обработки изделий; режимы интенсивного нагрева (охлаждения) компонентов наноэлектроники и нанотехники; нагрев, плавление и абляция вещества при воздействии сверхкоротких лазерных импульсов и др. Интенсификация тепловых процессов потребовала для их описания уточнения гипотезы Фурье, что и было проделано в рамках учета локальной неравновесности, заложенной в соотношении

Я (М, t) = -к grad Т(М, t) - тг , (4)

дt

которое учитывает конечную скорость распространения теплоты. Здесь хг — время, являющееся мерой инерции теплового потока, связанной

со скоростью распространения теплоты соотношением vт = -у/ а / хг. На необходимость учета влияния ограниченности скорости переноса теплоты (массы) указывали Дж. Максвелл в теории газодинамики [13],

А.В. Лыков при исследовании тепло- и влагопереноса в капиллярно-пористых телах [14], К. Каттанео [15] и П. Вернотт [16] в теории теплопроводности. Уравнение энергии и соотношение (4) приводят к уравнению переноса теплоты гиперболического типа

er (м, t) а 2т (м, t) ^

= aAT (M, t) -xr----ъ

dt dt2

(5)

cp

&(M ■ t) + ± р (м, t)

dt x

r

и соответствующим краевым задачам нестационарной теплопроводности обобщенного типа. При математической постановке указанных задач следует применять соответствующие локально-неравновесные граничные условия. Использование стандартных локально-равновесных граничных условий (3) (довольно часто имеет место в публикациях по аналитической теплофизике) может привести к физически противоречивым результатам (например, к появлению отрицательных решений для температуры [10]). Эти вопросы детально рассмотрены в [17]. Сформулированы корректные обобщенные граничные условия на основе соотношения (4) в интегральной и эквивалентной дифференциальной формах. Так, в первом случае для условия теплового нагрева (охлаждения) граничное условие второго рода имеет вид

— } дТ(И Т) I Ие з ехр [ - —] йх = +(Я0 / (г), г >0. (6)

хг о ОП V Тг )

В случае нагрева средой следует записать

г

Здесь

1 t dT (M, т) | ( t-тЛ — J---1MeS expI--Idx =

Xr Q ОП v Tr J

h {T(M, t)\Me S-[To + S+ (t)(Tc - Tq ]} , t > Q. (7)

f1, t > 0,

S+ (t) = \ (8)

I 0, t < 0.

При (1 / к) ^ 0 (к ^ го) из (7) имеем граничное условие при температурном нагреве

Т(м, г)| Ие Б = Т + (г)(ТС - Т0), г > 0. (9)

Следует отметить, что гиперболическое уравнение (5) для описания локально-неравновесных процессов тепломассопереноса получено впервые в работах И.А. Фока [18] и Б.И. Давыдова [19] на основании предположения о конечном значении скорости частиц, переносящих энергию или массу. Уравнение (5) также получил А.С. Предводителев [20], исходя из анализа скоростей перемещения изотермических поверхностей с использованием представлений Римана, т. е. при полном отказе от релаксационной формулы (4).

Далее рассмотрим практические случаи нестационарного теплообмена, которые затем объединим в обобщенную модель.

В первом случае в декартовых координатах (х, у, г) рассмотрена область г > Я, г > 0, ограниченная плоской поверхностью, температурное состояние которой описывается функцией Т (г, г) (г = 1,2,3), удовлетворяющей условиям нагрева (охлаждения):

дТ (г, г) д2Т (г, г) д2Т- (г, г) о

—--= а----%г-;—, г > Я, г > 0, (10)

дг дг2 дг2

Т(г, г)|г=0 = Т), дТ^, г) | г = 0 = 0, г > Я, |Т(г, г)\ <ю, г > Я, г > 0, (11)

Т1(г, г)| г=Я = Тс, г > 0, (12)

1 г дТг(г, т), ( г-IV _ 1

— ]---1 г = Я ехр |--| йх = + — qo, г > 0, (13)

%г 0

dz 1 v тг ) X

— } дТз(г, Т) I г = я ехр {-0 = к\Тз(г, г)\г = Я - Тс ], г > 0. (14)

%г 0 дг 1 ^ %г ) 1 J

При Тс > Т0 происходит нагрев, при Тс < Т0 — охлаждение.

Во втором случае в сферических координатах (р, ф, 9) рассмотрена область р > Я, г >0 с внутренней сферической полостью при нагреве в условиях центральной симметрии Т = Т (р, г):

dt

■ = a

v ар2 р dp j

DTi(p, t) _ i d2Ti(p, t) , 2 DTi(p, t) 1 a2Ti (p, t)

-Tr-^^, P> R, t> 0, (15)

dt2

Ti (p, t) t=0 = T0, M-Bl^ t=0 = 0, p> R, | Ti (p, t)\ <«>, p> R, t > 0,(16) dt

H(p, t)| p = r = Tc, t > 0, (17)

1 jöT2(p, т)

Zr о

öp

p = R

exp

t-тЬ -1

- I dт = + — qo, t > 0,

ir 1 к

(18)

1 j 5Ts(p, T)

Zr 0

5p

p = R

exp

dт = h [Ts(p, t)| p = r - Tc ], t > 0. (19)

В третьем случае в цилиндрических координатах (г, ф, г) рассмотрена область г > Я, t > 0 с внутренней цилиндрической полостью в условиях радиального изменения температуры Т = Т (г, t):

STi (r, t) dt

= a

f д2Ti (r, t), i dTi (r, t) 1 a2Ti (r, t) 0

-Tr-;—, r > R, t> 0, (20)

dr2

r dr

dt2

Ti(r, t)| t=0 = T0, r > R,

dTi (r, t)

t = 0

= 0, r > R, \Ti(r, t)\ <ю, r > R, t > 0, (21)

Ti(r, t) r = R = Tc, t > 0,

JöT2(r, T)

r 0

j. 5Ts(r, T) Tr 0 dr

dr

exp

r = R

t-Tb 1

- I dт = + — q0, t > 0,

zr j X

(22)

r = R

exp

t-T

dт = h [Tj(r, t)|

r = R - Tc |, t > 0. (23)

Подобные случаи интенсивного нагрева (охлаждения) поверхности области (реального тела) представляют значительный практический интерес [21]. Так, в случаях поверхностного диэлектрического нагрева, расчета термических напряжений в стенках цилиндров паровых машин и двигателей внутреннего сгорания, в теории автоматических систем регулировки температуры, при исследовании области звуковых частот металлов при высоких или очень низких значениях температуры поверхности, в многочисленных случаях резкой смены температуры поверхности космических авиационных объектов, в машиностроительной отрасли при работе на различных экспериментальных установках для определения температурного состояния образцов и других представляет интерес охватить одновременно три случая во всех системах координат в рамках обобщенной модели. Это может представлять несомненную практическую значимость, например, в теории теплового удара. Во избежание излишней громоздкости введем следующие безразмерные переменные.

В области г > Я, г > 0:

5=У£, х=±, ^, р=;>, Bi•=ккa, а а а уТ V р

щ = Т(г,г)-[(1 -у)Т0 +утс], . = 1,3, (24)

^ ' (Тс -Т0)(1 -2у) ' ' ' ' '

Т) = Т2(г,г)-Т0 , г = 2. ^0/ Х)(а / V р)

Здесь Vр =у]20(1 - V)/ [р(1 - 2у)] — скорость распространения волны расширения в упругой среде, близкая к скорости звука. В области р > Я, г > 0:

а а а ут ур

х) =Т (р;!>-'I1 )Т2 +) ^' .( = 1,3, (25)

(Тс - Т0)(1 - 2у)

ЩЬ т) = Т2(р,г)~Т0 ,,= 2. (^0 / А, )(а / V р)

В области г > Я, г > 0:

^, , ^=М, .вi-=iкa,

а а а уТ V р

шп ^ Т(г, г) -[(1 -у )Т0 +уТс ] ,

т) =-№ 0 ч-1, г = 1,3, (26)

(Тс - Т0)(1 - 2у)

т) = Т2(г,г)" Т0 ,,= 2. ^0/ Х)(а / V р)

В (24)-(26) у = 0 при нагреве, у = 1 при охлаждении.

Обобщенная математическая модель задачи и ее аналитическое решение. Запишем обобщенную модель задачи, используя (10)-(26):

дWi(Ъ т) х) + 2т +1 х) 2 дт) = 0

дх д^ Эх2 ,

^>^0, т>0, (27)

У dWtä, х)

0,

W&т)|т = 0= —Ц-, i = 1,3, ^

V-v — 1

2у -1 дх

= 0, i = 1,2,3, (28)

т = 0

0,i = 2, ^0, |Wi&т)|<«>, ^>^0, T>0,

1 "У

Mfe t)|

1 - 2y

, т> 0,

1 Tfaw2(^, T')|

i—ж—N=^0 exP ß2 0 ^

( x-x'^

ß2

d x' = -

2y-1

т> 0,

1 iw&X) I

j-ж-N = ^0 exP

ß2 0 ^

( X~X^

d x' =

= Bi*

W3& x)|

5 = ^0

1 "У

, x> 0.

1 - 2 у_

Здесь г = 1 при т = —1/2; г =2 при т = 1/2; г = 3 при т = 0. В пространстве изображений по Лапласу

ж

Щ & р) = { Щ (£, т)ехр(-рх)й х

о

решение тепловой задачи (27)-(32) имеет вид

W (£, р) = ■

■fi(^0, р)Кт р)] + Xi

(1 - 2 у)

где Ц(р) = >/ß2р2 + р; Х1 = 1; Х2 = 0; Хэ = 1;

1 . . - Ш р)

р(2 у-1)

_ , i = 1; fi- = _

рКт [№р)Г ' i р2Кт+1 [Шр)]

_Bi* Ц( р)

i = 2;

р {рКт+1 [Мр) ] + Bi* |Д(р)Кт р) ]}

, i = 3.

(29)

(30)

(31)

(32)

(33)

i = 1,2,3, (34)

(35)

Здесь Кт (г) — модифицированные функции Бесселя. Обратное преобразование в (34), (35) весьма трудоемко и приводит к сложным и трудно обозримым выражениям. Если учесть, что в практических задачах наибольший интерес представляет начальное время теплового процесса, при котором проявляются, например, инерционные эффекты (в соответствующих термических напряжениях), то в (34), (35) рассмотрим большие значения р, что соответствует в оригиналах малому значению времени т. Для этого используем в (34), (35) асимптотическое представление цилиндрических функций для больших значений р:

1~т2

р) = ехр(-р) (V > 0).

л/р

В результате операционное решение обобщенной задачи примет вид

У

Wi(5' p)=p)exp JP (pЬ>(2 у-1)

Ф»(Scp) = -, i = 1, Ф»(Scp) p2 + p-, i = 2, P P2

Ф. (Sc p) = ■

Bi*j!( p)

(36)

(37)

i = 3.

р [ р + ВГ |д(р)]

Переходя к оригиналам, получаем искомые функциональные конструкции аналитических решений гиперболических задач локально-неравновесного теплообмена одновременно в трех системах координат в условиях нагрева (охлаждения):

Х) = ^У -1/2 „

exp

2ß

^0

2ß ' eXP

2р (4-40)Р

2ß2

V У

(1 - 2 у)

Ц 71 [^2 -ß2 L

Vx'2 -(^-^с)2ß2

d х' х

2 у-1

W2(5. х) = х

(1 - 2 у)

(38)

f [х-(S-S^ßjexp

2ß J

S"S0

0R i f (^')exP

(^)ß

f x' ^

. 71 f-(S-Ы2ß2

v 2ß2 ,

2ß2

Vx'2 - ß2

-dx' X

(39)

*0

f (т) = ß exp

W3& т) =

^ J f (T-x')exp

2p (4-40)ß

1

. 2ß2

(W ^)т+1/2 (1 - 2 у)

2ß2) ß

+ — J exp

2ß2

I0

rjL ,2ß2

d x';

f [t-(^-^0)ß]exp

2ß

. I1 T'2 -(^0)2ß2

2ß

2ß2

Vx'2 -(M0)2ß2

■dx'Y

(2y-1)

1 ю I

f (x) = 1 --^ J Ш^) j exp(-ß2x')

Bi*

Bi*

exp

Г x'

Ф*

Ш*2

V Bi'y

(40)

d x'.

Здесь 1У (г )(у = 0,1) — модифицированные функции Бесселя; ]0(г) — функция Бесселя; г|(г) — функция Хэвисайда; Ф*(г) = 1 -Ф(г), Ф(г) — функция Лапласа.

Кривые температурного, теплового нагрева (у = 0) и нагрева средой в сечении £ = 2, рассчитанные по (38)-(40) при т = -1/2, = 0,

Р = 1, В1 * = 0,3, приведены на рисунке.

Wtä, т)

0,8 \ ... т

0,6 - V

0,4 -

0,2 -

0 1

Изохроны, рассчитанные по формулам (38) (а), (39) (б) и 40 (в), Bi* = 0,3

Согласно аналитическим решениями (38)-(40), особенностью процесса теплопроводности в рамках математических моделей для уравнения (27) является волновой характер, что выражается наличием ступенчатой функции Хэвисайда в аналитическом решении задач. В любой момент времени существует область теплового следа < т и невозмущенная область (см. рисунок). В точках области да), находящихся на

расстоянии от границы более чем = т, изменение температуры не происходит и ее значение равно начальному. На поверхности фронта распространяющейся волны <^ = т и на фронте температурный профиль имеет разрыв, что приводит к задержанию оттока теплоты за границу разрыва. Амплитуда разрыва на фронте волны быстро затухает с увеличением времени нагрева, несмотря на это поступающий в область тепловой поток от границы увеличивается. В области за фронтом тепловой волны в направлении текущей пространственной переменной наблюдается существенное различие решений уравнений гиперболического и параболического типов: из решений (38)-(40) получаем начальное условие, а решения для параболического уравнения приводят к гладким решениям, существенно большим начального условия.

Аналитические решения гиперболических моделей переноса в нецилиндрических областях. Изложенные здесь вопросы относятся к числу новых в аналитической теплофизике и до настоящего времени практически не рассматривались. Здесь приведены гиперболические модели в области с движущейся во времени границей (нецилиндрические области). Вопросов, при рассмотрении которых приходится сталкиваться с краевыми задачами нестационарной теплопроводности в нецилиндрических областях типа [0, у(г)], г > 0, или [у(г), ю], г > 0, где у(г) — непрерывная функция, весьма много. Подобные вопросы возникают при теоретическом изучении процессов переноса энергии, связанных с изменением агрегатного состояния вещества, в теории прочности, теории плотин, механике почв, термике нефтяных пластов, электродинамических задачах, в задачах фильтрации, в теории зонной очистки материалов, кинетической теории роста кристаллов, в термомеханике при изучении теплового удара и т. д. [3, 22, 23]. С математической точки зрения краевые задачи переноса в области с движущимися границами принципиально отличны от классических. Вследствие зависимости границы области от времени к этому типу задач неприменимы классические методы уравнений математической физики, так как, оставаясь в рамках этих методов, не удается согласовать решение уравнения теплопроводности

с движением границы области. Этим объясняется тот факт, что в аналитической теплофизике к настоящему времени рассмотрены лишь простейшие случаи с равномерно движущейся границей или частично с корневой зависимостью и только для случаев на основе классической феноменологии Фурье.

Рассмотрим тепловые модели вида

дт(г, о д2Т(г, о д2Т(г, г)

8t Ti (z, t)

• = a-

— tr

t=0 - T0

dz2 dt2

ST(z, t)

-, z > R + vt, t > 0,

dt

t=0 = 0, z > R, i = 1,2,3,

|Ti(z, t)|<да, z > R + vt, t >0, Ti(z, t)| z=r+vt = Tc, t >0,

j ÖT2(z, T)

Tr 0 8z

z = R +vt

exp

t-тЬ -1

- I dт = + — q0, t > 0,

1r 1 к

(41)

(42)

(43)

(44)

(45)

j ÖT3(z, T)

Tr 0 dz

z = R + vt

exp

t-T

d т =

= h

T$(z, t)

- T

z = R + vt 1c

, t > 0.

(46)

Последовательно проведем следующие преобразования. Введем безразмерные переменные

ц =

V pz

, Т =

^ ^0, , ы*=ha

V0 =■

a

a

vT

_ ( ) Ti(z,t)-[(1 -y)Tp + YTc] . 13

©i (л, T) =---1, i = 1,3,

' (Tc -Tc)(1 - 2y)

T2(z, t) - T0

, i = 2.

(^0 / Ша / V р)

Здесь у = 0 при нагреве, у = 1 при охлаждении. Далее введем подвижную систему координат ^ = л^0т, полагая (л, т) = Щт). Соотношения (41)-(46) преобразуются следующим образом:

= (1 _у0Р2) +У0 ВЩх) + 22 йо_

Эг 0 Я;2 0 сЧ д£дт

"ß2

d2Wi & т)

Wi & т)|

т = 0-"

ÖX2

У

, i = 1,3

£>£0, Т>0, i = 1,2,3, dWi ß, т)

2у -1 5т

1 "У

0, i = 2, W1& т)| = 1

т = 0

т> 0,

= 0, ^>^0,

1 TfdW2ß, т')

ß2 0

5 = ^0

exp

т-т

- 2у

1

d т' = -

± Tfaws(^, то

ß2i ^

exp

5 = ^0

т-т

2y-1

d т' =

-, 0,

(47)

(48)

(49)

(50)

= Bi*

W3& x)

1 "У

£ = £0 1 - 2y |W& t)| ^>^0, T> 0.

, T> 0,

(51)

(52)

Переходим в пространство изображений по Лапласу (33) и последовательно находим операционные решения модельных представлений (47)-(52):

W1&р) =

1

х exp

р(1 - 2 у)

V0ß2(^~^0) 1 "V(2ß2

exp

(р + 2a*)(р + 2ß*)

1 "V(0ß2

р

exp

(v0/2)(^-^0)

1 -v2ß2

У

W2& р) =

(1 - 2 у)

exp

1 -v0ß2

V( р + 2a*)(р + 2ß*)

exp

(1 + ß2 р)(1 -v0ß2)

ßV (р + 2a*)(р + 2ß*) + (V0/2) + V0ß2j

exp

р(2у-1)

V0ß2(^-^p)

1 -v0ß2

(У0/2Х5-Ы 1 -v0ß2

(53)

р

, (54)

W3&р) =

Bi*(1 ~V()ß2) 1 - 2y

exp

(р + 2a*)(р + 2ß*)

1 - v00ß2

exp

V0ß2(^~^0)

1-v0ß2

р

(1 + ß2 p)

p {ßV(p + 2a*)(p + 2ß*) + (vc / 2) + voß2p + Bi*(1 + ß2p)(1 - v2ß2)}'

(vo/2)(S-So)'

x exp

2r2

1 -v|Cß

p(2y —1)

Здесь

1+ У1 -y0ß2 ; 1 -yj1 -v0ß2

Uv л J P _

4ß2 4ß2

Переходим последовательно к оригиналам. Температурный нагрев (охлаждение):

(55)

(56)

Ш, т) =

1 - 2 у

exp

(vo/2)(5-5o) (1 -V(2ß2)

х ( exp

2ß(1 -v2ß2)

voß2(^4o)

1-voß2

2ß2

.,2 (Mo)2ß2

(1 -v^ß2)2

x-

(S-So)ß

1 -Voß

(57)

Тепловой нагрев (охлаждение):

^ x) = ^^p

(1 - 2 y)

(vo/2)(^-^o)

2r2

1 -voß

x <

f

x-

(S-^o)ß 1 + Voß

exp

2ß(1 -v2ß2)

+ ^ x 2r2

W1 -voß

1-V2ß2

- i f

(4-4o)ß

1-V2ß2

т voß2Q;4o) .

т---T

1 -v0ß2

exp

rdL ,2ß2

т,2 _ (Mo)2 ß2

V (i -v0ß2)2

Здесь f (x) = 1 - (1/2)exp(-x /(2ß2)). Нагрев (охлаждение) средой:

Л

т--

(Mo)ß 1 -Voß

W3& T) =

1 - 2y

Mo

/ (

W f x-

V

voß2(^4o)

1-vgß2

i f

(4-4o)ß

1-vgß2

(JLM'

1 + Voß .

exp

(Mo)

2ß(1 -voß2)

voß2(Mo) 1 -vooß2

-x

exp

2ß2

.,2 (Mo)2ß2

i (1 -v2ß2)2

dx'

T-

(Mo)ß'

1 -voß ,

Здесь

f * (T) = 7~~T~ + A2 exp ; Ao = 2ßBi* (1 - voß);

1 + Ao ^ A1)

A1 =

(Ao + 2)ß2 1 + Ao

; A2 =

An

1 + Ao

'E 1'

A1

(58)

(59)

При нахождении оригиналов изображений (54), (55) использовано приближенное равенство ^(р + 2а*)(р + 20*) = [р / (2р2)] (2р2 +1/ р).

В соответствии с соотношениями (57)-(59) в области с движущимися границами нахождение аналитических решений гиперболических моделей локально-неравновесного теплообмена предполагает длительную процедуру вычислений на основе операционного исчисления. При этом все преобразования доводятся до конца. Конечные результаты достаточно удобны для численных расчетов, учитывая хорошо развитое в настоящее время программное обеспечение. Здесь просматривается волновой характер распространения теплоты с задержкой начала процесса (о чем свидетельствует наличие функции Хэвисайда в аналитических решениях). Таким образом, функциональные конструкции (57)-(59) в качестве аналитических решений гиперболических краевых задач локально-неравновесного теплообмена могут быть использованы в качестве приложений для решения многих задач в науке и технике.

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

ЛИТЕРАТУРА

[1] Карташов Э.М. Новые операционные соотношения для математических моделей локально-неравновесного теплообмена. Российский технологический журнал, 2022, т. 10, № 1, с. 7-18. DOI: https://doi.org/10.32362/2500-316X-2022-10-1-68-79

[2] Карташов Э.М. Аналитические методы в теории теплопроводности твердых тел. М., Высш. шк., 2001.

[3] Карташов Э.М., Кудинов В.А. Аналитические методы теории теплопроводности и ее приложений. М., URSS, 2018.

[4] Лыков А.В. Теория теплопроводности. М., Высш. шк., 1967.

[5] Зарубин В.С. Инженерные методы решения задач теплопроводности. М., Энергоатомиздат, 1983.

[6] Тихонов А.Н., Самарский А.А. Уравнения математической физики. М., Наука, 1977.

[7] Формалев В.Ф. Уравнения математической физики. М., URSS, 2021.

[8] Sobolev S.L. On hyperbolic heat-mass transfer equation. Int. J. Heat Mass Transf., 2018, vol. 122, pp. 629-630.

DOI: https://doi.org/10.1016/j.ijheatmasstransfer.2018.02.022

[9] ^динов И.В., ^динов ВА. Maтемaтическaя модель локально-неравновесного теплопереноса с учетом пространственно-временной нелокальности. Инженерно-физический журнал, 201S, т. 88, № 2, с. 398-40S.

[10] Kудинoв ВА., Еремин A^., Kудинoв И.В. Разработка и исследование сильно неравновесной модели теплообмена в жидкости с учетом пространственно-временной нелокальности. Теплофизика и аэромеханика, 2017, № б, с. 929-93S.

[11] ^рсанов ЮА., Kиpсaнoв A^. Об измерении времени тепловой релаксации твердых тел. Известия РАН. Энергетика, 201S, № 1, с. 113-122.

[12] Синкевич O.A., Семенов A.M. Решение уравнения Больцмана методом разложения функции распределения в ряд Энскога по параметру ^удсена в случае наличия нескольких масштабов зависимости функции распределения от времени и координат. Журнал технической физики, 2003, т. 73, № 10, с. 1-S.

[13] Maxwell J.C On the dynamical theory of gases. Phil. Trans. Royal. Soc., 1967, vol. 1S7, no. 1, pp. 49-88. DOI: https://doi.org/10.1098/rstl.1867.0004

[14] Лыков A^. Теплопроводность и диффузия в производстве кожи, заменителей и других материалов. M., Гизлегпром, 1941.

[15] Cattaneo C. Sulla Conduzione de Calore. Atti del. Seminario Matematico e fisico dela Universita di Modena, 1948, vol. 3, pp. 3-21.

[16] Vernotte P. Les paradox de la theorie continue de l'equation de la chaleur. C. R. Acad. Sci., 19S8, vol. 246, no. 22, pp. 31S4-31SS.

[17] ^рташов ЭЖ. Aнaлитические решения гиперболических моделей теплопроводности. Инженерно-физический журнал, 2014, т. 87, № S, с. 1072-1081.

[18] Фок ИА. Решение одной задачи теории диффузии по методу конечных разностей и приложение его к диффузии света. Л., ГОИ, 192б.

[19] Давыдов Б.И. Диффузионное уравнение с учетом молекулярной скорости. ДАН СССР, 193S, № 2, с. 474-47S.

[20] Предводителев A.C Проблемы тепло- и массопереноса. M., Энергия, 1970.

[21] Зарубин В.С., Kувыpкин Г.Н. Maтемaтические модели термомеханики. M., ФИЗMATЛИT, 2002.

[22] Aттеткoв A^., Волков И.К Температурное поле анизотропного полупространства с подвижной границей при локальном тепловом воздействии в условиях теплообмена с внешней средой. Тепловые процессы в технике, 2018, т. 10, № 1-2, с. S6-61.

[23] Arre^m A^., Беляков Н.С., Волков И.К Влияние подвижности границы на температурное поле твердого тела с цилиндрическим каналом в нестационарных условиях теплообмена с внешней средой. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2006, № 1 (20), с. 31-40.

Карташов Эдуард Михайлович — д-р физ.-мат. наук, профессор кафедры высшей и прикладной математики ИТХТ им. М.В. Ломоносова РТУ МИРЭА (Российская Федерация, 119454, Москва, пр-т Вернадского, д. 78).

Просьба ссылаться на эту статью следующим образом:

Карташов Э.М. Новые соотношения для аналитических решений локально-неравновесного теплообмена. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2023, № 6 (111), с. 4-24. БОТ: https://doi.org/10.18698/1812-3368-2023-6-4-24

NEW RELATIONS FOR ANALYTICAL SOLUTION OF THE LOCAL NON-EQUILIBRIUM HEAT TRANSFER

E.M. Kartashov

RTU MIREA, Moscow, Russian Federation

professor.kartashov@gmail.com

Abstract

The paper develops generalized model representations of the locally non-equilibrium heat transfer in terms of the theory of non-stationary heat conduction for the hyperbolic type equations. The model simultaneously includes three coordinate systems:

1) Cartesian coordinates-a massive body bounded by a flat surface (one-dimensional case); 2) spherical coordinates-a massive body with an internal spherical cavity (central symmetry); 3) cylindrical coordi-nates-a massive body with an internal cylindrical cavity (radial heat flow). Three types of intensive heating (cooling) are considered: 1) temperature;

2) thermal; 3) media heating. Examples of locally non-equilibrium heat transfer of a wave nature are provided taking into account the finite speed of heat propagation. The wave character is expressed by the Heaviside step function in analytical solution to the basic problems for the partially bounded regions. Isochrones for the temperature functions were constructed. It is shown that the temperature profile is having discontinuity on the surface of the traveling wave front. This leads to a delay in the heat outflow beyond the discontinuity boundary, which is a characteristic feature in analytical solutions to the wave equations. Functional constructions of analytical solutions to basic problems of the locally non-

Keywords

Generalized mathematical models, locally non-equilibrium heat transfer, regions with a moving boundary

equilibrium heat transfer are presented in a region

with the moving boundary, it is a practically new Received 25.01.2023

case in analytical thermophysics. The presented Accepted 13.06.2023

material has wide practical application © Author(s), 2023

REFERENCES

[1] Kartashov E.M. New operational relations for mathematical models of local nonequi-librium heat transfer. Russian Technological Journal, 2022, vol. 10, no. 1, pp. 7-18 (in Russ.). DOI: https://doi.org/10.32362/2500-316X-2022-10-1-68-79

[2] Kartashov E.M. Analiticheskie metody v teorii teploprovodnosti tverdykh tel [Analytical methods of heat conduction theory of solid body]. Moscow, Vysshaya shko-la Publ., 2001.

[3] Kartashov E.M., Kudinov V.A. Analiticheskie metody teorii teploprovodnosti i ee prilozheniy [Analytical methods of heat conduction theory and its applications]. Moscow, URSS Publ., 2018.

[4] Lykov A.V. Teoriya teploprovodnosti [Theory of thermal conductivity]. Moscow, Vysshaya shkola Publ., 1967.

[5] Zarubin V.S. Inzhenernye metody resheniya zadach teploprovodnosti [Engineering methods for solving heat conduction problems]. Moscow, Energoatomizdat Publ., 1983.

[6] Tikhonov A.N., Samarskiy A.A. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1977.

[7] Formalev V.F. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, URSS Publ., 2021.

[8] Sobolev S.L. On hyperbolic heat-mass transfer equation. Int. J. Heat Mass Transf., 2018, vol. 122, pp. 629-630.

DOI: https://doi.org/10.1016/j.ijheatmasstransfer.2018.02.022

[9] Kudinov I.V., Kudinov V.A. Mathematical simulation of the locally nonequilibrium heat transfer in a body with account for its nonlocality in space and time. J. Eng. Phys. Thermophy, 2015, vol. 88, no. 2, pp. 406-422.

DOI: https://doi.org/10.1007/s10891-015-1206-6

[10] Kudinov V.A., Eremin A.V., Kudinov I.V. The development and investigation of a strongly non-equilibrium model of heat transfer in fluid with allowance for the spatial and temporal non-locality and energy dissipation. Thermophys. Aeromech., 2017, vol. 24, no. 6, pp. 901-907. DOI: https://doi.org/10.1134/S0869864317060087

[11] Kirsanov Yu.A., Kirsanov A.Yu. About measuring the thermal relaxation time of solid body. Izvestiya RAN. Energetika [Proceedings of the Russian Academy of Sciences. Power Engineering], 2015, no. 1, pp. 113-122 (in Russ.).

[12] Sinkevich O.A., Semenov A.M. Solution of the Boltzmann equation by expanding the distribution function with several time and coordinate scales in the Enskog series in Knudsen parameter. Tech. Phys., 2003, vol. 48, no. 10, pp. 1221-1225.

DOI: https://doi.org/10.1134/L1620111

[13] Maxwell J.C On the dynamical theory of gases. Phil. Trans. Royal. Soc., 1967, vol. 157, no. 1, pp. 49-88. DOI: https://doi.org/10.1098/rstl.1867.0004

[14] Lykov A.V. Teploprovodnost i diffuziya v proizvodstve kozhi, zameniteley i dru-gikh materialov [Heat conduction and diffusion in the production of leather, substitutes and other materials]. Moscow, Gizlegprom Publ., 1941.

[15] Cattaneo C. Sulla Conduzione de Calore. Atti del. Seminario Matematico e fisico dela Universita di Modena, 1948, vol. 3, pp. 3-21.

[16] Vernotte P. Les paradox de la theorie continue de l'equation de la chaleur. C. R. Acad. Sci., 1958, vol. 246, no. 22, pp. 3154-3155.

[17] Kartashov E.M. Analytical solutions of hyperbolic heat-conduction models. J. Eng. Phys. Thermophy., 2014, vol. 87, no. 5, pp. 1116-1125.

DOI: https://doi.org/10.1007/s10891-014-1113-2

[18] Fok I.A. Reshenie odnoy zadachi teorii diffuzii po metodu konechnykh raznostey i prilozhenie ego k diffuzii sveta [Solving the problem of diffusion theory by the finite difference method and its application to light scattering]. Leningrad, GOI Publ., 1926.

[19] Davydov B.I. Diffusion equation with molecular velocity. DAN SSSR, 1935, no. 2, pp. 474-475 (in Russ.).

[20] Predvoditelev A.S. Problemy teplo- i massoperenosa [Problems of heat and mass transfer]. Moscow, Energiya Publ., 1970.

[21] Zarubin V.S., Kuvyrkin G.N. Matematicheskie modeli termomekhaniki [Mathematical models of thermomechanics]. Moscow, FIZMATLIT Publ., 2002.

[22] Attetkov A.V., Volkov I.K. Temperature field of anisotropic half-space with movable boundary at local thermal impact in conditions of heat exchange with outer ambient. Teplovye protsessy v tekhnike [Thermal Processes in Engineering], 2018, vol. 10, no. 1-2, pp. 56-61 (in Russ.).

[23] Attetkov A.V., Belyakov N.S., Volkov I.K. Influence of boundary mobility on temperature field of solid body with cylindrical channel under non-stationary conditions of heat exchange with environment. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2006, no. 1 (20), pp. 31-40 (in Russ.).

Kartashov E.M. — Dr. Sc. (Phys.-Math.), Professor, Department of Higher and Applied Mathematics, RTU MIREA (Vernadskogo prospekt 78, Moscow, 119454 Russian Federation).

Please cite this article in English as:

Kartashov E.M. New relations for analytical solution of the local non-equilibrium heat transfer. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2023, no. 6 (111), pp. 4-24 (in Russ.). DOI: https://doi.org/10.18698/1812-3368-2023-6-4-24