Эволюция пространственно локализованных тепловых возмущений

DOI: 10.18698/1812-3368-2015-6-16-24 УДК 517.958

ЭВОЛЮЦИЯ ПРОСТРАНСТВЕННО ЛОКАЛИЗОВАННЫХ ТЕПЛОВЫХ ВОЗМУЩЕНИЙ

Л.К. Мартинсон, О.Ю. Чигирёва

МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: mathmod@bmstu.ru

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

Ключевые слова: краевая задача нелинейной теплопроводности, эффект пространственной локализации тепловых возмущений.

EVOLUTION OF SPATIALLY LOCALIZED THERMAL PERTURBATIONS L.K. Martinson, O.Yu. Chigireva

Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: mathmod@bmstu.ru

The article considers a one-dimensional boundary problem of the nonlinear heat conduction equation in a two-dimensional layer filled with a medium with bulk heat absorption. Numerical solutions to this problem for different characteristic values, while using a standard difference scheme, confirm the theoretical conclusions about the propagation mode of thermal perturbations. Thermal perturbations from the heated walls prove to propagate through the nonlinear medium with bulk heat absorption at the finite rate of the front propagation. An effect of spatial localization of thermal perturbations, which can even reach the finite depth for nonterminating time, is also observed at the defined values of the problem characteristics.

Keywords: boundary problem of nonlinear heat conduction, effect of spatial localization of thermal perturbations.

Введение. Изучение особенностей задач для нелинейных дифференциальных уравнений параболического типа [1-3] показало, что нелинейности не только изменяют количественные характеристики процессов, описываемых такими уравнениями, но и качественно изменяют характер их протекания. В нелинейных процессах теплопроводности могут наблюдаться эффекты конечной скорости распространения и пространственной локализации тепловых возмущений [4, 5]. Именно эти нелинейные эффекты обусловливают необычный ход эволюции температурных полей в задачах, рассматриваемых ниже.

Физическая постановка задачи и математическая модель процесса. Рассмотрим задачу о нахождении нестационарного температурного поля и (х,£) в плоском слое конечной толщиной /, заполненном средой, коэффициент теплопроводности Л которой изменяется в зависимости от температуры и по степенному закону Л (и) = Л0иа, где Ло > 0; а > 0 — параметр нелинейной среды.

В начальный момент времени £ = 0 температура слоя равна нулю. При £ > 0 на поверхности х = 0 этого слоя поддерживается температура и0, а на поверхности х = I — нулевая температура. Во внутренних точках слоя происходит поглощение теплоты, удельная мощность которой f (и) = р0иа, где а > 0; р0 > 0 — коэффициент поглощения.

Математическая модель процесса имеет вид

дм д

pCdt = Ло дХ

u

.дм дх

- р0ма, 0 < х < /, t> 0;

и (х, 0) = 0, 0 < х < /; и (0,£) = Ц, и (М) = 0, £ > 0.

Здесь р, с — плотность и удельная теплоемкость среды.

Согласно общей теории нелинейных тепловых процессов, при а > 0 тепловое возмущение от нагретой поверхности х = 0 распространяется в виде тепловой волны с конечной скоростью перемещения ее фронта х* (£). Кроме того, при некоторых значениях параметра а в задаче (1) также наблюдается эффект пространственной локализации, когда тепловое возмущение от нагретой поверхности проникает в слой на конечную глубину Ь и не достигает поверхности х = I даже при £ ^ ж.

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

Ust (х) =

