CC BY
45
УДК 517.94
Н.В. Дилигенский, А.П. Ефимов
ПРИМЕНЕНИЕ ОПЕРАТОРОВ ДРОБНОГО ДИФФЕРЕНЦИРОВАНИЯ ДЛЯ ПОСТРОЕНИЯ АППРОКСИМАТИВНЫХ РЕШЕНИЙ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ1
Статья посвящена проблеме построения приближенных аналитических решений уравнения теплопроводности на основе использования операции дробного дифференцирования.
Ключевые слова: уравнение теплопроводности, аппроксимации, приближенные решения, дробное дифференцирование.
N. V. Diligenskiy, A.P. Efimov
APPLICATION OF FRACTIONAL DIFFERENTIATION OPERATORS FOR CONSTRUCTION OF APPROXIMATE SOLUTIONS OF HEAT TRANSFER EQUATION.
The paper is devoted to the problem of construction of approximate analytical solutions of heat transfer equation based on operations of fractional differentiation.
Keywords: heat transfer equation, approximation, approximate solution, fractional differentiation.
Введение
Эффективность применения аналитических решений нестационарных задач теплопроводности для целей расчета, анализа, управления и оптимизации в определяющей степени зависит от формы получаемых решений. Во многих случаях точные решения оказываются неэффективными в прикладных задачах вследствие сложности представлений, неудовлетворительной сходимости рядов, трудоемкости вычислительных операций. В статье рассматривается подход к построению приближенных аппроксимационных решений параболических задач, опирающийся на использование операторов дробного дифференцирования.
1. Постановка задачи
Сущность подхода изложим на примере построения приближенных решений канонической линейной задачи теплопроводности в пластине конечной толщины с граничными условиями 1 -го рода
дв д 2в
дт dX2
в x=1 = U (т) , (1)
дв dX
в
=о
x=0
= 0 ,
т=0
где в,т, X - соответственно безразмерные температура, время (число Фурье), пространственная координата, и (т) - известная граничная температура.
Решение задачи (1) методом разделения переменных дает представление в виде сходящихся рядов по собственным функциям. Однако такое представление является неэффективным при малых временах вследствие низкой скорости сходимости рядов.
Запишем решение задачи (1) в трансформантах Лапласа в операторной форме [1]
в( X) = W (р, х)й (р), (2)
где
е4рх+Р-^рх
W (р) = —Г---------------------------------------------Г~ (3)
в^р + е V р
1 Работа выполнена в рамках целевой программы «Научные и научно-педагогические кадры инновационной
России на 2009-13 гг.» (государственный контракт №П1448).
оператор перехода (передаточная функция) от U к в, символ " " означает трансформанту Лапласа, р
- оператор преобразования Лапласа
Переход к оригиналам по времени на основе решения в транформантах (2) (3) приводит к представлениям в виде рядов, эффективных при малых временах (больших p ), но неработоспособных
при больших. Они соответствуют внутренним асимптотическим разложениям решения (2) в1, удобных при p — ж .
Построение на основе (2), (3) формальных асимптотических рядов при больших временах (малых p) приводит к сингулярным внешним разложениям в0 , представимым в оригиналах в классах обобщенных функций, использование которых в прикладных решениях затруднительно.
Получаемые таким образом внутренние и внешние разложения не стыкуются друг с другом, не имеют общей области работоспособности, и построение на их основе известными методами сращивания составных асимптотических решений невозможно.
2. Принцип построения аппроксимативных решений
Сформулируем сущность подхода построения аппроксимативных решений задачи (2), (3) в виде следующих положений:
1. Примем в качестве базовых составляющих приближенных решений N, М конечные суммы внутреннего etN и внешнего в0м асимптотических разложений.
2. На основе etN и воМ сконструируем аддитивную конструкцию
в a = (1 -O.K ) • ^t,N + Qk • во,М , (4)
где Qк - гладкая (положим K - дифференцируемая) функция-мультипликатор, обеспечивающая сращивание разнородных разложений.
Мультипликатор Qk (р) обладает следующими предельными свойствами
ltmQ.K (р) = 1, lim Qk (р) = 0 . (5)
р—— 0 р—Ж
Мультипликатор взвешивает ряды вtN и в ОМ с единичным весом в своих - внутренней и внешней -областях асимптотического описания и монотонно уменьшает их вклад до нуля вне своих областей. В соответствии с этим аппроксимативное решение (4) будет иметь внутренние и внешние асимптотики, соответствующие точному решению (2), (3), плавно сопрягающиеся в промежуточной области.
Построим явную форму мультипликатора Qк на основе анализа точного решения задачи (3). Видно, что выражение (3), определяющее сущность и закономерности протекания динамических процессов, зависит только от оператора у[р , отвечающего операции дробного дифференцирования степени Уг [2]. В соответствии с этим положим, что для сохранения структурных свойств решений мультипликатор также целесообразно строить как функцию от у[р. Далее учтем, что решение (3) определяется
экпоненциальной зависимостью от дробного оператора вида exp^^^). На основании этого
мультипликатор будем конструировать как функцию от expY^/^), где у - некоторая константа.
В соответствии с этими двумя вышеизложенными соображениями построим мультипликатор в следующей форме
Q K (4¥) = exp к (Г\[р) • exp(-^/р), (6)
K Z
где expк (z) = многочлен Тейлора степени K для функции exp(z), У - численный параметр,
l=0 1!
определяемый из условий удовлетворения решения (4) граничным условиям с достаточной степенью точности, соответствующей точности использованных при построении решения (4) внутреннего и внешнего разложений. Очевидно, что асимптотическое поведение мультипликатора (6) удовлетворяет требованиям (5), и структура (6) близко соответствует структуре точного решения (3).
За счет выбора количества членов N,М, K соответствующих рядов вtN, вом , Qk и параметра у получаем возможность конструирования аналитических аппроксимативных решений параболических задач теплопроводности с различной степенью точности и описывающих температурные распределения
во всей области определения решений (для те [0,го] , Xе [0,1]) с предельным поведением внутренних и внешних асимптотических разложений.
3. Приближенные решения задачи теплопроводности
Запишем явную форму аппроксимативных решений для оператора (3):
- трёхчленное внутреннее асимптотическое разложение
^,2 (р, х) = ехр(-7Р (1 - х)) + ехр(-7р (1 + х) + ехр(-7Р (3 - х)) +
+0(ехр(-з7Р)) (7)
- двухчленное внешнее асимптотическое разложение
Wо,2( р, X) = 1 + р (х2 +1) + 0( р 2) (8)
Видно, что И^'2( Р, х) и Wо 2( р, х) являются существенно различными зависимостями,
несогласующимися между собой, причем внутреннее разложение (7) имеет так же, как и (3), структуру, соответствующую операциям дробного дифференцирования, а представление (8) - существенно иную, отвечающую операциям дифференцирования в пространствах обобщенных функций.
Используя конструкцию (4) и учитывая различное количество членов во внутреннем в'ы и внешнем
вош асимптотических разложениях, получим на основе (7), (8) следующие аппроксимирующие решения
Шц(р,х) = И1,2(4Р) = ехр(-^Р(1 -х)) + ^[р(1 -х)ехр(д/р), (9)
(р, х) = ехр(-^(1 - х)) + ехр(-д/р(1 + х)) - ехр(-2^/р), (10)
^2,2 (р, х) = ехр(^7р (1 - х)) + ехр(-^/р (1 + х)) -1 + р(2х--ехр(-^тр), (11)
Wзl( р, х) = ^з,2( р, х) = ехр(-,/р (1 - х)) + ехр(^^р (1 + х)) - ехр{-Гр (3 - х)) -
= Гр (1 - х)ехр(-37р). (12)
Видно, что все члены аппроксимирующих решений (9) - (12) имеют аннотированный вид экспоненциальных функций от оператора дробного дифференцирования у[р и асимптотики аппроксимаций (9)-(12) при р ^ 0 и р ^ го совпадают с соответствующими представлениями (7), (8).
Решения (9) - (12) имеют форму, допускающую в большинстве важных для практики случаев переход к стандартным аналитическим функциям в пространстве оригиналов.
Так, для и (Г) = 1 решения задачи (2), (3) в пространстве оригиналов при W(р, х), отвечающих приближениям (9), (10), (11) (12), имеют вид
1 — х 1 — х 1
в 1 (х,т) = в2(х,т) = еф(—.=) + ~Г= ехр( - ), (13)
2л/т л\жт 4т
021(х, т) = ег/с(1~1х) + ег/с(1^) - ег/с(-^), (14)
2л/ т 2у1 т у/т
1 - х 1 + х 1 х2 -1 1
в22(х,т) = ег/с(—!=) + ег/с(-~1=) - ет/с(-тд-у= ехр(—), (15)
2>/т 2\ т у/т 2^жт3 т
п / \ п / \ гЛ - хч ГЛ + х 1 3 - х 1 - х 9 /1/-\
в31 (х,т) = 032 (х,т) = ет/с^—/=) + ег/с(—!=) - -ег/с(—1=) -^= ехр( - ) . (16)
2л1 т 2у1 т 2 2л/т у/пт 4т
На рисунке 1 представлены результаты сопоставления точного и приближенных, соответствующих
(13) - (16), решений для температурных полей и на рисунках 2, 3 - величины максимальных
относительных ошибок для температур и градиентов температур (тепловых потоков). В целом видна
удовлетворительная точность аппроксимативных решений и локально-экстремальный по времени
характер погрешностей составных разложений.
Максимальная ошибка ё Безразмерная температура 0
Координата х
1. Зависимость температуры от координаты для некоторых моментов времени
Безразмерное время т
Рис. 2. Зависимость максимальных по координате ошибок расчётов температур
от времени
При
0.5 1 2 3 4 567 10
Безразмерное время т Рис. 3. Зависимость максимальных по координате ошибок расчётов тепловых потоков
от времени
учете дальнейших членов приближенных представлений погрешность аппроксимаций
уменьшается.
Заключение
Применение предложенного подхода позволяет получать аналогичные приближенные решения для граничных условий 1, 2 и 3 рода для пластин, цилиндрических тел [3], [4], нелинейных процессов, задач Стефана [5] и иных краевых задач [6].
Литература
1.Бутковский А.Г. Характеристики систем с распределенными параметрами. М.: Наука, 1979. 224 с.
2.Маслов В.П. Теория возмущений и асимптотические методы. М.: МГУ, 1965. 549 с.
3. Дилигенский Н.В., Ефимов А.П. Структура аппроксимационных решений нестационарных задач теплообмена в телах конечных размеров // Труды III Минского международного форума. Минск: ИТМО им А.В.Лыкова, НАН Беларуси, 1996. Т. 9. С. 62-67.
4.Ефимов А.П. Метод построения равномерно-пригодных аппроксимаций решений нестационарных задач теплопроводности в телах конечных размеров // Вестник СамГТУ. 2008. Вып. 2(24). с.196-200.
5.Дилигенский Н.В., Ефимов А.П., Лившиц М.Ю. Применение метода возмущений для решения задачи Стефана в процессах промышленной теплофизики // Труды IV Минского международного форума. Минск: ИТМО им А.В. Лыкова, НАН Беларуси, 2000. Т. 3. С. 14-20.
6.Дилигенский Н.В., Ефимов А.П. Использование принципа длительности для конструирования систем математических моделей задач теплопроводности с требуемыми аппроксимативными свойствами // Труды 3-й Российской национальной конференции по теплообмену. М.: МЭИ, 2002. Т. 7. С. 111-114.
Дилигенский Николай Владимирович, доктор технических наук, профессор, заведующий кафедрой «Управление и системный анализ в теплоэнергетике» Самарского государственного технического университета, 443100, Самара, ул.Молодогвардейская 244, тел. (846) 332-42-34, e-mail: usat@ samgtu.ru
Ефимов Александр Порфирьевич, кандидат технических наук, доцент, докторант кафедры «Управление и системный анализ в теплоэнергетике» Самарского государственного технического университета.
Diligenskiy Nikolay Vladimirovich, Doctor of Technical Sciences, Professor, Samara State Technical University.
Efimov Alexandr Porfirjevich, Candidate of technical science, Associate Professor, Samara State Technical University.
