Распространение тепла в бесконечном твердом теле от цилиндрической полости при конвективном теплообмене

М. А. Осипенко и др. Распространение тепла в бесконечном твердом теле от полости

УДК 622.4:536.21

РАСПРОСТРАНЕНИЕ ТЕПЛА В БЕСКОНЕЧНОМ ТВЕРДОМ ТЕЛЕ ОТ ЦИЛИНДРИЧЕСКОЙ ПОЛОСТИ ПРИ КОНВЕКТИВНОМ ТЕПЛООБМЕНЕ

М. А. Осипенко, О. И. Дударь, Е. С. Дударь

Пермский национальный исследовательский политехнический университет,

кафедра теоретической механики E-mail: olegdudar@yandex.ru

Методом функций Грина в сочетании с преобразованием Лапласа получено решение задачи о распространении тепла в бесконечном твердом теле от движущегося в цилиндрической полости газа.

Ключевые слова: цилиндрическая полость, нестационарная теплопроводность, конвективный теплообмен, функция Грина, преобразование Лапласа, функции Бесселя и Неймана.

The Heat Conductivity in the Infinite Solid of the Convection in a Cylindrical Cavity

M. A. Osipenko, O. I. Dudar, E. S. Dudar

Perm National Research Politechnical University,

Department of Theoretical Mechanics E-mail: olegdudar@yandex.ru

The problem of heat conductivity in the infinite solid from gas moving in a cylindrical cavity is solved by the combination of the Green function method and the Laplace transformation,

Key words: cylindrical cavity, nonstationary heat conductivity, convective heat transfer, Green’s function, Laplace transformation, Bessel’s and Neumann’s functions.

ВВЕДЕНИЕ

Задача распространения тепла в бесконечном пространстве от движущегося в цилиндрической полости газа является одной из основных в горной теплофизике и теплофизике подземных сооружений. Ее решение позволяет рассчитывать тепловой режим горных выработок, тоннелей, метрополитенов и других подземных сооружений. Н. Б. СагБ1а-№ [1], А. Н. Щербань и О. А. Кремнев [2] получили решение данной задачи с помощью преобразования Лапласа для случая неизменной во времени температуры движущегося в полости газа. В работах А. С. Галицына и А. Н. Жуковского [3], а также А. С. Галицына [4] для получения решения использовалось модифицированное преобразование Вебера при произвольном законе изменения температуры движущегося в полости газа и произвольном начальном условии. А. Ф. Галкин и Ю. А. Хохолов [5] задачу решали численно в сопряженной постановке, то есть одновременно для твердой и газовой фаз.

В настоящей статье задача решена для случая произвольного изменения температуры движущегося в полости газа и произвольного начального условия, то есть в той же постановке, что и в работах [3,4]. Однако для решения применялся метод функций Грина в сочетании с преобразованием Лапласа. Найденное решение позволило выявить допущенную в работе [4] неточность.

1. ПОСТАНОВКА ЗАДАЧИ

Рассмотрим задачу передачи тепла от газа (воздуха), движущегося в цилиндрической полости, окружающему её бесконечному однородному изотропному твердому телу (горному массиву). Полагаем, что влияние осевой координаты отсутствует, поэтому изменение температуры в горном массиве описывается нестационарным уравнением теплопроводности, записанным в полярной системе координат:

дт 1 д / дт\

Ь > 0,

