Получение аналитического решения уравнения теплового пограничного слоя на основе введения дополнительных граничных условий

УДК 532.526.4

ПОЛУЧЕНИЕ АНАЛИТИЧЕСКОГО РЕШЕНИЯ УРАВНЕНИЯ ТЕПЛОВОГО

ПОГРАНИЧНОГО СЛОЯ НА ОСНОВЕ ВВЕДЕНИЯ ДОПОЛНИТЕЛЬНЫХ

ГРАНИЧНЫХ УСЛОВИЙ

© 2010 Е. В. Стефанюк, И. В. Кудинов Самарский государственный технический университет

С использованием математической модели теплового пограничного слоя, представленной в виде интегрального уравнения Г. Н. Кружилина, на основе введения дополнительных граничных условий разработана методика получения приближённого аналитического решения исходного дифференциального уравнения (уравнения Польгаузена), позволяющая получать решения с заданной степенью точности. Проведено уточнение критериальной зависимости для определения коэффициента теплоотдачи на границе жидкость-стенка, уточнены формулы для определения теплового пограничного слоя. Проведён анализ распределения изотерм в пределах теплового пограничного слоя, исследованы скорости движения изотерм по поперечной координате в зависимости от величины продольной переменной.

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

Известно, что теплообмен между жидкостью и стенкой в значительной степени определяется процессами, протекающими в пограничном слое вблизи стенки. Интенсивность теплообмена характеризуется величиной коэффициентов теплоотдачи на границе жидкость-стенка. Для их определения необходимо иметь аналитическое решение задачи о распределении температуры в тепловом пограничном слое. Точные аналитические решения подобных задач в настоящее время не получены. Приближённые аналитические решения через решение интегральных уравнений Кармана и Польгаузена найдены лишь в первом приближении - их отличие от решений, полученных численным

интегрированием исходных дифференциальных уравнений, находится в пределах 8 - 10 %.

При обтекании тела потоком жидкости, имеющей скорость и (скорость невозмущённого потока), вблизи

поверхности тела образуется слой ¿(х) (гидродинамический пограничный слой), в пределах которого скорость течения изменяется от нуля на стенке до скорости невозмущённого потока (рис. 1).

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

среды изменяется от 4т до температуры невозмущённого потока 1ср (ст < ¿ср), где 4т -температура стенки (рис. 2).

Рис. 1. Схема гидродинамического пограничного слоя

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

Теплообмен в тепловом пограничном слое (рис. 2) описывается уравнением энергии (известным под названием уравнения Польгаузена [1] - [3]) и соответствующими граничными условиями

п

х у — ^

йг ду 5у2

(0 < у < 4(х); 0 < х <¥); *(х,0)_ *ст; *(х, 4)_ *Ср;

5*(х, 4) / ду _ 0,

иу -

координатным

п -

составляющие осям; х, у -коэффициент

где * - температура; их, скорости по координаты; температуропроводности.

Соотношения (3), (4) представляют условия сопряжения прогретой и непрогретой зон. Условие (3) означает, что температура на границе теплового пограничного слоя (на границе фронта температурного возмущения) равна температуре невозмущённого потока *ср. Согласно условию (4) тепловой поток не распространяется за пределы пограничного слоя. Математическое доказательство условий (3), (4) дано в [7].

Для получения ещё одного граничного условия запишем уравнение

применительно к точке у = 0. Так как в этом случае их

52*(х,0) _ 0

ду2 '

пределах толщины теплового пограничного слоя А(х):

/и _ п/

д*

дх

ду

ду2

4у.

Определив интегралы в (6), после некоторых преобразований получим следующее интегральное уравнение, впервые полученное Польгаузеном [1]:

