CC BY
101
ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1476 УДК 517.958:536.2
ОБ ОДНОМ МЕТОДЕ РЕШЕНИЯ НЕСТАЦИОНАРНЫХ ЗАДАЧ ТЕПЛОПРОВОДНОСТИ С НЕСИММЕТРИЧНЫМИ ГРАНИЧНЫМИ УСЛОВИЯМИ
И. В. Кудинов, Е. В. Стефанюк, М. П. Скворцова, Е. В. Котова, Г. М. Синяев
Самарский государственный технический университет, Россия, 443100, Самара, ул. Молодогвардейская, 244.
Аннотация
С использованием дополнительных граничных условий и дополнительной искомой функции в интегральном методе теплового баланса, получено приближенное аналитическое решение нестационарной задачи теплопроводности для бесконечной пластины при несимметричных граничных условиях первого рода. Решение имеет простой вид тригонометрического полинома с коэффициентами, экспоненциально стабилизирующимися во времени. С увеличением числа членов полинома получаемое решение приближается к точному. Введение зависящей от времени дополнительной искомой функции, задаваемой в одной из граничных точек, позволяет свести решение дифференциального уравнения в частных производных к интегрированию обыкновенного дифференциального уравнения. Дополнительные граничные условия находятся в таком виде, чтобы их выполнение искомым решением было эквивалентно выполнению исходного дифференциального уравнения в граничных точках. Показано, что выполнение уравнения лишь в граничных точках приводит к его выполнению и внутри области, минуя интегрирование по пространственной переменной, заменяемого выполнением искомым решением интеграла теплового баланса (осредненного дифференциального уравнения в частных производных). Отсутствие необходимости интегрирования исходного уравнения по пространственной переменной позволяет применять данный метод к решению нелинейных краевых задач с переменными начальными условиями и физическими свойствами среды и др.
© 2016 Самарский государственный технический университет. Образец для цитирования Кудинов И. В., Стефанюк Е. В., Скворцова М. П., Котова Е. В., Синяев Г. М. Об одном методе решения нестационарных задач теплопроводности с несимметричными граничными условиями // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. Т. 20, № 2. С. 342-353. аок 10.14498/vsgtu1476. Сведения об авторах
Игорь Васильевич Кудинов (к.т.н.; [email protected]), доцент, каф. теоретических основ теплотехники и гидромеханики.
Екатерина Васильевна Стефанюк (д.т.н., доц; [email protected]), профессор, каф. теоретических основ теплотехники и гидромеханики.
Марина Петровна Скворцова (аспирант; [email protected]), ассистент, каф. теоретических основ теплотехники и гидромеханики.
Евгения Валериевна Котова (к.т.н.; [email protected]; автор, ведущий переписку), доцент, каф. теоретических основ теплотехники и гидромеханики.
Геннадий Михайлович Синяев (к.т.н., доц; [email protected]), доцент, каф. теоретических основ теплотехники и гидромеханики.
Ключевые слова: нестационарная теплопроводность, бесконечная скорость распространения теплоты, неоднородные граничные условия, интеграл теплового баланса, дополнительная искомая функция, дополнительные граничные условия, тригонометрические координатные функции.
В теории теплопроводности известны методы, в которых используется понятие фронта температурного возмущения — глубины термического (прогретого) слоя [1—9]. То есть несмотря на то, что решению подлежит параболическое уравнение теплопроводности, в котором заложена бесконечная скорость распространения теплоты, вводится допущение о её конечности. При использовании данного метода процесс теплопроводности формально разделяется на две стадии по времени. Первая стадия заканчивается после достижения фронтом температурного возмущения противоположной стенки пластины. Во второй стадии изменение температуры происходит по всему объёму тела. Эта модель теплопроводности используется во многих методах: Т. Гуд-мена [3], М. Био [4], А. И. Вейника [5], М. Е. Швеца [6], Ю. С. Постоль-ника [7] и др. Их важным преимуществом является возможность получения простых по форме приближённых аналитических решений многих сложных краевых задач — нелинейных, с переменными коэффициентами, динамических и тепловых задач пограничного слоя и др. Однако к их существенным недостаткам относится малая точность получаемых решений. С целью повышения точности требуется увеличивать число членов аппроксимационного ряда, для определения неизвестных коэффициентов которого используются различные методы. В частности, в работах [1, 10-14] используются дополнительные граничные условия, задаваемые в граничных точках и на фронте температурного возмущения, позволяющие существенно повысить точность получаемых решений.
В работах [1, 12-14] показано, что с увеличением числа приближений n наряду с возрастанием точности решений происходит уменьшение времени продвижения фронта температурного возмущения Foi от поверхности до центра пластины (симметричная задача). И в пределе при n ^ œ>, Foi ^ 0, что подтверждает бесконечную скорость распространения теплоты. При этом было отмечено, что с увеличением числа приближений в связи с уменьшением времени первой стадии процесса её значение в оценке температурного состояния конструкции уменьшается, а второй — возрастает. В связи с этим в настоящей работе рассматривается метод решения, в котором исключается первая стадия процесса. Основную его идею рассмотрим на примере решения следующей краевой задачи (рис. 1):
д 6(£, Fo) д 26({, Fo) . n п Л ,,
д{2 ) (Fo>0<^<1); (1) 6(е, °) = (2)
6(0, Fo) = 1; (3)
6(1, Fo) = 0; (4)
где 6 = (T — T0)/(Tc1 — То) — безразмерная температура; T0 — начальная температура; Tc1 — температура стенки при x = 0; Tc2 = Т0 — температура стенки
[Figure 1. A calculation scheme for a heat transfer problem]
при х = 5; х — координата; 5 — толщина пластины; £ = х/5 — безразмерная координата; Ро = а£/52 — число Фурье (безразмерное время); £ — время; а — коэффициент температуропроводности.
В точке £ = 1 введем дополнительную искомую функцию вида
. , дв(£, Ро) 9(ро)= (££ ) ?=1 =tg(5)
где а — зависящий от времени угол наклона температурной кривой к оси £. Минимальное значение этот угол (а = 0) будет иметь при Ро = 0, а максимальное (а = п/2) — при Ро ^ те.
Ввиду бесконечной скорости распространения теплоты, описываемой параболическим уравнением теплопроводности (1), угол а будет увеличиваться сразу после приложения граничного условия первого рода в точке £ = 0 и далее будет возрастать во времени, устремляясь к конечному значению а = п/2 при Ро ^ те, то есть при установлении стационарного режима теплообмена. Следовательно, процесс изменения производной от искомой функции по координате £ в точке £ = 1 будет включать весь диапазон времени нестационарного процесса 0 < Ро < те.
Решение задачи (1)—(4) принимается в виде
в(£, Fo) = 1 + £ bk(q)<pk(£), (6)
k=i
где bk(q), k = 1, n, — неизвестные коэффициенты;
,-k (£)=«.(ri (i - 2))
— координатные функции.
Очевидно, что соотношение (6) благодаря принятой системе координатных функций удовлетворяет граничному условию (3). Для определения неизвестных коэффициентов будем использовать основное граничное условие (4), соотношение (5) и некоторые дополнительные граничные условия, задаваемые в точках £ = 0 и £ = 1 и определяемые таким образом, чтобы их
выполнение искомым решением (6) было эквивалентно выполнению уравнения (1) в этих точках [1, 12-14].
Найдём дополнительные граничные условия применительно к точке £ = 0. Для этого продифференцируем граничное условие (3) по переменной Ро:
дв(0, Ео) дРо
0.
Сравнивания полученное соотношение с уравнением (1) применительно к точке £ = 0, получаем дополнительное граничное условие вида
д2в(£, Ро)
д£2
?=о
0.
(7)
Продифференцируем соотношение (7) по переменной Ро и запишем полученное в виде
д2 / дв(£, Ро) \ д£П дРо )
0.
?=о
(8
Соотношение (8) с учётом уравнения (1) приводится к следующему дополнительному граничному условию:
д4в(£, Ро)
д£4
?=о
0.
(9)
Аналогично дифференцируя соотношение (9) по переменной Ро, с учётом уравнения (1) получаем ещё одно дополнительное граничное условие:
д6в(£, Ро)
д£6
?=о
0.
(10)
На основе соотношений (7), (9), (10) можно записать общую формулу дополнительных граничных условий в точке £ = 0:
д*в(£, Ро)
д£г
?=о
0, г = 2, 4, 6, ....
(11)
Отметим, что благодаря принятой системе координатных функций все дополнительные граничные условия, определяемые формулой (11), решением (6) выполняются при любом числе приближений. Для каких-либо других, более сложных краевых задач (нелинейных, с переменными коэффициентами и др.), дополнительные граничные условия в точке £ = 0 принятым решением могут и не выполняться заранее. Поэтому их выполнение следует проводить через определение неизвестных коэффициентов Ъ^(д) искомого решения вида (6).
Для нахождения дополнительных граничных условий в точке £ = 1 продифференцируем соотношения (4), (5) по переменной Ро:
дв(1, Ро)
дРо
(12) 345
йд(Ро) д2в(£, Ро)
й Ро
д£дРо
?=1
(13)
Сравнивая соотношение (12) с уравнением (1), находим первое дополнительное граничное условие вида
д2в(£, Ро)
д£2
Перепишем соотношение (13) в виде
?=1
0.
йд(Ро) _ д / дв(£, Ро) йРо = V дРо
?=1
(14)
(15)
Соотношение (15) с учётом уравнения (1) приводится к следующему дополнительному граничному условию:
д 3в(£, Ро)
д£3
?=1
йд(Ро) й Ро '
(16)
Продифференцируем соотношения (14), (16) по переменной Ро и запишем полученное в виде
д2 / дв(£, Ро) д£П дРо д3 / дв(£, Ро) V дРо
?=1
?=1
0; 0.
(17)
(18)
Соотношения (17), (18) с учётом уравнения (1) приводятся к следующим дополнительным граничным условиям:
д4в(£, Ро)
д£4
д 5в(£, Ро)
?=1
д£5
?=1
0;
й2 д(Ро) й Ро2 '
(19)
(20)
Аналогично, дифференцируя соотношения (19) и (20) по переменной Ро, с учётом уравнения (1) получаем дополнительные граничные условия вида
д6в(£, Ро)
д£6
д 7в(£, Ро)
?=1
д£7
?=1
0;
й3 д(Ро) йРо3 '
(21) (22)
Анализируя соотношения (14), (19), (21), а также (16), (20), (22), можно записать общие формулы дополнительных граничных условий применительно к точке £ = 1:
дгв(£, Ро)
де
?=1
0, г = 2, 4, 6, . . . ,
(23)
в2-+1в(£, Р°) _(Лг(Го) , = 1, 2, 3 . .. (24)
д£2г+1
£=1 агРо
Очевидно, что граничное условие (4) и условие (5), а также все дополнительные граничные условия, определяемые по общим формулам (23), (24), решением (6) не выполняются. Все эти условия будут использованы для нахождения неизвестных коэффициентов (9) решения (6).
Для получения решения в первом приближении подставим (6) (ограничиваясь двумя членами суммы) в соотношения (4) и (5):
1 + Ы^1(1) + Ь2(^2 (1) = 0; 1 91 (Ро) - (1)/д£ - Ь2(д)д^2(1)/д£ = 0. / ( )
Соотношения (25) относительно неизвестных коэффициентов 61(9) и 62(9) представляют собой систему двух алгебраических линейных уравнений, из решения которой находим
6,(9)= 3П), 62(9) = Л2^.
4п 4п
Подставляя найденные значения 61(9) и 62(9) в (6), получаем
/2
в(£, Ро) = 1 + 42 [(4^ - 3п) (1 - |) ) + (49 - п) еов (ЗП (1 - £ ) ) ] • (26
Потребуем, чтобы соотношение (26) удовлетворяло осредненному в пределах толщины пластины уравнению (1), то есть интегралу теплового баланса вида
Г ^ = Г ^-е- (27)
Подставляя (26) в (27), относительно неизвестной функции 9(Ро) получаем следующее обыкновенное дифференциальное уравнение:
А А29(Ро) + Аз = 0, (28)
аРо
где А1 = 128 - 64/2, А2 = 12п2(2 + /2), А3 = 9п3/2. Интегрируя уравнение (28), находим
9(Ро) = С1 ехр(^Ро) - (29)
где С1 — постоянная интегрирования,
_ 3п2(2 + /2) 3п/2
" 16(л/2 - 2) ' ^ 4(л/2 + 2)' Подставляя (29) в (26), получаем
в(£, Ро) = 1 + (/^в"Ро - + М) еов(2 (1 - 2)) +
Для определения постоянной интегрирования С1 составим невязку начального условия (2) и потребуем ортогональности невязки к координатной функции ^1(£):
[ в(£, 0)^(£К = 0. (31)
о
Подставляя (30) в (31), получаем
Соотношение (32) ввиду ортогональности косинусов приводится к виду
да« - "Н3^2^ - |))]й£ = ^
=- Ю Ч2(1 - 2Ж (33)
После вычисления интегралов в (33) относительно неизвестного коэффициента С1 получаем алгебраическое уравнение, из решения которого находим С1 = п(ЗП(2 + ^2) - 16)/16.
После определения постоянной интегрирования решение задачи (1)-(4) в первом приближении находится из соотношения (30). Это решение точно удовлетворяет граничным условиям (3), (4) и интегралу теплового баланса (27). Уравнение (1) и начальное условие (2) в данном случае выполняются лишь приближенно в (первом приближении). Отметим, что число приближений будем определять не числом членов суммы решения (6), а порядком дифференциального уравнения относительно неизвестной функции д(Ро) или, что то же самое, числом констант интегрирования этого уравнения.
Анализ результатов расчётов по формуле (30) позволяет заключить, что при 0.1 ^ Ро < те их расхождение с точным решением [15] не превышает 4%. Повышение точности решения связано с увеличением числа членов ряда (6), неизвестные коэффициенты которого будем определять из условий (4), (5), (23), (24).
Решение задачи (1)-(4) во втором приближении принимается в виде
в(£ Ро)=1+Ч П (1 - 2))+ъ Ч | (1 - 2))+
+ Ь3-(| (1 - 2)) + Ъ4со8(7П(1 - 2)). (34)
Подставляя (34) в условия (4), (5), (23) (при г = 2), (24) (при г = 1), относительно неизвестных коэффициентов Ъ^ будем иметь систему четырёх
алгебраических уравнений, из решения которой получаем
Ь1 = £к/ + 29В2д - 105Вз, 62 = Вд' + 37В2д + 35Вз,
1 11 1 13 (35)
Ьз = ^В1д' - 11- 7Вз, 64 = 3+ 13В2д + 5Вз.
Здесь д' = ^д/^Ро; В = ^2/(2пз); В2 = ^2/(32п); Вз = ^2/128.
Подставляя (34) (с учётом (35)) в интеграл теплового баланса (27), относительно неизвестной функции д(Ро) будем иметь обыкновенное дифференциальное уравнение второго порядка
Ад'' + Ад' - Взд - В4 = 0, (36)
в котором
В1 = 92/6^2 - 16384, В2 = п2(10624^2 - 34816), Вз = п4(6860^2 + 6720), В4 = 3675п5 ^2, д'' = ^2д/^Ро2.
Интегрируя уравнение (36), находим
д(Ро) = С1 ехр(^1Ро) + С2 ехр(,г2Ро) + В, (37)
где
, п2(2^22938 - 14368^2 ^ 83^2 ± 272) 105п^2 ^1,2 = ±-_ „,„ ^-—-, В — •
16(9л/2 - 16) ' 196\/2 +192'
Ci, C2 — константы интегрирования, определяемые из начального условия (2). Для этого составляется невязка начального условия и требуется ортогональность невязки к координатным функциям <^i(£) и ^>2(£):
Г в(е, о)^^(Ode = 0, j = 1,2. (38)
Jo
Подставляя (34) (с учётом (35), (37)) в (38), после определения интегралов относительно Ci и C2 получаем систему двух алгебраических линейных уравнений, из решения которой находим C1 = -1.490559; C2 = 2.022993.
Результаты расчётов по формуле (34) в сравнении с точным решением [15] приведены на рис. 2. Их анализ позволяет заключить, что во втором приближении решение существенно уточняется. И, в частности, при 0.09 ^ Fo < те расхождение с точным решением уменьшается с 6 % (в первом приближении) до 2 % — во втором.
Результаты расчётов в третьем и четвёртом приближениях в сравнении с точным решением приведены на рис. 2. Их анализ позволяет заключить, что при 0.04 ^ Fo < те расхождение четвёртого приближения с точным решением не превышает 1.5 %. Невязки уравнения (1) и начального условия (2) с увеличением числа приближений уменьшаются, что свидетельствует о сходимости предлагаемого метода решения.
1.0 0
0.8
0.6
0.4
0.2
0
О
0.2
0.4
0.6
0.8
£ 1.0
Рис. 2. Графики изменения температуры: сплошная линия — точное решение; □ —второе приближение; д —третье приближение; ◦ —четвёртое приближение; метки: 1— Fo = 1.0;
2 — Fo = 0.2; 3 — Fo = 0.09; 4 — Fo = 0.04 [Figure 2. Graphs of the temperature change: the exact solution (solid curve), the second approximation of the solution (□), the third approximation of the solution (д), and the fourth approximation of the solution (◦); Labels: 1 — Fo = 1.0; 2 — Fo = 0.2;
3 — Fo = 0.09; 4 — Fo = 0.04]
Таким образом, основным преимуществом рассмотренного метода является возможность сведения решения исходного дифференциального уравнения в частных производных к интегрированию обыкновенного дифференциального уравнения относительно дополнительной искомой функции, зависящей лишь от времени, что позволяет применять данный метод к решению сложных краевых задач (нелинейных, с переменными коэффициентами и др.), дифференциальные операторы которых не допускают разделения переменных, а также для задач с переменными во времени граничными условиями, с переменными начальными условиями, с источниками теплоты и др.
Выводы. Получено приближённое аналитическое решение нестационарной задачи теплопроводности для бесконечной пластины с несимметричными граничными условиями первого рода. Введение дополнительной искомой функции q(Fo), представляющей изменение во времени угла наклона температурной кривой к оси пространственной переменной, позволяет свести решение уравнения в частных производных к интегрированию обыкновенного дифференциального уравнения. Использование функции q(Fo) обосновывается описываемой параболическим уравнением теплопроводности бесконечной скорости распространения теплоты, согласно которой угол наклона температурной кривой к оси £ в точке £ = 1 начинает возрастать сразу после приложения граничного условия первого рода в точке £ = 0, и, следовательно, диапазон изменения функции q(Fo) включает весь диапазон времени нестационарного процесса 0 < Fo < те.
Дополнительные граничные условия находятся в таком виде, чтобы их выполнение искомым решением было эквивалентно выполнению исходного дифференциального уравнения (в частных производных) в граничных точках
области. Показано, что выполнение уравнения на границах приводит к его выполнению и внутри рассматриваемой области с точностью, зависящей от числа приближений принятого решения.
Благодарности. Работа выполнена при поддержке Российского фонда фундаментальных исследований (проект № 16-38-00059-мол_а).
ORCIDs
Игорь Васильевич Кудинов: http://orcid.org/0000-0002-9422-0367 Екатерина Васильевна Стефанюк: http://orcid.org/0000-0003-1494-8161 Марина Петровна Скворцова: http://orcid.org/0000-0003-1812-4178 Евгения Валериевна Котова: http://orcid.org/0000-0003-4423-2738 Геннадий Михайлович Синяев: http://orcid.org/0000-0001-9094-2306
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кудинов В. А., Стефанюк Е. В. Задачи теплопроводности на основе определения фронта температурного возмущения// Известия РАН. Энергетика, 2008. №5. С. 141-157.
2. Лыков А. В. Методы решения нелинейных уравнений нестационарной теплопроводности // Известия АН СССР. Энергетика и транспорт, 1970. №5. С. 109-150.
3. Goodman T. R Application of integral methods to transient nonlinear heat transfer // Advances in Heat Transfer, 1964. Т. 1. С. 51-122. doi: 10.1016/S0065-2717(08)70097-2.
4. Biot M. A. Variational Principles in Heat Transfer: A Unified Lagrangian Analysis of Dis-sipative Phenomena/ Oxford Mathematical Monographs. Oxford: Clarendon Press, 1970. x+185 pp.
5. Вейник А. И. Приближенный расчет процессов теплопроводности. М., Л.: Госэнерго-издат, 1959. 184 с.
6. Швец М. Е. О приближенном решении некоторых задач гидродинамики пограничного слоя // ПММ, 1949. Т. 13, №3. С. 257-266.
7. Тимошпольский В. И., Постольник Ю. С., Андрианов Д. Н. Теоретические основы теплофизики и термомеханики в металлургии. Минск: Белорусская наука, 2005. 560 с.
8. Глазунов Ю. Т. Вариационные методы. М., Ижевск: НИЦ "Регулярная и хаотическая динамика"; Институт компьютерных исследований, 2006. 470 с.
9. Беляев Н. М., Рядно А. А. Методы нестационарной теплопроводности. М.: Высш. шк., 1978. 328 с.
10. Федоров Ф. М. Граничный метод решения прикладных задач математической физики. Новосибирск: Наука, 2000. 220 с.
11. Кудряшев Л. И., Меньших Н. Л. Приближенные решения нелинейных задач теплопроводности. М.: Машиностроение, 1979. 232 с.
12. Кудинов В. А., Стефанюк Е. В. Аналитический метод решения задач теплопроводности на основе введения фронта температурного возмущения и дополнительных граничных условий // Инженерно-физический журнал, 2009. Т. 82, №3. С. 540-558.
13. Стефанюк Е. В., Кудинов В. А. Получение приближенных аналитических решений при рассогласовании начальных и граничных условий в задачах теории теплопроводности// Изв. вузов. Матем., 2010. №4. С. 63-71.
14. Кудинов В. А., Кудинов И. В., Скворцова М. П. Обобщенные функции и дополнительные граничные условия в задачах теплопроводности для многослойных тел // Ж. вычисл. матем. и матем. физ., 2015. Т. 55, №4. С. 669-680. doi: 10.7868/ S0044466915040080.
15. Лыков А. В. Теория теплопроводности. М.: Высш. шк., 1967. 600 с.
Поступила в редакцию 12/II/2016; в окончательном варианте — 29/III/2016; принята в печать — 27/V/2016.
KyguHOB M. B., CTetfcaHMK E. B., Ckbo pn,OB a M. n., KoTOBa E. B., CuHaeB r. M.
Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki
[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 342-353 ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1476
MSC: 80A20
ON ONE METHOD FOR SOLVING TRANSIENT HEAT CONDUCTION PROBLEMS WITH ASYMMETRIC BOUNDARY CONDITIONS
I. V. Kudinov, E. V. Stefanyuk, M. P. Skvortsova,
E. V. Kotova, G. M. Sinyaev
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.
Abstract
Using additional boundary conditions and additional required function in integral method of heat-transfer we obtain approximate analytical solution of transient heat conduction problem for an infinite plate with asymmetric boundary conditions of the first kind. This solution has a simple form of trigonometric polynomial with coefficients exponentially stabilizing in time. With the increase in the count of terms of a polynomial the obtained solution is approaching the exact solution. The introduction of a time-dependent additional required function, setting in the one (point) of the boundary points, allows to reduce solving of differential equation in partial derivatives to integration of ordinary differential equation. The additional boundary conditions are found in the form that the required solution would implement the additional boundary conditions and that implementation would be equivalent to executing the original differential equation in boundary points. In this article it is noted that the execution of the original equation at the boundaries of the area only (via the implementation of the additional boundary conditions) leads to the execution of the original equation also inside that area. The absence of direct integration of the original equation on the spatial variable allows to apply this method to solving the nonlinear boundary value problems with variable initial conditions and variable physical properties of the environment, etc.
Keywords: transient heat conductivity, infinite speed of heat propagation, integral method of thermal balance, inhomogeneous boundary conditions, approximate analytical solution, additional required function, additional boundary conditions, trigonometric coordinate functions.
© 2016 Samara State Technical University. Please cite this article in press as:
Kudinov I. V., Stefanyuk E. V., Skvortsova M. P., Kotova E. V., Sinyaev G. M. On one method for solving transient heat conduction problems with asymmetric boundary conditions, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 342-353. doi: 10.14498/vsgtu1476. (In Russian) Authors Details:
Igor V. Kudinov (Cand. Tech. Sci.; [email protected]), Associate Professor, Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics.
Ekaterina V. Stefanyuk (Dr. Techn. Sci.; [email protected]), Professor, Dept. ofTheoret-ical Fundamentals of Heat-Engineering and Hydromechanics.
Marina P. Skvortsova (Postgraduate Student; [email protected]), Assistant, Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics. Eugeniya V. Kotova (Cand. Techn. Sci.; [email protected]; Corresponding Author), Associate Professor, Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics.
Gennady M. Sinyaev (Cand. Techn. Sci.; [email protected]), Associate Professor, Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics.
Acknowledgments. This work was supported by the Russian Foundation for Basic Research
(project no. 16-38-00059-mol_a).
ORCIDs
Igor V. Kudinov: http://orcid.org/0000-0002-9422-0367 Ekaterina V. Stefanyuk: http://orcid.org/0000-0003-1494-8161 Marina P. Skvortsova: http://orcid.org/0000-0003-1812-4178 Eugeniya V. Kotova: http://orcid.org/0000-0003-4423-2738 Gennady M. Sinyaev: http://orcid.org/0000-0001-9094-2306
REFERENCES
1. Kudinov V. A., Stefanyuk E. V. A heat conduction problem by determining the temperature perturbation front, Izvestiia RAN. Energetika, 2008, no. 5, pp. 141-157 (In Russian).
2. Lykov A. V. Methods for solving nonlinear equations of nonstationary heat conduction, Izvestiia AN SSSR. Energetika i transport, 1970, no. 5, pp. 109-150 (In Russian).
3. Goodman T. R Application of integral methods to transient nonlinear heat transfer, Advances in Heat Transfer, 1964, vol.1, pp. 51-122. doi: 10.1016/S0065-2717(08)70097-2.
4. Biot M. A. Variational Principles in Heat Transfer: A Unified Lagrangian Analysis of Dis-sipative Phenomena, Oxford Mathematical Monographs. Oxford, Clarendon Press, 1970, x+185 pp.
5. Veinik A. I. Priblizhennyi raschet protsessov teploprovodnosti [An approximate calculation of heat conduction processes]. Moscow, Leningrad, Gosenergoizdat, 1959, 184 pp. (In Russian)
6. Shvets M. E. The approximate solution of some problems of the hydrodynamics of the boundary layer, Prikladnaia matematika i mekhanika, 1949, vol. 13, no. 3, pp. 257-266 (In Russian).
7. Timoshpol'skii V. I., Postol'nik Yu. S., Andrianov D. N. Teoreticheskie osnovy teplofiziki i termomekhaniki v metallurgii [Theoretical foundations of thermal physics and metallurgy thermomechanics]. Minsk, Belorusskaia nauka, 2005, 560 pp. (In Russian)
8. Glazunov Yu. T. Variatsionnye metody [Variation Methods]. Moscow, Izhevsk, Reguliarnaia i khaoticheskaia dinamika; Institut komp'iuternykh issledovanii, 2006, 470 pp. (In Russian)
9. Belyaev N. M., Ryadno A. A. Metody nestatsionarnoi teploprovodnosti [Methods of Transient Heat Conduction]. Moscow, Vyssh. shk., 1978, 328 pp. (In Russian)
10. Fedorov F. M. Granichnyi metod resheniia prikladnykh zadach matematicheskoi fiziki [The boundary method for solving applied problems of mathematical physics]. Novosibirsk, Nauka, 2000, 220 pp. (In Russian)
11. Kudryashev L. I., Men'shikh N. L. Priblizhennye resheniia nelineinykh zadach teploprovod-nosti [Approximate solution of the nonlinear heat conduction problems]. Moscow, Mashinos-troenie, 1979, 232 pp. (In Russian)
12. Kudinov V. A., Stefanyuk E. V. Analytical solution method for heat conduction problems based on the introduction of the temperature perturbation front and additional boundary conditions, Journal of Engineering Physics and Thermophysics, 2009, vol. 82, no. 3, pp. 537555. doi: 10.1007/s10891-009-0223-8.
13. Stefanyuk E. V., Kudinov V. A. Approximate analytic solution of heat conduction problems with a mismatch between initial and boundary conditions, Russian Math. (Iz. VUZ), 2010, vol.54, no. 4, pp. 55-61. doi: 10.3103/S1066369X10040079.
14. Kudinov V. A., Kudinov I. V., Skvortsova M. P. Generalized functions and additional boundary conditions in heat conduction problems for multilayered bodies, Comput. Math. Math. Phys., 2015, vol.55, no. 4, pp. 666-676. doi: 10.1134/S0965542515040089.
15. Lykov A. V. Teoriia teploprovodnosti [Theory of heat conduction]. Moscow, Vyssh. shk., 1967, 600 pp. (In Russian)
Received 12/II/2016;
received in revised form 29/III/2016;
accepted 27/V/2016.