U0 (1 - L

0, х > L.

<г+1-а

х < L;

Здесь

L =

2U

a+1-a 2

0

а

(а + 1 + а)! 2

р (а + 1 — а)2

= Ло/ (pc), p = ро/ (pc).

Отметим, что формула (3) имеет смысл при а < а + 1. Форма профиля стационарного решения зависит от соотношения параметров а и а. Качественный вид стационарных температурных профилей (2) в слое при заданном значении а и различных значениях а представлен на рис. 1.

2

2

а

Далее будем исследовать задачи, для которых Ь < /, когда в (1) наблюдается эффект пространственной локализации тепловых возмущений.

С учетом теорем сравнения, приведенных в работе [6], существование стационарного решения

(2) означает следующее: для любо-Рис. 1. Качественный ввд стационар- го £ е [0, решение задачи (1) ных температурных профилей в слое

при а = 5/2 и значениях а = 1/2 (1), мажорируется стаци°нарньш ре-3/2 (2), 5/2 (5) шением, т.е. и (х, £) < и^ (х). Это

позволяет записать приближенное решение задачи (1) в форме тепловой волны с конечной скоростью перемещения ее фронта

u (x, t) =

U 1-

x

x*

a+1-а

, 0 < x < x* (t);

0, x* (t) < x < l.

Для нахождения закона движения фронта тепловой волны запишем интегральное условие теплового баланса

x*

=a4 t

0 0

.du

u

x*

_ ¿х — р иа (х,£) ¿х. дх У

о

Подставим предполагаемую форму решения (4) в соотношение (5) и вычислим интегралы в обеих частях равенства (5). Учитывая, что на фронте х = х* (£) тепловой поток равен нулю, приходим к следующему дифференциальному уравнению:

d

dt

а + 3 — а

= b0 - bix* (t),

_1 а + 3 — а

где Ъо = 4Ща2^—^-^; Ъг = р: +

(а + 1 — а)2 а + 1 + а

Очевидно, что искомая функция должна удовлетворять начальному

условию

х*

Интегрируя дифференциальное уравнение (6) и учитывая начальное условие (7), находим

x* (t) = Уb0 (1 — exp(—bit)).

2

Принимая во внимание (3), вычисляем отношение Ь0 /&1, в результате получаем

Рис.2. Характер изменения положения фронта тепловой волны

х*

Переходя в равенстве (8) к пределу при £ ^ ж, видим, что х* стремится к конечному значению Ь при а < а + 3. Очевидно, что этому неравенству удовлетворяют такие значения параметров а и а, при которых имеет смысл формула (3).

Решение (4) описывает тепловое возмущение, распространяющееся в виде тепловой волны, фронт х* (£) которой перемещается с конечной скоростью. Изменение положения фронта тепловой волны подчиняется закону (8) (рис. 2). Следует отметить, что согласно (8) характерное время установления стационарного распределения температуры в слое можно оценить как £уст = Ь-1.

Разностная схема. Рассмотрим задачу (1) в постановке, когда начальное распределение температуры в слое задает функция ^ (х): ^ (х) > 0, ^ (х) е С [0, /], таким образом

k (м) дм

— f (u), 0 < х < l, t> 0;

( дм д д£ дх

u (х, 0) = ^ (х), 0 < х < l;

, u (0, t) = Uo, u (l,t) = 0, t > 0,

где k (u) = а2ua, f (u) = pua.

Для численного решения задачи (9) применим разностный метод. В области D = [0 < х < l] х [0 < t < T] зададим равномерную сетку

uhT = |(хп, tm): хп = nh, tm = тт, n = 0, N, m = 0,m| с шагами h = l/N и т = T/M по переменным х и t.

Обозначим Уп = У (хп ,tm), Уп = У (хп,tm+l), ^п = ^ (хп),

fn = f (Уп), = k ^Уп-1 Уп^, используя четырехточечный шаблон (рис. 3), запишем неявную разностную схему [7, 8]:

Уп — Уп

1

h2

&п+1 (Уп+1 - Уп) - (Уп - Уп-1) - Лп, п = 1, N - 1.

Из граничных условий определяем у0 = Ц0, у^ = 0.

При решении системы

Рис. 3. Четырехточечный шаблон разностной схемы

монотонна, имеет погрешность аппроксимации О (Н2 + т) и нелинейна относительно уп.

Для отыскания решения полученной системы нелинейных уравнений (10) на временном слое £т+1 применим метод последовательных приближений

VnJ

Уп

т

1

h2

k

'n+1 l yn+l

уП°

-kns-1)

где s = 1, 2,... — номер итерации.

В качестве начального приближения уП0) примем значение уп с предыдущего временного слоя и вычислим величины кП0) и /П0). Решая систему (11) методом прогонки [9], находим уП^ — первое приближение решения уп и т.д. Условие окончания итерационного процесса

< г, где г — заданное значение абсо-

имеет вид max

1< n<N-1

уП)

-

лютной погрешности.

Выбор величин Н и т осуществляется на основе методики, предложенной в работе [10].

