CC BY
14
УДК 66.01.011
Б01: 10.17277/уе81тк.2016.04.рр.596-601
СПОСОБ ПОВЫШЕНИЯ КАЧЕСТВА МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ НЕСТАЦИОНАРНЫХ ТЕМПЕРАТУРНЫХ ПОЛЕЙ
Е. Н. Туголуков1, А. В. Непрокин2, А. В. Горбунов2, В. М. Нечаев3
Кафедры: «Техника и технологии производства нанопродуктов» (1); tugolukov.en@mail.ru; «Биомедицинская техника» (2);
«Технологические процессы, аппараты и техносферная безопасность» (3), ФГБОУ ВО «ТГТУ», г. Тамбов, Россия
Ключевые слова: аналитическое решение; математическое моделирование; метод конечных интегральных преобразований; сходимость рядов; температурные поля.
Аннотация: На примере задачи нестационарной теплопроводности для многослойной неограниченной пластины с произвольными начальными и неоднородными граничными условиями показан способ улучшения результатов компьютерной реализации аналитических решений задач математической физики. Представлен алгоритм использования метода конечных интегральных преобразований по координате, вдоль которой теплофизические характеристики рассматриваемого объекта меняются ступенчато. Проведено сравнение видов решений, полученных с выделением стационарной составляющей и без него.
Применение современных средств компьютерной техники позволяет использовать задачи математической физики в частных производных для решения широкого круга инженерных и научно-исследовательских задач. Достоинство аналитических подходов при математическом моделировании сложных процессов переноса, помимо отсутствия систематических погрешностей методов решений, -относительно малый объем вычислений и, как следствие, малое время получения результатов расчетов, что позволяет использовать аналитические подходы в системах реального времени. Для решения рассматриваемого класса задач наиболее целесообразным представляется метод конечных интегральных преобразований, теория которого разработана Н. С. Кошляковым [1].
Компьютерная реализация аналитических решений задач математической физики в частных производных имеет свои особенности, в частности, связанные с необходимостью суммирования рядов, входящих в решения. Мероприятия, направленные на повышение сходимости рядов, позволяют повысить качество результатов расчетных работ.
Рассмотрим задачу теплопроводности для п-слойной неограниченной пластины с произвольным начальным и неоднородными граничными условиями [2]:
д ti (xi, т) 2 д2ti (xi, т) 1У 1 = а~—^ 2 , i = 1,2,..., N 0 < xi <Я1, т> 0; (1)
St г 5х?
t, (х, ,0) = f (x); (2)
д t1 (0, т) / , . ч
_ai ((0,T)-tcl) = 0;
ХN дtN(RN,т) +а^ ( (rn,T)-tcN ) = 0;
д rN v '
t ( т) = t ( т). х дtJ (,т) =х дtJ +l (0т)
tj ((,т)=tj+i((,т); х j —д7— " х j+1 д r. . '
(4)
] = 1,2,..., Ж -1, (5)
Ч -гз+1
где tj (х{ ,т) - температурное поле 1-й области; х - пространственная координата; т - время; , а2 - соответственно коэффициенты теплопроводности и температуропроводности /-й области; N - число слоев многослойной области; - координата границы /-й области; а^ аж - коэффициенты конвективной теплоотдачи от внешних поверхностей в окружающую среду; tc-í, ^ - температуры окружающей среды на внешних границах области.
Для исключения координаты х, вдоль которой свойства пластины изменяются ступенчато, используем формулу перехода к изображениям
N ^ Ят
и (Ц Т)=Х I ^ (хт, т)т (хт, Ц)Схт. (6)
т = 1 ат о
Ядро интегрального преобразования Wm (хт, ц) - решение задачи Штурма-Лиувилля с соответствующими однородными граничными условиями, определяемое с точностью до постоянного множителя (здесь ц - параметр):
С! ^ ((т, ц) + Wm (хт, ц) = 0, т = 1,2,..., Ж, 0 < хт < Ят; (7)
Сг а
Iл Iт "т
dW1 (Л0, ц) Xi--a1 W1 (R0, ц) = 0;
d r
-N'
dWN ( Rn , ц)
d r
N
+aN Wn (Rn , ц)= 0;
/ \ / ч dW. (0, ц) dW,+1 (0, ц)
WJ (j, ц) = WJ+1(0, ц);
Общее решение уравнения (7) имеет вид
Wm (xm, ц) = Ст sin
ц xm
+ Фт
(8) (9)
j = 1,2,..., N -1. (10) (11)
Из граничных условий (8) - (10) определим числа фт, Ст, ц, причем
С = 1:
Ф j+1 = arctg
Ф1 =- arctg
(^,+1 a, ( ц Л.
Х1 ц |.
a1 a1
V+1" J
vX jaj+1
tg
a.
V J
- + Ф j
W
j j
j = 1, 2,..., N -1;
81И
CJ +1 = CJJ
ц R,
aj
V_J
-+Ф,
sin (ф j+1)
(12)
(13)
(14)
ц - последовательные положительные корни уравнения
^ cos (^Л
aN I aN
(
+ aN sin
— rn
+ ФN
= 0.
(15)
Для перехода к изображениям необходимо формулу (6) применить почленно к уравнению (1) и начальному условию (2). В изображениях получаем задачу:
+ ц2 и (ц, т)=а^ (К, ц)^ ^ ((?!, ц )tc!; (16)
X
N — Rm
U (—,0)=Z—Г i fm (xm )Wm (xm ,—)dxm.
a
(17)
т = 1 т р ,
Решение задачи (16) - (17) в изображениях имеет вид
U(—, т) = exp (——2т) U(—,0) +—2—(anW R, —)tCN + aW (R, —) C) (exp (—2т) -1)
I — —
Обратный переход для задачи (1) - (5) осуществляется по формуле
и (ц т)^т (т, ц)
1 2 ' п=1
tm (rm,т) £"
(18) (19)
где
N Л Rm — „
N
Z = £ Ol J Wm2 (xm, —)dxm = 0,5 £ ^
m = 1 am 0 m = 1 am
r —
( ( sin
-
- Rm
V V am
- + Ф2
A (
008
— Rm
\
- + Ф2
— sin (Ф2 )COs (Ф2 )
(20)
J J
Суммирование в выражении (19) проводится по значениям собственных чисел д. Компьютерная реализация решения (19) задачи (1) - (5) показывает,
что наличие слагаемого
—2— (a nW (Rn , —)tcN +a1 W (R,, —)t^)) exp (—2 т) —1)
в формуле (18) заметно ухудшает сходимость ряда в (19).
Улучшить сходимость ряда в решении задачи (1) - (5) и, как следствие, повысить качество компьютерной реализации решения можно, предварительно выделив стационарную составляющую решения
^ (х1, т) = ^ (х1) + Р1 (х,, т), 1 = 1,2,..., N, (21)
где (х1) - решение стационарной задачи с неоднородными граничными усло-
виями:
d 2 S, ( x,)
-= 0, i = 1,2,..., N, 0 < x, < R,;
d xj
—1 ^ _a1 (S1 (0) —td )=0;
dx1 4 '
— N-N ( N ) +aN (SN (RN ) — tcN ) = 0;
(22)
(23)
(24)
, (,) = ,„ (0); ,
тогда решение стационарной задачи (22) - (25) примет вид:
$ (х,) = Ai х, + Bi;
Ai = (
tcN - tc1
1 ^ Лк 1
— + ^+-
a1 к=1 Х к a N (tcN -tc1 )fR + £ ^ + aL
j = 1,2,..., N -1, (25)
(26) (27)
Bi = ^N
VXi к=i+1 Хк a N
N
1 Rk 1
—+ £+—
a1 к=1 Хк a N
(28)
Р (х,, т) - решение нестационарной задачи с однородными граничными условиями:
5 Р (х,, т)= 2 52Р (х,, т)
дт
д x2
i = 1,2,..., N, 0 <xi <Ri, т>0;
Pi (хг ,0)=f (x)-S, (x);
Х
д P1 (0, т)
д x1
-a1 P1 (0, т) = 0;
Хn 5Rn,т) +aNpv (RN,т) = 0;
д x
N
PJ (RJ, т) = PJ+1 (00, т); хJ
дpj (Rj,т) =х ^j+1 О0,т)
] = 1,2,..., N -1. (33)
" ЛУ+1
Решение задачи (29) - (33) может быть получено методом конечных интегральных преобразований с использованием формулы перехода к изображениям (6) относительно функции р (х,, т)
N л Ят
и (ц,т)=Е"Т I Рт (хт,т)Wm (хт, ц)Схт. (34)
т=1ат 0
Обратный переход для задачи (29) - (33) осуществляется аналогично формуле (19)
д х,
J +1"
д х,
P (х,,т)=£
Ü1 (ц, т) Wi (х, , ц)
n = 1
где Ü1 (ц, т) - решение задачи в изображениях:
+ ц2и1 (ц, т) =0;
d т
N х Л
Ü1 (ц,0)=! ~mt j (fm (хт ) - Sm (хт )) Wm (хт, ц) <К.
т = 1 am 0
(35)
(36)
(37)
Решением задачи (36) - (37) является функция
U1 (—, t) = U1 (—,0 )exp (——2т)
Таким образом, решение исходной задачи (1) - (5), полученное при выделении стационарной составляющей, имеет вид
Результатом выделения стационарной составляющей при решении задачи (1) - (5) является то, что ряд в решении (39) имеет существенно лучшую сходимость. Благодаря этому повышается качество компьютерной реализации решения исходной задачи теплопроводности, так как снижаются объем вычислений и общая вычислительная погрешность. Кроме того, процедура получения решения Р1 (х1 , т) нестационарной задачи теплопроводности значительно упрощается, а решение (х1 ) стационарной задачи (22) - (25) может иметь самостоятельную ценность при расчетах стационарных процессов.
Аналогичный подход при получении аналитических решений других задач нестационарной теплопроводности в ходе разработки математических моделей температурных полей использован в ряде работ [3, 4].
Список литературы
1. Кошляков, Н. С. Уравнения в частных производных математической физики / Н. С. Кошляков, Э. Б. Глинер, М. М.Смирнов. - М. : Высшая школа, 1970. -
2. Туголуков, Е. Н. Математическое моделирование технологического оборудования многоассортиментных химических производств : монография / Е. Н. Туго-луков. - М. : Машиностроение, 2004. - 100 с.
3. Коновалов, В. И. О возможностях использования точных, интервальных и приближенных аналитических методов в задачах тепло- и массопереноса в твердых телах. Часть 1: Постановка проблемы. Точные аналитические методы / В. И. Коновалов, Е. Н. Туголуков, Н. Ц. Гатапова // Вестн. Тамб. гос. техн. ун-та. -1995. - Т. 1, № 1-2. - С. 75 - 90.
4. Туголуков, Е. Н. Моделирование теплопередачи в биметаллических аппаратах с каналами охлаждения в стенках / Е. Н. Туголуков, В. А. Богуш, А. Г. Ткачев // Вестн. Тамб. гос. тнхн. ун-та. - 2003. - Т. 9, № 1. - С. 42 - 49.
Method for Improving the Quality of Mathematical Modeling of Unsteady Temperature Fields
E.N. Tugolukov1, A.V. Neprokin2, A.V. Gorbunov2, V. M. Nechaev3
Departments: "Technology and Manufacturing ofNanoproducts" (1); tugolukov.en@mail.ru; "BiomedicalEngineering" (2); "TechnologicalProcesses, Devices and Technosphere Safety" (3), TSTU, Tambov, Russia
Keywords: analytical solution; mathematical modeling; method of finite integral transformations; convergence of series; temperature fields.
Abstract: Using the example of transient heat conduction problem for the multilayer infinite plate with arbitrary initial and inhomogeneous boundary conditions, we showed how to improve the results of the computer implementation of analytical solutions of mathematical physics problems with inhomogeneous boundary conditions. The algorithm for the method of integral transformations in the coordinate along which the thermal characteristics of the object change in steps. A comparison of the types of solutions obtained without isolation and separation of the stationary component is made.
n = 1
(39)
712 с.
References
1. Koshlyakov N.S., Gliner E.B., Smirnov M.M. Uravneniya v chastnykh proizvodnykh matematicheskoi fiziki [Partial differential equations of mathematical physics], Moscow: Vysshaya shkola, 1970, 712 p. (In Russ.)
2. Tugolukov E.N. Matematicheskoe modelirovanie tekhnologicheskogo oboru-dovaniya mnogoassortimentnykh khimicheskikh proizvodstv [Mathematical modeling of technological equipment of chemical manufactures mnogoassortimentnyh], Moscow: Mashinostroenie, 2004, 100 p. (In Russ.)
3. Konovalov V.I., Tugolukov E.N., Gatapova N.Z. [Possibilities for the Use of Some Exact, Interval and Approximate Analytical Methods for Heat and Mass Transfer Problems in Solids. Part 1: Statement of Problem. The Exact Analytical Methods], Transactions of the Tambov State Technical University, 1995, vol. 1, no. 1-2, pp. 75-90. (In Russ., abstract in Eng.)
4. Tugolukov E.N., Bogush V.A., Tkachyov A.G. [Modeling of Heat Transfer in Bimetal Apparatuses with Canals of Cooling in Walls], Transactions of the Tambov State Technical University, 2003, vol. 9, no. 1, pp. 42-49. (In Russ., abstract in Eng.)
Weise der Erhöhung der Qualität der mathematischen Modellierung der nichtstationären Temperaturfelder
Zusammenfassung: Am Beispiel der Aufgabe der nichtstationären Wärmeleitfähigkeit für die mehrschichtigen unbeschränkten Platte mit den willkürlichen Anfangsbedingungen und mit den verschiedenartigen Randbedingungen ist die Weise der Verbesserung der Ergebnisse der Computerrealisierung der analytischen Lösungen der Aufgaben der mathematischen Physik mit den ungleichartigen Randbedingungen vorgeführt. Es ist der Algorithmus der Nutzung der Methode der endlichen Integralumgestaltungen nach der Koordinate, die entlang die wärme-physikalischen Charakteristiken des betrachtenden Objektes sich stufig ändern, angeführt. Es wird der Vergleich der Arten der Lösungen, die ohne Absonderung und mit der Absonderung der stationären Komponente erhalten sind, durchgeführt.
Moyen de l'augmentation de la qualité de la modélisation mathématique des champs de température non stationnaires
Résumé: A l'exemple de la tâche non-stationnaire de la conductivité thermique pour une plaque multicouche illimitée avec les conditions limites des initiales facultatives et hétérogènes est montré le moyen de l'amélioration des résultats de la réalisation informatique des solutions analytiques des problèmes de la physique mathématique avec les conditions limites hétérogènes. Est présenté l'algorithme de l'utilisation de la méthode des transformations finales intégrales sur une coordonnée le long de laquelle les caractéristiques thermophysiques de l'objet examiné changent en étapes. Est effectuée une comparaison des solutions obtenues.
Авторы: Туголуков Евгений Николаевич - доктор технических наук, профессор кафедры «Техника и технологии производства нанопродуктов»; Непрокин Алексей Владимирович - студент; Горбунов Алексей Викторович - доктор медицинских наук, доцент кафедры «Биомедицинская техника»; Нечаев Василий Михайлович - кандидат технических наук, доцент кафедры «Технологические процессы, аппараты и техносферная безопасность», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
Рецензент: Гатапова Наталья Цибиковна - доктор технических наук, профессор, заведующая кафедрой «Технологические процессы, аппараты и техносферная безопасность», ФГБОУ ВО «ТГТУ», г. Тамбов, Россия.