дТ _ 1 д ( дТ\ dt ar dr V dr /’

r > r0.

коэффициент температуропроводности горного массива, ro

радиус

где Т — температура, а полости.

Будем полагать, что в начальный момент времени температура в массиве с полостью является произвольной функцией радиуса

Т |г=о = То

но при этом на бесконечности как в начальный, так и в любой другой момент времени она равна Т^:

0 I r ——^с

_ Т„

© Осипенко М. А., Дударь О. И., Дударь Е. С., 2012

89

Изв. Сарат. ун-та. Нов. сер. 2012. Т. 12. Сер. Математика. Механика. Информатика, вып. 1

где Тс — естественная температура горных пород, то есть температура пород при отсутствующей полости.

В качестве граничных условий примем условие конвективного теплообмена на стенке полости

дТ \

-дг + Н(Т - Та

= 0,

где Н = а/Л, а — коэффициент теплоотдачи, Л — коэффициент теплопроводности горного массива, Та(*) — закон изменения температуры воздуха в полости.

По аналогии с работами [1, 6, 7] решение задачи (1)-(4) можно записать с использованием функции Грина д(г', г,*):

Т(г, *) = Тс + аНго д(го,г,* - т)[Та(т) - Тж]^т + д(г',г,*)[То(г') - Тж]г'¿г'.

Тогда решение исходной задачи сводится к определению функции Грина, которая может быть найдена из решения следующей сопряженной задачи:

= ^я7' г'я¡7 + Г6(г' - гМ), г' > го, * > 0

д\г=о = 0,

- 0+кя д\г' — с 0;

= 0,

здесь через 6 обозначена дельта-функция Дирака.

Решение задачи в такой постановке позволяет не конкретизировать вид функций То(г) и Та(*), входящих в начальное и граничное условия.

3. РЕШЕНИЕ ЗАДАЧИ ДЛЯ ИЗОБРАЖЕНИЯ ФУНКЦИИ ГРИНА

Применив к (6) преобразование Лапласа [8] по переменной *, получаем следующую задачу для изображения функции Грина С(г',г,р):

1 д / 'дС\ 1 '

рС = а—— г’— + -6(г' - г), г' > го,

' ' '

г' дг' V дг' / г

= 0,

д£

- от + НС

С\г' —— с 0;

где р — комплексный параметр.

Вводя новые переменные

х = вг, х' = вг', хо = вго, в = а/р/а.

и учитывая в (7) дельта-функцию «явно», приходим к следующей системе уравнений:

,д2 С ,дС

+ х'—— - х С = 0, \хо\ < \х'\ < \х\, \х'\ > \х\.

дх'2 дх'

С\х' = х+о - С\х' =х—о = °;

дС

дх'

1С\:

х'=х+о

дС 1, „ - а;7 + ;НС = 0.

дх'

х' =х—о = 0,

1

ах

7" = Т'о

ос

Г =Го

М. А. Осипенко и др. Распространение тепла в бесконечном твердом теле от полости

Первое уравнение в системе (9) — это модифицированное уравнение Бесселя с индексом нуль. Его решение имеет вид [9,10]:

Є =

С/о(х/) + БКо(х;), ж'| > |ж|.

где /о, Ко — функции нулевого порядка: модифицированная Бесселя и Макдональда соответственно.

Так как /о(х') ^ го при х' ^ го [7], то из последнего уравнения системы (9) сразу следует, что С = 0. Далее, используя известные свойства рассматриваемых цилиндрических функций [11,12]

/о (х) = /х(х), ко (х) = -Кх(х),

где /1, К — функции первого порядка: модифицированная Бесселя и Макдональда соответственно, приходим к следующей системе линейных алгебраических уравнений относительно неизвестных А, В, Я:

- А/о(х) - ВКо (х) + ЯКо(х) = 0,

-А/1 (х) + ВК1 (х) - ЯК1 (х) = -—,

ах

, А [-/1 (хо) + (Н/в)/о(хо)] + В [К1(хо) + (Н/в)Ко(хо)] = 0. Учитывая значение определителя Вронского [10]

/о(х)К1(х) + /1 (х)Ко (х) = -,

решение системы (10) можем записать в виде

А =

Ко

В = А

/1 (жо) - (^/а)/о(жо) Кі (жо) + (ЛУз)Ко(жо).

Б =

/о

+ В.

Возвращаясь к старым переменным, получим следующее решение для изображения функции Грина:

Є(г/, Г, 5) = <

г Ко

Ко

, Го < Г/ < Г,

, Г < Г/,

где К(вго) = вК1 (вго) + НКо(вго), /(вго) = -в/1(вго) + Н/о(вго). Если принять обозначения: т = ш1п(г',г), М = шах(г',г), то решение можно записать одной формулой

Є(г/, г, 5) =

Ко(5М)[/о(5т)К(5Го) - Ко(5Ш)/(5Го)]

аК

4. РЕШЕНИЕ ЗАДАЧИ ДЛЯ ФУНКЦИИ ГРИНА

Для определения оригинала по изображению (11) воспользуемся формулой обращения преобразования Лапласа [8]:

д(г/, г, £) = 2~т І Є(г/,г,р)вр^ ф.

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

В формуле (12) интеграл берется по мнимой оси. Чтобы привести его к более удобному виду, не содержащему комплексных чисел, перейдем к новой переменной интегрирования £ = ^¿в = \/р/а,

где нижний знак относится к нижней мнимой полуоси, а верхний — к верхней. При этом для переменной £ нижняя полуось перейдет в полуось Г1, а верхняя — в полуось Г2 (рисунок).

а

а

— 1Ж

Механика

91

Изв. Сарат. ун-та. Нов. сер. 2012. Т. 12. Сер. Математика. Механика. Информатика, вып. 1

У ‘ г;

0 П

Деформирование контура интегрирования Гі и Г2 в контур Г1 и Г2

По аналогии с работой [8] можно показать, что знаменатель функции Є(г/,г,р) не имеет нулей, поэтому выбранная для Є(г/,г,р) ветвь аналитична на всей комплексной плоскости за исключением начала координат с разрезом по отрицательной вещественной полуоси. Для аналитичной функции контуры Г1 и Г2 можно деформировать в контуры Г1, и Г2 (см. рисунок), делая тем самым £ вещественным положительным числом. Тогда формулу (12) можно привести к виду

д(г'.г.і) = -( І Є(г',г,8)|,=_і{е а{2*£й£ + І С(г',г,8)|,=і{е £й£

Подставив (11) в (13) и используя формулы связи [10-12]:

/о(±гх) = 7о (х), /1 (±гх) = ±г71(х),

¿п п

Ко (±гх) = у[т^о (х) + ¿N3 (х)], К1 (±гх) = ^[-^1 (х) ± ¿N1 (х)],

где х — вещественное положительное число, 7о, 71 — функции Бесселя, а N, N1 — функции Неймана нулевого и первого порядка соответственно, получим следующее выражение для оригинала:

сс

г

’ ’ У .72(£го)+ №'2(£го) ’

о

где 7(£го) = £Л(£го) + Н7о(£го), N(£го) = £^(£го) + НNо(£го).

Тем самым определена функция Грина для задачи (6). В частности, если принять г' = го, то т = шт(го,г) = го, М = шах(го,г) = г. Тогда, учитывая значение определителя Вронского [10]

2

^з(х)71 (х) - 7о(х)^(х) = —,

пх

получим следующее выражение для функции Грина при условии г' = го:

д(го-г'*) = *(еЖ+- дат е_*2 *

о

Подставляя формулы (14) и (16) в (5) получим решение исходной задачи (1)-(4) распространения тепла в бесконечном горном массиве от движущегося в цилиндрической полости воздуха:

с ( I Л

Т (г,*) = Тс + 2ПН / ^ £((£У£П’) ^ е—а^2 <‘—т > [Та(т) - Тс] ¿Т [ £ +

. . ДО)(£ГУ(£го) - (£го)]ДО(СГУ(£го) - ^о(СГЖ(£го)] е_а£2*£,£[т (г/) т ]г/ ,/

+ // 72(Сго) + N2 (£го) е £^[і°(Г) іо]ГгіГ'

Го о

о

ОС

А. С. Устинова. Влияние проскальзывания на вискозиметрическое течение

ЗАКЛЮЧЕНИЕ

Рассмотрена задача о распространении тепла в бесконечном пространстве (горном массиве) от цилиндрической полости, по которой движется поток газа (воздуха). Методом функций Грина в сочетании с преобразованием Лапласа получено решение этой задачи при произвольном законе изменения температуры воздуха и произвольном начальном условии. Решение может быть использовано для расчета теплового режима подземных сооружений с формой, близкой к цилиндрической.

Библиографический список

1. Carslaw H. S., Jaeger J. C. Conduction of Heat in Solids. Oxford : Clarendon Press, 1959. 450 p.

2. Щербань А. Н, Кремнев О. А. Научные основы расчета и регулирования теплового режима глубоких шахт : в 2 т. Киев : Изд-во АН УССР, 1959. Т. 1. 425 с.

3. Галицын А. С., Жуковский А. Н. Интегральные преобразования и специальные функции в задачах теплопроводности. Киев : Наук. думка, 1976. 286 с.

4. Галицын А. С. Краевые задачи теплофизики подземных сооружений. Киев : Наук. думка, 1983. 236 с.

5. Галкин А. Ф., Хохолов Ю. А. Теплоаккумулирующие выработки. Новосибирск : ВО «Наука», 1992. 133 с.

6. Будак Б. М., Самарский А. А., Тихонов А. Н. Сборник задач по математической физике. М. : Физматлит, 2003. 688 с.

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

8. Лаврентьев М. А., Шабат Б. В. Методы теории функций комплексного переменного. М. : Наука, 1987. 688 с.

9. Бицадзе А. В., Калиниченко Д. Ф. Сборник задач по уравнениям математической физики. М. : Наука, 1985. 312 с.

10. Коренев Б. Г. Введение в теорию бесселевых функций. М. : Наука, 1971. 288 с.

11. Никифоров А. Ф., Уваров В. Б. Специальные функции математической физики. М. : Наука, 1984. 344 с.

12. Попов Б. А., Теслер Г. С. Вычисление функций на ЭВМ. Киев : Наук. думка, 1984. 600 с.

УДК 539.374

ВЛИЯНИЕ ПРОСКАЛЬЗЫВАНИЯ НА ВИСКОЗИМЕТРИЧЕСКОЕ ТЕЧЕНИЕ УПРУГОВЯЗКОПЛАСТИЧЕСКОГО МАТЕРИАЛА МЕЖДУ ЖЕСТКИМИ КОАКСИАЛЬНЫМИ ЦИЛИНДРАМИ

А. С. Устинова

Институт автоматики и процессов управления Дальневосточного отделения РАН, Владивосток, лаборатория механики деформируемого твердого тела E-mail: ustinova@iacp.dvo.ru

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

Ключевые слова: упругость, вязкопластичность, большие деформации, остаточные напряжения.

Influence of Slippingon Viscosimetric Flow of a Elastoviscoplas-tic Material Between Rigid Coaxial Cylinders

A. S. Ustinova

Institute of Automation and Control Processes,

Far-Eastern Branch of RAS, Vladivostok,

Laboratory of a Deformable Solid Mechanics E-mail: ustinova@iacp.dvo.ru

The viscoplastic flow of an incompressible elastoviscoplastic material between two rigid coaxial cylindrical surfaces is considered when slipping of a material is possible at one of them. The solution is constructed using the model of large elastoviscoplastic deformations, Reversible deformation, development and braking of a viscoplastic flow, unloading and deformation under rotation in the opposite direction are considered. Laws of movement of elastic-plastic boundaries are received.

Key words: elasticity, viscoplasticity, large deformations, residual stresses.

ВВЕДЕНИЕ

Вискозиметрические и прямолинейные течения неньютоновских материалов в рамках модели Шведова - Бингама рассматривались достаточно подробно [1-7]. Упругие свойства деформируемых материалов при этом не учитывались. В настоящее время представляет интерес исследование эффектов,

© Устинова А. С., 2012

93