Результаты численных расчетов. Приведем примеры численного расчета температуры в слое при следующих значениях параметров: и0 = 1; I = 1; а2 = 1; а = 5/2. Начальное распределение температуры ^ (х) = 0.

Пример 1. Для значений параметров а = 1/2 и р = 32/9, используя формулу (3), находим Ь = 1/2. Следовательно, стационарное распределение температуры в слое имеет вид

2 1 (1 - 2x)3 , 0 < x< -;

uSt) (x) ^ 1 2

0, - < x < 1. 2 < <

Распределения температуры в фиксированные моменты времени ¿1 = 0,01 с, Ь2 = 0,1 си Ь3 = 0,5 с представлены на рис.4, а. Расчеты показали, что при Ь > 1,0 с процесс выходит на стационарный режим (12).

Пример 2. При значениях а = 3/2 и р =10 получаем Ь = 1/2. В этом случае стационарное распределение температуры представляет

Рис. 4. Распределения температуры в слое в фиксированные моменты времени = 0,01 с (1), ¿2 = 0,1 с (2) и = 0,5 с (3) и стационарные профили (ж) и

<42*

собой "линейный" профиль

1 - 2х, 0 < х < -; и2}

0, - < х < 1. 2 < <

Распределения температуры в моменты времени ¿1 = 0,01 с, ¿2 = 0,1 си = 0,5 с приведены на рис.4, б. При £ > 1,2 с процесс выходит на стационарный режим (13).

В рассмотренных примерах тепловое возмущение, распространяющееся от левой поверхности х = 0, проникает в слой на конечную глубину Ь = 1/2 (за пределами этой глубины температура равна нулю) и не достигает правой поверхности х = 1 даже при £ ^ ж.

Устойчивость стационарного решения. Для значений параметров задачи, рассмотренных в примере 1, выясним, как будет вести себя решение краевой задачи (9), если начальное распределение температуры в слое имеет вид

2

V (х) = < V1 - ' 0 < х<Ь»;

0, Ь0 < х < I,

где Ь0 > Ь. В рассматриваемом примере полагаем Ь0 = 3/4.

Расчеты показали (рис. 5), что в этом случае фронт х* (£) тепловой волны движется от точки х* (0) = Ь0 в противоположном направлении, и решение и (х,£) также стремится к стационарному решению и^} (х) в процессе эволюции теплового возмущения. Этот важный результат доказывает, что стационарный профиль (2) является устой-

и

Рис. 5. Эволюция теплового возмущения в слое, когда ф

чивым предельным распределением температуры в слое при любых начальных условиях.

Выводы. Численное решение краевой задачи (1) подтверждает теоретический вывод о том, что при некоторых значениях параметров задачи в процессе эволюции тепловых возмущений от нагретой стенки наблюдаются нелинейные эффекты конечной скорости распространения фронта тепловой волны и пространственной локализации тепловых возмущений.

ЛИТЕРАТУРА

1. Мартинсон Л.К., Малое Ю.И. Дифференциальные уравнения математической физики. М.: Изд-во МГТУ им. Н.Э. Баумана, 2002. 368 с.

2. Мартинсон Л.К. Исследование математической модели процесса нелинейной теплопроводности в средах с объемным поглощением. В кн.: Математическое моделирование. Процессы в нелинейных средах. М.: Наука, 1986. С. 279-309.

3. Маслое В.П., Данилов В.Г., Волосое К.А. Математическое моделирование процессов тепломассопереноса. М.: Наука, 1987. 362 с.

4. Мартинсон Л.К., Чигирёеа О.Ю. Пространственная локализация тепловых возмущений в нелинейном процессе теплопроводности // Вестник МГТУ им.Н.Э.Баумана. Сер. Естественные науки. 2013. № 4. С. 27-33.

5. Мартинсон Л.К., Чигирёеа О.Ю. Краевые задачи для квазилинейных уравнений параболического типа // Необратимые процессы в природе и технике: Труды Седьмой Всероссийской конференции. В 3 ч. М., 2013. Ч. II. С. 32-33.

6. Самарский А.А., Галактионое В.А., Курдюмое С.П., Михайлое А.П. Режимы с обострением в задачах для квазилинейных параболических уравнений. М.: Наука, 1987. 480 с.

