CC BY
108
УДК 517.958+536.2(075)
МОДЕЛИРОВАНИЕ ТЕПЛОВЫХ ВОЛН В НЕЛИНЕЙНОМ АНИЗОТРОПНОМ ПРОСТРАНСТВЕ
В. Ф. Формалёв, И. А. Селин, Е. Л. Кузнецова
Московский авиационный институт (государственный технический университет),
125993, Москва, А-80, ГСП-3, Волоколамское шоссе, 4.
E-mails: [email protected], [email protected]
Получено аналитическое решение краевой задачи с граничным условием первого рода в начале координат для квазилинейного уравнения параболического типа в анизотропном пространстве со степенной зависимостью компонентов тензора теплопроводности от температуры. Анализ решения показал волновой характер распространения тепла в анизотропном пространстве с конечной скоростью в отличие от бесконечной скорости для линейного уравнения параболического типа. Установлено, что фронт тепловой волны в анизотропном пространстве имеет форму эллипсов на плоскости. Выполнен анализ результатов.
Ключевые слова: квазилинейное параболическое уравнение, тепловая волна, тензор теплопроводности, температура, тепловой поток, фронт.
Введение. Известно, что распространение тепла, описываемое линейным уравнением параболического типа, осуществляется с бесконечной скоростью [1—5], что следует из любого фундаментального решения линейного уравнения диффузии. Этот факт приводит к парадоксу распространения тепла, поскольку все опытные данные показывают хотя и большую, но конечную скорость распространения тепла. Для преодоления этого парадокса ещё Максвелл выдвинул гипотезу о том, что распространение тепла имеет не только диффузионный, но и волновой характер.
В данной статье рассматривается краевая задача для квазилинейного уравнения теплопроводности в анизотропном пространстве и устанавливается волновой характер распространения тепла, то есть в каждый момент времени изотермическая поверхность разделяет конечное пространство с ненулевой температурой и полубес-конечное пространство с нулевой (начальной) температурой.
1. Постановка задачи. Рассматривается нестационарное распределение температуры Т(x,y,t) в двумерном анизотропном пространстве при действии источника температуры То в начале координат 0(0, 0), то есть рассматривается следующая начально-краевая задача при (x, у) G R2, t > 0:
Владимир Федорович Формалёв (д.ф.-м.н., проф.), профессор, каф. вычислительной математики и программирования. Илья Александрович Селин, аспирант, каф. вычислительной математики и программирования. Екатерина Львовна Кузнецова (к.ф.-м.н., доц.), докторант, каф. вычислительной математики и программирования.
(1)
Т (0,0,t) = Т0, Т(x, у, 0) = 0.
(2)
(3)
Компоненты тензора теплопроводности определяются выражениями [6]
Ахх (Т) = А? (Т) cos2 у + An (Т) sin2 у>, Ayy (Т) = А? (Т) sin2 у + А^ (Т) cos2 у>,
Axy(Т) = Аух(Т) = (А£(Т) — An, (Т)) sin V cos ^,
где А^(Т) = к^Тст, Ап(Т) = knТст; к^, kn — const; V — угол, ориентирующий главную ось 0£ и ось On относительно декартовой оси Ox.
Таким образом, необходимо найти нестационарное распределение температур Т(x,y,t) под действием точечного источника (2), приложенного в точке О.
2. Метод решения. Поворотом вокруг начала координат на угол V декартовой системы координат x = £ cos v — n sin V, У = С sin V + П cos v задача (1)-(3) сводится к задаче для уравнения, не содержащего смешанных дифференциальных операторов при (£, n) G R2, t > 0:
<9£ У + a??
дп / ’
Т (0,0,і) = То, Т (£,п,0)=0.
Переходя к новой системе координат ж і = £(Ь/к^)1/2, Х2 = п(Ь/^п )1/2, где Ь — любое число (например, Ь =1), получаем
д
дТ
д
дТ \
St <9жі V дх\) дх2 V дх2) ’ Т (0,0,t)= То, Т (xl,x2,0) = 0.
(4)
(5)
Будем искать решение задачи (4), (5) в автомодельном виде с использованием подстановки
Т (Ж1,Ж2,£)= £“0(^1,92), (6)
где
ql = xi
/te,
q2 = x2
/te.
(Т)
В (б), (Т) а и в подлежат определению. Для их определения подставим (б), (Т) в (4), (5), получим
л/ \ of дУ дУ \
а%Ъ<й) - p\qi^-----h <?2- 1
д92Л
t“(CT+1)-2e а
д91 V дд^ дд2 V дд2'
дУ \
д
дУ
t“0(0,0)= То, t“0(ql,q2)=0.
Из (В), (9) получаем
а = 0, в =1/2.
С учётом (10) задача (В), (9) запишется при (q1, q1) Є R2 в виде
(В)
(9)
(10)
д / дУ \ д / дУ \ 1 / дУ дУ \
-— б»*7-— 4- а-— 0а-— + - яі-------І-ач-— =
дУ
+ а
д
дУ
+
дУ
V" ' “9q2 V" ' 2 V?1 dqi У(0,0) = То.
+ q2
дУ
(11)
(12)
Пусть функция 0(^1, 92) является радиально симметричной, т. е. 0(^1, ?2)=0(?1(г),
?2(г)) = 0(г), тогда задача (11), (12) трансформируется в следующую задачу Коши для нелинейного обыкновенного дифференциального уравнения:
d f ~ dy л dy
2а^г ^)+г^
(13)
а— 1
t
0
а
0
0(0) = Т0, Г 0'(0) = 0. Первый интеграл уравнения (13) записывается в виде
0<т0' + — е = с1г
2а
причём в силу симметрии при г = 0 постоянная Сі = 0. Следовательно,
_ ~ г ~
0<т0' + —0 = 0;
2а
0(0) = То.
(14)
(15)
Уравнение (14) является обыкновенным дифференциальным уравнением с разделяющимися переменными. Его общее решение имеет вид
/ аг2 \ 1/ст
ад = (<?,--) ,
в котором постоянная интегрирования С2 может быть определена из краевого условия (15):
м = (т»-^
аг2 \ 1/ст
(16)
Возвращаясь к декартовым координатам, из (16) получаем решение исходной задачи (1)-(3):
Т(ж, у,і)
Т£ - -^—((хсовір + увіпір)2^ + (-хвіпір + усовір)2^} 4аі V кс кп)
1/0
(17)
3. Анализ результатов. На рис. 1-3 приведены результаты расчётов температурных полей по формуле (17) со следующими данными: а € {0,5; 1; 1,5}, ^ = 0, к£ = 5, кп = 1, Ь =1, ср = 2000 Дж/м3 • К, ТО = 5.
Из представленных рисунков видно, что областями ненулевого решения в различные моменты времени являются области, при (р = 0 ограниченные эллипсами
с полуосями 1Х = по оси Ох, 1у = по оси Оу. Температурное
возмущение существует только внутри этих эллипсов.
Из рисунков следует, что распределение температур имеет волновой характер с фронтом волны на изотерме Т = 0, описываемой эллипсами с полуосями 1х, 1у, а
Т(х,0,і), К
Т(0,у,Ґ), К
Рис. 1. Распределение температур при а = 0,5: а — сечение у = 0, б — сечение х = 0
а
Т(х,0,$), К
Г(0,у,*),К
Рис. 2. Распределение температур при а = 1: а —сечение у = 0, б — сечение х = 0 Т(х, 0,4), К Г(0,»,4),К
Рис. 3. Распределение температур при а = 1: а —сечение у = 0, б — сечение х = 0
фронт тепловой волны распространяется со скоростями Ух = вдоль оси Ох,
уу = \!то^Тл ВД°ЛЬ оси °У-
Компоненты вектора теплового потока вычисляются по формулам
_ дТ дТ _ дТ дТ
Чх~ ххд^ + Хху~д^' ухд^+ ууд^'
В соответствии с (17) компонента теплового потока дх по оси Ож будет определяться выражением
</ж(ж,0,і) = Ххх— (Тп--------ж2^) ^ (-------жУ
4 У ’ сЛ 0 4) V 2а% )
При 0 < а < 1 тепловой поток на границе фронта равен нулю. При а > 1 тепловой поток на границе фронта не существует. При увеличении параметра а поверхности, образуемые изотермами, не равными нулю, имеют большую область определения.
В точках фронта при а ^ 1 первые производные по времени и первые и вторые производные по пространственным переменным не существуют.
При 0,5 ^ а < 1 первые производные по времени и по пространству существуют внутри областей, ограниченных эллипсами с полуосями 1х, 1у, ив точках фронта. При 0,5 < а < 1 вторые производные в точках фронта не существуют.
При 0 < а < 0,5 существуют производные по времени и вторые производные по пространственным переменным внутри областей, ограниченных эллипсами с полуосями 1х, 1у, и в точках фронта.
а
а
Выводы. Впервые получено и проанализировано аналитическое решение краевой задачи в нелинейном анизотропном пространстве для квазилинейного уравнения теплопроводности с тензором теплопроводности, компоненты которого являются степенными функциями температуры. Это решение показало волновой характер распространения тепла. Волновыми фронтами в анизотропном случае в различные моменты времени являются эллипсы на плоскости.
Проведён анализ существования решения на фронте волны и производных разных порядков по времени и пространству на всей плоскости в зависимости от показателя степени температуры, от которой зависят компоненты тензора теплопроводности.
Работа выполнена при поддержке РФФИ (проекты 08-08-00880-а, 08-08-00896-а) и Президента Российской Федерации (проекты МК-646.2008.8, МК-1184.2009.8, МК-1669. 2009.8).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кудинов В. А. Анализ нелинейной теплопроводности на основе определения фронта температурного возмущения // Теплофизика высоких температур, 2006. — Т. 44, №5. — C. 577-585.
2. Лыков А. В. Теория теплопроводности. — М.: Высш. шк., 1967. — 600 с.
3. Шашков А. Г., Бубнов В. А., Яновский С. Ю. Волновые явления теплопроводности. — М.: Эдиториал УРСС, 2004. — 290 с.
4. Самарский А. А., Галактионов В. А., Курдюмов С. П., Михайлов А. П. Режимы с обострением в задачах для квазилинейных параболических уравнений. — М.: Наука, 1987. — 478 с.
5. Соболев С. Л. Процессы переноса и бегущие волны в локально-неравновесных системах // УФН, 1991. — Т. 161, №3. — C. 5-29.
6. Формалёв В. Ф. Тепломассоперенос в анизотропных телах. Обзор // Теплофизика высоких температур, 2001. — Т. 39, №5. — C. 810-832.
Поступила в редакцию 11/X/2009; в окончательном варианте — 15/III/2010.
MSC: 80A17, 80M25
SIMULATION OF HEAT WAVES IN AN NONLINEAR ANISOTROPIC SPACE
V. F. Formalyov, I. A. Selin, K. L. Kuznetcova
Moscow Aviation Institute (State University of Aerospace Technologies),
4, Volokolamskoe Shosse, Moscow, 125993.
E-mails: [email protected], [email protected]
For the first time analytical solution of the problem with boundary conditions in the non-linear anisotropic space for the quasilinear parabolic heat equation where heat conductivity tensor’s components are temperature functions is obtained. This solution has shown a wave type of the heat diffusion. Wave’s fronts are ellipse in the anisotropic space at different time moments. Existence of the solution in wave’s front and time derivative and space derivatives different orders has been analysed dependence from power of temperature which influences heat conductivity tensor’s components.
Key words: quasilinear parabolic equation, heat wave, heat conductivity tensor, temperature, heat flow, front.
Original article submitted 11/X/2009; revision submitted 15/III/2010.
Vladimir F. Formalyov (Dr. Sci. (Phys. & Math.)), Professor, Dept. of Computational Mathematics & Programming. Iliya A. Selin, Postgraduate Student, Dept. of Computational Mathematics & Programming. Ekaterina L. Kuznetcova (Ph. D. (Phys. & Math.)), Doctoral Candidate, Dept. of Computational Mathematics & Programming.