й /иX [*ср - *

йх

_ п-

ду

Уравнение (7) можно также получить из рассмотрения балансовых соотношений по тепловым потокам в пределах теплового пограничного слоя [3].

Введём избыточную температуру по соотношению Т = * - *ст. Тогда Тср = *ср - *ст. Интегральное уравнение (7) и граничные условия (2) - (4) для избыточной температуры примут вид

4

х [тср - Т(x, у)] Ф

й

йх / х1ср Т(х,0)_ 0; Т(х, 4)_ ТСр; дТ(х,4) _ 0;

5у '

д2Т(х,0) _ 0

ду2 '

_ п дТ(х,0) ;

ду '

Найдём решение задачи (8) - (12), используя формулы для определения скорости их в пределах гидродинамического пограничного слоя [4]. Формула для их в первом приближении имеет вид

их _ и

Получение аналитического решения задачи (1) - (5) представляет значительные математические трудности [1, 3]. Поэтому потребуем, чтобы искомое решение удовлетворяло не уравнению (1), а некоторому осредненному уравнению (1), проинтегрированному по переменной у в

3 у -1 [ у 2 8 21 8

где и - скорость невозмущённого набегающего потока; ¿(х) - толщина гидродинамического пограничного слоя, определяемая по формуле

4

4

3

S(x) = 4,64л/vx / u ,

где v - коэффициент кинематической вязкости.

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

интенсивностью соответствующего процесса переноса. Так как характеристикой интенсивности переноса импульса является кинематический коэффициент вязкости, а теплоты - коэффициент

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

Ввиду того, что толщины гидродинамического и теплового

пограничных слоёв должны подчиняться условию А(х)<^(х) [1] - [3], то величина критерия Прандтля должна удовлетворять условию Pr > 1. Это условие приближённо выполняется для газов (Рг ~ 0,75) и для неэлектропроводных жидкостей (Рг > 1) и не выполняется для жидких металлов ввиду высокого значения коэффициента

температуропроводности (10-3 < Pr < 10-2).

Решение задачи (8) - (12) принимается в виде

п

T (x, y ) = £ ak (D У ,

k=0

t _зy iœys

T

ср

2 D 2 è D

где А(х) - толщина теплового пограничного слоя.

Подставив (13) и (16) в интегральное уравнение (8), найдём

3 (D 2ôl S

= ad7(x,0) dy

= a

Дифференцируя (16) по переменной y, применительно к y = 0 будем иметь

dT(x,0) = 3 ТСр dy 2 D

Если принять А < ¿, то вторым членом в левой части уравнения (17) можно пренебречь. Тогда соотношение (17) с учётом (18) примет вид

d_

dx

и—(p2ô) = 10 —,

где в = Ы0.

Ввиду независимости величины в от координаты x, de/dx = 0. Отсюда получаем

1 из? dd

— uB3ô— = a . 10 dx

D(x ) =

Подставив (14) в (20), найдём 4,64x

л/rÎ-Vpr

где а^(А) - неизвестные коэффициенты, определяемые из граничных условий (9) -(12). После определения коэффициентов а^(£ = 0,3) соотношение (15) принимает вид

где Явх = и XIV, Рг = vla.

Соотношения (16), (21) определяют решение задачи (8) - (12) в первом приближении. Результаты расчётов относительных избыточных температур 0 = Т1ТСр = ( - I (4р - 4т) по формуле (16) в сравнении с точным решением [1] представлены на рис. 3. Их анализ приводит к заключению о том, что расхождение полученных по формуле (16) значений температур от точных их значений находится в пределах 5 - 11 %. Причём максимальное расхождение наблюдается вблизи верхней границы теплового пограничного слоя.

4

3

1,0 0

0,9

0,8

0,7

0,6

0,5

0,4

0,3

0,2

0,1

,0

Рис. 3. Распределение безразмерных температур & = Т/Тср в зависимости от безразмерной координаты ц = Ул/ и/пх . 1, 2, 3 - соответственно первое, второе и третье приближения;

4 - точное решение [1] Соотношение (16) благодаря

полиномиальной зависимости температуры от координаты У позволяет построить линии изотерм в пределах толщины теплового пограничного слоя в координатах "у - х" (рис. 4). Задавая любые постоянные значения безразмерной температуры Т/Тср = &, для различных значений координаты х можно найти такие У, которые удовлетворяют соотношению (16).

Анализ распределения изотерм позволяет заключить, что все они (0 < & < 1) возникают на поверхности стенки в точке х = 0, у = 0. Изотерма нулевой температуры & = 0 совпадает с осью х. Изотерма единичной температуры & = 1 совпадает с линией теплового пограничного слоя. Отмечается сгущение изотерм вблизи стенки и их разрежение вблизи границы, о том, что максимальную скорость отделяющей возмущённый поток от перемещения имеет единичная изотерма. невозмущённого. Скорость перемещения нулевой изотермы

На основе графиков рис. 4 по равна нулю. Все изотермы возникают на соотношению Ж = Ду/Дх определяются поверхности стенки в точке х = 0, у = 0, имея скорости перемещения изотерм по при этом бесконечно большие начальные координате у в зависимости от координаты х скорости. Затем по мере продвижения (рис. 5). Их анализ позволяет сделать вывод изотерм по координате у в зависимости от

координаты х их скорости существенно

0,0015 У, м 0,0010

5-10"'

0 = 1,0 '0,97-' 0,9  0,8 .

0,6 . 0,5 0,4

0,3. 0,2

0,1

- 0,001

0 0,2 0,4 0,6 0,8

х, м

Рис. 4. График распределения изотерм & = Т/Тср в тепловом пограничном слое. Рг = у/а = 1, и = 5 м/с

уменьшаются с последующей стабилизацией изменения по закону, близкому к линейному. Отметим, что наибольший градиент скорости имеют изотермы малого потенциала на относительно небольшом расстоянии по координате х.

0,0030

ш

0,0025 0,0020 0,0015

0,0010 5-10"4

0,1

0,2 а

0,3

0,4

0,004 №

0,003

0,002

0,001

вновь возникающих при этом неизвестных коэффициентов а^(А) будем привлекать дополнительные граничные условия. Принцип их получения заключается в следующем [4] - [6]. Для получения первых из них уравнение (1) применяется в точках у = 0 и у = А. Именно таким путём было получено дополнительное граничное условие (5). Для получения второго дополнительного граничного условия применим уравнение (1) к точке у = А(х), которое с учётом (4) будет иметь вид

д ^х, Л)_ а д2t

д х

Ьх

д у2

Дифференцируя (2) по переменной х и сравнивая полученное соотношение с (22), находим второе дополнительное граничное условие (первым из них является соотношение (5)) вида

ду2

_ 0.

\\ X \о,4 40-35 "" 0,45

\о ХД25 ЧД2 3

\ \Д1 ,06 5 """

Чаи чО.ОШ^'

Для получения третьего

дополнительного граничного условия продифференцируем уравнение (1) по переменной у и применим полученное соотношение для точки у = А(х):

дих дt(х, Л) + ^ д2^х, Л) + дЬу дt(х, Л) +

ду дх

дхду ду ду

+и,

д2t

ду2

_а

ду3

Соотношение (24) с учётом (4), (5), а также того, что (дих / ду)у_8(х) _ 0, примет

вид

х. м

О 2 4 6 х-103, м 10

б

Рис. 5. Графики изменения скоростей движения изотерм Ж = Ау/Ах по координате у в зависимости от координаты х

Для повышения точности решения задачи (1) - (5) необходимо увеличивать степень полинома (15). Для определения

д2t(х, Л) _ а_ дХх,Л) дхду и х ду3

Дифференцируя (4) по переменной х и сравнивая полученное соотношение с (25), находим третье дополнительное граничное условие:

д3/

~~ду3

= 0.

1 2 А

Т 3 Т

, ср -3 х ср

а4 = 5 —г; а5 =---г.

А 3

А

45

2 А 5

Подставляя найденные значения коэффициентов ак в (15), находим

^ 2 У - •(А)'-Ш-Й <27.

Формула для распределения скорости во втором приближении имеет вид [4]:

их = и

где

2 У 5

д(х ) = 7,035^л/ПхТй.

Подставив (27) и (28) в интегральное уравнение (8) и определив интегралы относительно неизвестной функции Д(х), получаем следующее обыкновенное дифференциальное уравнение:

и ё

(

5544ёх

.А2

„А4

,А5

.А

6^

660--165-г+77— -12—

5 53 54 55

=5 а. (30) 2А

Ввиду того, что Д<^, то всеми членами в левой части соотношения (30), кроме первого, можно пренебречь:

165 ё V

5а

(в2 ) = 1388 ёхуи ' 2 рд

Так как величина в = Д/^ не зависит от х, то последнее уравнение приводится к виду

Подставляя (15) при п = 6 в граничные условия (2) - (5), (23), (26), относительно

неизвестных коэффициентов ак (к = 0,5) будем иметь систему шести алгебраических линейных уравнений (отметим, что все дополнительные граничные условия для / и Т идентичны). Её решение имеет вид

5 Тс Тс

а„ = 0; а, = ——; а2 = 0; а3 = -5 ;

1

ёд

—ив3 д — = а. 21 ёх

А(х ) =

Подставляя (29) в (31), находим 6,661х

Соотношения (27), (32) представляют решение задачи (1) - (5) во втором приближении. Результаты расчётов безразмерных температур & = Т/Тср по формуле (27) показаны на рис. 3. Их анализ позволяет сделать вывод о том, что уточнение решения во втором приближении по сравнению с первым составляет около 3 % для 0,5 < щ < 3,0, а вблизи границы теплового пограничного слоя (п > 5,0) полученное во втором приближении решение практически совпадает с точным.

Соотношение (27) точно удовлетворяет всем основным (2) - (4) и дополнительным (5), (23), (26) граничным условиям, а также интегральному уравнению (8). По сравнению с первым приближением в данном случае наблюдается более точное выполнение уравнения (1) ввиду выполнения

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

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

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

соотношение для точки = 0:

дих д/(х,0) + и д2^х,0) + диУ д/(х,0) +

дУ дх

дхд д д

+ V У

д2/

дУ2

= а-

дУ3

Продифференцируем условие (2) по переменным х и У:

д/ (х,0) = 0; д/ (х,0) = 0.

дх ' дУ

Соотношение (33) с учётом Vх|У=0 = VУ | 0 = 0 приводится к виду

д3/(х,0) = 0

дУ3

(34) (34) и

Соотношение (35) представляет первое дополнительное граничное условие третьего приближения.

Аналогично, путём двух" и трехкратного дифференцирования

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

д 4/(х, А ) = 0; д5/(х, А) = 0.

дУ4

дУ

Граничные условия (2) - (5), (23), (26), (35), (36) позволяют найти уже девять неизвестных коэффициентов ак(Д) (к = 0,8) ряда (15). Подставив (15) в перечисленные граничные условия, относительно неизвестных коэффициентов ак(Д) получим систему девяти алгебраических линейных уравнений. Её решение имеет вид:

8 Тср

а0 =0; ° =3а ;

Тср Тср

ср

а =0; а =0; а4 =——;

70Тр

3 А4

Тср

—ср _ ср а =56А^' а =-56А^' ^ 3 А7"4 "А

80 ^ср ^ср

а ~т; а8 =-5^г.

Подставляя найденные коэффициенты ак(Д) в (15), находим

М =- 70 (у I4+56 (у'

ср

3 А 3 (А

граничное -56

у

- 80 (1 3 (А.

- 5

.У

Формула для распределения скорости их в третьем приближении имеет вид [4]:

и = и

8У

3 5

—561

70 ( У

3 (5

\6 У 1 ,

+ 56

80 ( у

3 и.

- 5

.У

где

д(х) = 7,7931л/пх /V .

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

пограничных слоев, а также всех основных и дополнительных граничных условий.

Подставляя (37) в интегральное уравнение (8), относительно неизвестной функции Д(х) приходим к следующему обыкновенному дифференциальному

уравнению:

ё_ ёх

и

196911

19448

А1

3672

А

8а.

3 А

А9 А5 +9^—5- + 6644-7-

ё8 ё4

А6 А8 -7344—- 935^=-

ё5 ё7

Так как Д<ё, то всеми членами в левой части уравнения (40), кроме первого, можно пренебречь.

19448 196911

V

ё_ ёх

8 а 3 р5

Учитывая, что в = А/д не зависит от х, находим

1 О3о й8

—ир 8— _ а . 27 йх

Подставив в (42) соотношение (39), получим

А(х)_

7,491х

д4t(х,0) _ 0 ; д6t(х, Л) _ 0 ; д7t(х, Л) _ 0

ду4

ду

ду7

После определения коэффициентов ак(А) (к соотношение (15) примет вид

неизвестных

0,11)

5 о3с й8

-ир 8 — _ а.

156 йх

Л(х)_

Подставляя (46) в (47), находим 8,017 х

Соотношения (37), (43) представляют решение задачи (1) - (5) в третьем приближении. Результаты расчётов по формуле (37) в сравнении с точным решением [2] приведены на рис. 3. Их анализ позволяет заключить, что отклонение полученных по формуле (37) безразмерных температур от их точных значений не превышает 2 %.

Дополнительные граничные условия, необходимые для получения решения задачи (1) - (5) в четвёртом приближении, имеют вид

Соотношения (45), (48) представляют решение задачи (1) - (5) в четвёртом приближении. Результаты расчётов по формуле (45) показывают, что расхождение с точным решением не превышает 0,5 %. Сравнение результатов расчётов

безразмерных температур в первом и четвёртом приближениях показано на рис. 6.

На рис. 7 даны результаты расчётов по формуле (45) для различных значений критерия Прандтля (Рг = у/а = 0,6; 1,0; 3,0; 15,0) в сравнении с точным решением [1], [2]. Как видно из рисунка, для Рг = 1; 3; 15 полученные по формуле (45) значения безразмерных температур 0 = Т/Тср = ( -4т)/(4р - 4т) практически совпадают с их точными значениями. Расхождение результатов, составляющее около 3 %, наблюдается для Рг = 0,6, т. е для газов.

На основе полученных выше зависимостей для скорости и температуры в пределах соответствующих пограничных слоёв и использования дифференциального уравнения конвективной теплоотдачи

Т 11 у 231 у5 у6 у7 у8 - у у + 462 - 825 Аг + 825Аг-

Т

ср

4 Л 2 Л 5

Л 6

Л 7

>925 у9 + ,54¿0-21 у"

4Л

9

Л

10

Л

11

Л 8

4(§)у _ О- ^ ) •

можно определить

теплоотдачи

Формула для скоростей их в четвёртом приближении записывается в виде, аналогичном формуле (45), где д(х) определяется следующим образом:

а _

1 дТ

(тср - Тст) ду '

коэффициент

8(х) _ 8,2572л/йхТи .

где X - коэффициент теплопроводности.

Если использовать решение в первом приближении (16), то для коэффициента теплоотдачи получим формулу

Обыкновенное дифференциальное уравнение относительно А(х) имеет вид

а _ 1,5

4

Л

Для решения (45) в четвёртом приближении будем иметь

^ ггп 1

а = 2,75-

А(х )•

Подставив (48) в (52), получим

Ш х = 0,3430^^^.

Отсюда следует, что расхождение коэффициентов Нуссельта в критериальных

Подставляя (21) в (5l), определим уравнениях теплоотдачи в первом

и

критерии

Шх = — = 0,323^7^^7^Рг .

1,0

0

0,8

четвёртом приближениях составляет 1,98 %.

0,7 0,6

0,5 0,4

0,3 0,2 0,1

Д(0,0001) ___.

У /

1х = 0 0001 м />0,001 уу У 0,009

// I / д(000]) / У / У 0,004

// Д

II // //

II // // у у \

-—и- 1 --И-- // / - / у У

а // /у У У у'

1 // / / /V У

1 // / / / / / / / /

/ / 7 / /У У

// //

// // // х-

( А / / /у /

// // /о-

/ / /У

/ /\

/ /У

ШУ

0

0,45

0,9

1,35

1,8

2,25

2,7

3,15

3,6 у-104, м

4,5

Рис. 6. Изменение безразмерных температур 0 = Т/Тср в пределах теплового пограничного

слоя.------по формуле (16) (первое приближение);--по формуле (45) (четвёртое

приближение)

1,0 ©

О 1 2 3 4 5 6^7

Рис. 7. Распределение безразмерных температур 0 = Т/Тср в зависимости от безразмерной

координаты ] = ул/ и/пх для различных значений критерия Прандтля (Рг = у/а).--по

формуле (45); о - точное решение [1]

По результатам проведённых исследований можно сделать следующие выводы.

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

2. Дополнительные граничные условия находятся из дифференциального уравнения Польгаузена путём выполнения этого уравнения и производных от него в граничной точке (у = 0) и на фронте температурного возмущения (на границе пограничного слоя). Так как диапазон изменения фронта возмущения охватывает весь диапазон изменения поперечной координаты (0 < у < А(х)), то, следовательно, для всех значений переменной х, которым

соответствуют значения переменной у, обозначающие линию пограничного слоя, уравнение Польгаузена выполняется точно, в чём можно убедиться непосредственной подстановкой соотношений (16), (27), (37), (45) в (7). Расчёты показывают, что с увеличением числа приближений (степени аппроксимирующего полинома) решение всякий раз уточняется.

3. Выполнен анализ распределения изотерм в пределах теплового пограничного слоя и показано их сгущение вблизи поверхности стенки на весьма малом расстоянии по координате у. Проведён анализ изменения скоростей движения изотерм по поперечной координате у в зависимости от величины продольной координаты х и показано, что наибольшее изменение скорости их движения наблюдается в узком пристеночном слое и на незначительном отрезке вдоль координаты х.

4. На основе полученных уточнённых аналитических решений уравнения Польгаузена уточнены формулы для определения коэффициентов теплоотдачи на границе "жидкость-стенка" на величину 1,98 %. Уточнены формулы для определения толщины теплового пограничного слоя.

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

1. Шлихтинг, Г. Теория пограничного слоя [Текст] / Г. Шлихтинг. - М.: Наука, 1969.

2. Юдаев, Б. Н. Теплопередача [Текст]: учебник для вузов / Б. Н. Юдаев. - М.: Высшая школа, 1981.

3. Болгарский, А. В. Термодинамика и теплопередача [Текст] / А. В. Болгарский, Г. А. Мухачев, В. К. Щукин. - М.: Высшая школа, 1975.

4. Стефанюк, Е. В. Получение аналитического решения уравнений гидродинамического пограничного слоя на основе введения дополнительных граничных условий [Текст] / Е. В. Стефанюк, Б. В. Аверин, И. В. Кудинов // Известия Самарского научного центра Российской

академии наук. Специальный выпуск: "Актуальные вопросы тепло- и массообмена, энергоэффективность, исследование

вихревых закрученных потоков". - Самара. Самарский научный центр РАН, 2008. -С. 39-46.

5. Кудинов, В. А. Теплопроводность и термоупругость в многослойных конструкциях [Текст] / В. А. Кудинов, Б. В. Аверин, Е. В. Стефанюк. - М.: Высшая школа, 2008.

6. Кудинов, В. А. Задачи теплопроводности на основе определения фронта температурного возмущения [Текст] / В. А. Кудинов, Е. В. Стефанюк // Известия АН. Энергетика. - 2008. - №4. - С. 122-138.

References

1. Schliechting, G. Boundary layer theory / G. Schliechting. - Moscow: Nauka (Science), 1969.

2. Yudayev, B. N. Heat transfer: Textbook for higher school / B. N. Yudayev. - Moscow: Vysshaya shkola (Higher school), 1981.

3. Bolgarsky, A. V. Thermodynamics and heat transfer / A. V. Bolgarsky, G. A. Mukhatchyov, V. K. Shchukin. - Moscow: Vysshaya shkola (Higher school), 1975.

4. Stefanyuk, Ye. V. Obtaining an analytical solution of hydrodynamic boundary layer equations on the basis of introducing additional boundary conditions / Ye. V. Stefanyuk, B. V. Averin, I. V. Kudinov // News of the Samara Scientific Centre of the Russian Academy of Sciences. Special issue: "Topical questions of

heat and mass exchange, power efficiency, research of vortex swirled flows". - Samara. Samara Scientific Centre of the Russian Academy of Sciences, 2008. - PP. 39-46.

5. Kudinov, V. A. Heat conduction and thermoelasticity in multi-layer structures / V. A. Kudinov, B. V. Averin, Ye. V. Stefanyuk. -Moscow: Vysshaya shkola (Higher school), 2008.

6. Kudinov, V. A. Heat conductivity problems on the basis of defining the temperature perturbation front / V. A. Kudinov, Ye. V. Stefanyuk // Izvestiya (news) of the Academy of Sciences. Power engineering. -2008. - No. 4. - PP. 122-138.

OBTAINING THE ANALYTICAL SOLUTION OF THE HEAT BOUNDARY LAYER EQUATION BASED ON INTRODUCING ADDITIONAL BOUNDARY CONDITIONS

© 2010 Ye. V. Stefanyuk, I. V. Kudinov Samara State Aerospace University

A method of obtaining an approximate analytical solution of the initial differential equation (Polgauzen equation) has been developed using a mathematical model of the heat boundary layer presented in the form of Kruzhilin integral equation on the basis of introducing additional boundary conditions. The method developed makes it possible to obtain solutions with the specified degree of accuracy. Criterial dependence for defining the convective heat exchange coefficient at the fluid-wall boundary has been refined, the formulae for defining the heat boundary layer have been revised. Isotherm distribution within the heat boundary layer has been analysed, the velocity of isotherm motion along the transverse coordinate depending on the value of the longitudinal variable has been investigated.

Heat boundary layer, temperature perturbation front, additional boundary conditions, integral method of heat balance, isotherm, isotherm motion velocities.

Информация об авторах

Стефанюк Екатерина Васильевна, кандидат технических наук, старший преподаватель кафедры "Теоретические основы теплотехники и гидромеханика" Самарского государственного технического университета, stef-kate@yandex. ru. Область научных интересов: аналитические методы решения краевых задач теплопроводности, тепломассопереноса и термоупругости.

Кудинов Игорь Васильевич, аспирант кафедры "Прикладная математика и информатика" Самарского государственного технического университета, igor koudinov@mail.ru. Область научных интересов: аналитические методы решения краевых задач теплопроводности, тепломассопереноса и термоупругости.

Stefanyuk Yekaterina Vasilyevna, candidate of technical science, senior lecturer of the department "Theoretical foundations of heat engineering and hydromechanics", Samara State Aerospace University, stef-kate@yandex.ru. Area of research: analytical methods of solving boundary problems of heat conduction, heat and mass transfer and thermoelasticity.

Kudinov Igor Vasilyevitch, post-graduate student of the department "Applied mathematics and information science", Samara State Aerospace University, igor_koudinov@mail.ru. Area of research: analytical methods of solving boundary problems of heat conduction, heat and mass transfer and thermoelasticity.