7. Тихоное А.Н., Самарский А.А. Уравнения математической физики. М.: Изд-во МГУ, 2004. 798 с.

8. Самарский А.А. Теория разностных схем. М.: Наука, 1977. 656 с.

9. Амосое А.А., Дубинский Ю.А., Копченоеа Н.В. Вычислительные методы для инженеров. М.: Высш. шк., 1994. 544 с.

10. Матус П.П. О корректности разностных схем для полулинейного параболического уравнения с обобщенными решениями // Журнал вычислительной математики и математической физики. 2010. Т. 50. № 12. С. 2155-2175.

REFERENCES

[1] Martinson L.K., Malov Yu.I. Differentsial'nye uravneniya matematicheskoy fiziki [Differential equations of mathematical physics]. Moscow, MGTU im. N.E. Baumana Publ., 2002. 368 p.

[2] Martinson L.K. Issledovanie matematicheskoy modeli protsessa nelineynoy teploprovodnosti v sredakh s ob'emnym pogloshcheniem. V kn.: Matematicheskoe modelirovanie. Protsessy v nelineynykh sredakh [Research into the process mathematical model of nonlinear thermal conductivity in spatial absorption media. In the book: Mathematic Simulation. Processes in Nonlinear Media]. Moscow, Nauka Publ., 1986, pp. 279-309.

[3] Maslov V.P., Danilov V.G., Volosov K.A. Matematicheskoe modelirovanie protsessov teplomassoperenosa [Mathematic simulation of heat and mass transfer processes]. Moscow, Nauka Publ., 1987. 362 p.

[4] Martinson L.K., Chigireva O.Yu. Spatial Localization of Thermal Perturbations in Nonlinear Process of Heat Conduction. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2013, no. 4, pp. 27-33 (in Russ.).

[5] Martinson L.K., Chigireva O.Yu. Boundary-value problems for quasilinear equations of the parabolic type. Irreversible processes in nature and engineering: Proceedings of the Seventh Russian National Conference. Neobratimye protsessy v prirode i tekhnike: Tr. Sed'moy Vseross. konf. V 3 ch. Moscow, 2013. Part II, pp. 32-33 (in Russ.).

[6] Samarskiy A.A., Galaktionov V.A., Kurdyumov S.P., Mikhaylov A.P. Rezhimy s obostreniem v zadachakh dlya kvazilineynykh parabolicheskikh uravneniy [Blowup regimes in problems for quasilinear parabolic equations]. Moscow, Nauka Publ., 1987. 480 p.

[7] Tikhonov A.N., Samarskiy A.A. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, MGU Publ., 2004. 798 p.

[8] Samarskiy A.A. Teoriya raznostnykh skhem [Difference scheme theory]. Moscow, Nauka Publ., 1977. 656 p.

[9] Amosov A.A., Dubinskiy Yu.A., Kopchenova N.V. Vychislitel'nye metody dlya inzhenerov [Computational approaches for engineers]. Moscow, Vyssh. shk. Publ., 1994. 544 p.

[10] Matus P.P. Correctness of difference schemes for the semilinear parabolic equation with generalized solutions. Zh. Vychisl. Mat. Mat. Fiz. [Comput. Math. Math. Phys.], 2010, vol. 50, no. 12, pp. 2155-2175 (in Russ.).

Статья поступила в редакцию 22.06.2015

Мартинсон Леонид Карлович — д-р физ.-мат. наук, профессор кафедры "Физика" МГТУ им. Н.Э. Баумана.

МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

Martinson L.K. — D.Sc. (Phys.-Math.), Professor, Department of Physics, Bauman Moscow State Technical University.

Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

Чигирева Ольга Юрьевна — канд. физ.-мат. наук, доцент кафедры "Математическое моделирование" МГТУ им. Н.Э. Баумана.

МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

Chigireva O.Yu. — Ph.D. (Phys.-Math.), Associate Professor, Department of Mathematical Simulation, Bauman Moscow State Technical University.

Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.

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

Мартинсон Л.К., Чигирeва О.Ю. Эволюция пространственно локализованных тепловых возмущений // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2015. № 6. C. 16-24.

Please cite this article in English as:

Martinson L.K., Chigireva O.Yu. Evolution of spatially localized thermal perturbations. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2015, no. 6, pp. 16-24.