Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Жаров Г.Г. Судовые высокотемпературные газотурбинные установки

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.85 Mб
Скачать

его частные решения, которые удовлетворяют граничным условиям. Затем составляется ряд этих решений вида:

 

л

 

t =

- j - C%t% + • • • + C,J-,і = 2j CJn-

(138)

 

;i = l

 

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

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

из которых есть функция координат,

 

другая — времени.

Тогда

выражение для частного

решения

примет

вид

 

 

 

 

t

=

cf

(х, у,

г) ср (т),

 

(139)

где

с—произвольная

постоянная;

 

 

 

 

 

/ (х, у, z) — функция

от

координат;

 

 

 

 

 

 

ср (т) — функция

времени.

 

 

 

 

 

 

 

Подставим

уравнение

(139) в

уравнение

(137):

 

 

 

ф' М f (х, у,

z)

=

 

аср (т) Y2 f (х,

у,

z)

(140)

и,

произведя

разделение

переменных,

получим

 

 

 

 

 

Щ

 

 

у

'

г)

.

(141)

'

 

 

 

Ф (т)

 

/ (х, у,

 

к

Левая часть равенства (141) зависит только от времени, а правая — только от положения точки в пространстве.

Поскольку равенство (141) имеет место при любых значениях координат и времени, то каждая из его частей равна какой-то по­ стоянной величине. Обозначим эту величину через А. Тогда запишем

Щ.=А;

a t ! [ X

'

У

\ г )

=А.

(142)

ФМ

f{x,

 

У,

г)

 

у

>

 

 

 

 

 

 

 

Проинтегрировав первое выражение,

 

получим

 

 

 

Ф (-с) =

еА\

 

 

(143)

При этом, если тело стремится к температурному равновесию, то величина А будет отрицательной. Если температура тела увеличи­ вается, то величина А будет положительной. Если температура тела есть периодическая функция времени, то величина А будет мнимой величиной. В практике чаще всего встречается первый случай: А «£і 0. Рассмотрим его. Так как А выбирают из физических соображений, то можно положить

А = —а/с3 ,

(144)

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

 

к— постоянная, определяемая из граничных

условий.

Подставляя значение

А

в

(142),

находим

 

 

 

Ф (т) е~ак"-х;

(145)

V*f(x,

у,

г)

+K*f(x,

у, г) = 0 .

(146)

Если известно решение уравнения (146), то частное решение уравне­

ния теплопроводности имеет вид

 

і = ce-aK°-xf {х, у, г).

(147)

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

Общее решение можно представить как сумму частных решений по выражению (138). Решение объемной задачи таким методом яв­ ляется очень затруднительным. Чаще всего этим методом решают одномерные задачи, связанные с нахождением температурного поля в неограниченной пластине, цилиндре и шаре. Для охлаждаемых узлов такой метод применяют при определении температурного поля оболочковых лопаток. Приведем решение рассмотренным методом (методом разделения переменных) для оболочковой охлаждаемой лопатки. Считаем оболочку охлаждаемой лопатки как неограничен­

ную пластину. Тогда

уравнение теплопроводности примет вид

 

 

 

 

 

 

 

 

dt (х,

т)

_

а

дЧ

(я,

т)

 

 

 

^ 1 4 8 ^

 

 

 

 

 

 

 

дх

 

 

 

 

 

дх

 

 

 

 

 

 

Частное решение

этого уравнения

можно

представить

 

 

 

 

 

 

 

 

/

=

cf

(х) ср (т).

 

 

 

 

(149)

Подставляя

(149)

в

(148),

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф

=

а

Щ -

 

= -ак\

 

 

 

(150)

 

 

 

 

 

 

Ф(т)

 

-

Д А - )

 

 

 

 

 

 

 

Проинтегрировав

левую часть

уравнения

(150),

найдем

 

 

 

 

 

 

 

 

 

 

ср (т )

=

е - я к ' т .

 

 

 

 

(151)

Дифференциальное

уравнение

для

функции f (х)

имеет вид

 

 

 

 

 

 

 

 

Г

(х)

+

K*f (х)

=

0.

 

 

 

(152)

Частными

решениями

уравнения

(152)

являются

функции

sin

кх

и cos кх.

Эти решения — линейно-независимые.

Тогда общие

реше­

ния, уравнения

представим

 

 

 

 

 

 

 

 

 

 

 

 

/

(х)

=

cfx (х)

+ D / 2 (х)

=

с sin кх + D

cos

кх.

(153)

Частное

решение

уравнения

теплопроводности

выразим,

как

 

 

 

і (х,

х) — с sin хе-ак'х

 

-V D cos кхе~ак'х.

 

(154)

Общее

решение

уравнения

теплопроводности

запишем

 

 

 

00

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

t=

УІ

сп

sin кпх

ехр (—ак1%) -\-

S

Dm cos кт

 

ехр (—акт х),

(155)

n = l

т = 1

где с и D — произвольные постоянные, определяемые из начальных условий;

к — определяется из граничных условий.

Заметим, что даже в самом простом случае решение уравнения теплопроводности методом разделения переменных приводит к зна­ чительным трудностям. Поэтому им пользуются только в особых слу­ чаях. Существует и ряд других классических методов решения урав­ нения теплопроводности. Однако их применяют сравнительно редко. Из них заслуживает внимания метод источников, сущность которого состоит в том, что любой процесс распространения тепла в теле пред­

ставляют как

совокупность

процессов

выравнивания температуры

от множества

элементарных

источников

тепла, распределенных как

в пространстве, так и во времени. Задача в основном сводится к пра­

вильному выбору источников

и их

распределению

[42].

В целом следует отметить,

что

решение задач

теплопроводности

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

§ 46. Решение

дифференциального

уравнения

теплопроводности

методами

интегрального преобразования

К методам решения уравнения теплопровод­

ности путем

интегрального преобразования относятся операцион­

ные методы

и конечные интегральные преобразования.

При решении задач с внутренними стоками тепла операционные методы имеют большие преимущества по сравнению с методом разде­ ления переменных. Для охлаждаемых узлов газовой турбины, где происходит отвод тепла за счет внутреннего его стока, такие задачи являются основными. Разработка операционного метода и его обосно­ вание дано в работах [14, 37].

Операционный метод широко применяют в разных областях тех­ ники и рассматривают как самостоятельный для решения уравнения математической физики. По своей сущности он равнозначен инте­ гральному преобразованию Лапласа, который используют в теории функций комплексного переменного. В основу данного метода поло­ жено не изучение самой функции, а ее видоизменения. При этом рас­ сматриваемую функцию называют оригиналом, а ее видоизменение — изображением. Целесообразность перехода от исследуемой функции (оригинала) к ее видоизменению (изображению) заключается в том, что решается не дифференциальное уравнение для оригинала функ­ ции, а алгебраическое уравнение для ее изображения. При решении дифференциального уравнения преобразуют функцию в ее изобра­

жение (алгебраическое уравнение),

решают его, а затем переходят

от

изображения

обратно к функции. Преобразование от функции

к

изображению

осуществляется

путем умножения функции на

экспотенциальную функцию и интегрирования ее в определенных пределах.

Поскольку преобразование Лапласа является интегральным и, следовательно, обладает свойством операторов, то вместо дифферен­ циального уравнения для оригинала функции получают алгебраи­ ческое для нее изображение. Если имеется функция / (t) и ее изобра­ жение F (р), то зависимость между ними можно выразить так:

со

 

 

F(p) = l f(i)e~ptdt

= L[f(t)].

(156)

о

Для существования изображения оригинала необходимо, чтобы

интеграл

сходился; р—комплексное

число, причем

вещественная

часть в

нем — положительная. С

целью нахождения

изображения

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

жение представляет собой

 

 

 

.

F(k) = \ке-р* dt=4-

 

e~pt

f

(157)

о

 

 

 

 

при р > 0 .

 

 

Если оригинал мы имеем в виде к,

t,

то изображение

можно пред­

ставить

 

 

 

 

со

 

 

 

 

F (кі) = J Kte~pt

dt

= -~.

 

(158)

Получив таким образом изображение и решив его, переходим обратно к оригиналу. Оригинал по его изображению находят, как правило, по формуле обращения (обратное преобразование)

°~-Ь'со

?W = ikr

I F(p)ePtdP-

(159)

 

а—і со

 

Интегрирование происходит в комплексной плоскости

вдоль прямой

а = const, параллельной мнимой оси. Как правило, обратное пре­ образование осуществляют, используя таблицы и не прибегая к кон­ турному интегрированию. Вместо приведенной формулы можно использовать и другую

f (t) = lim

(—1)"

п

<16°)

 

(тГ"<р>(т)]-

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

13 Г. Г. Жаров

193

Дифференциальное уравнение температуры по высоте охлаждае­ мой лопатки теплоотводом в диск в общем виде запишем

 

 

 

 

 

£ + Л - 0 .

 

 

 

 

 

 

('61)

Пусть при х = О

 

 

 

 

tx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ (0) =

 

=

const;

 

 

 

 

 

 

 

 

 

 

 

f

(0) =

 

А

=

const.

 

 

 

 

 

Дифференциальное

уравнение представим

 

 

 

 

 

 

 

 

 

 

t" (х) — кН (х) =

0.

 

 

 

 

(162)

Применим

прямое

преобразование

Лапласа

 

 

 

 

 

 

 

 

 

t(p) = \

e~pxt(x)dx

 

=

L [t(x)\,

 

 

(163)

т. е.

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

U" (х)] — к2Ь

(*)] =

0.

 

 

(164)

Используя

основные

свойства

преобразования

Лапласа,

находим

 

 

р2 t (р) — tlP

А кЧ (р) = 0.

 

 

(165)

Последнее уравнение является алгебраическим, поэтому решаем

его, считая р

простым

числом:

 

 

 

 

 

 

 

 

 

 

 

 

t

=

Ч р +

^

=tj_

,

р

0

+

А—-0

к

 

.

(166)

 

р

 

р - — к-

 

1

р 2

— к 2

1

( р 2

— к 2 ) к

 

v

'

Пользуясь

обратным

преобразованием

Лапласа,

определяем

 

 

 

 

 

 

 

 

 

А

 

 

 

(167)

 

 

 

 

t (х) = tt ch кх

 

r

sli /ех.

 

 

 

Это и будет

температура

по

высоте

охлаждаемой

лопатки тепло-

отводом в

диск.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

метод

преобразований

Лапласа

при

решении

уравнений теплопроводности имеет ряд преимуществ перед методом

разделения переменных. Прежде

всего решение задачи однотипно

и не требует особого подхода в

каждом отдельном случае. Кроме

того, интегральные преобразования Лапласа позволяют в равной степени решать задачи при граничных условиях первого, второго, третьего и четвертого родов, а также быстро и легко решать задачи с простыми начальными условиями (по одной координате). Исполь­ зование этого метода в инженерной практике упрощается наличием таблиц изображений, в результате чего задача сводится к решению только алгебраического уравнения. Однако решение задач таким методом затруднено при задании начальных условий в виде функции пространственных координат. В таком случае прибегают к использо­ ванию преобразований Фурье и Хенкеля [42]. Выполняя эти пре­ образования, следует обращать большее внимание на сходимость интегралов, так как условия сходимости становятся более жесткими, чем условия сходимости в интегралах Лапласа.

Все рассмотренные преобразования можно применить для расчета тел полуограниченной протяженности. При решении задач с конеч­ ной областью изменения переменных создан ряд методов конечных интегральных преобразований. Эти методы мы специально рассма­ тривать не будем. С ними можно ознакомиться в специальной лите­ ратуре [14].

§ 47. Решение дифференциального уравнения теплопроводности приближенными методами

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

Метод конечных разностей или

метод сеток основан на замене про­ изводных их приближенным значе­ нием или замене непрерывной кривой

рассматриваемого

процесса

ломаной

и непрерывного во

времени

процесса

прерывным

с фиксированным

шагом.

Тем самым

дифференциальное

урав­

нение теплопроводности

сводится

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

Рис. 105. К замене частных производных конечно-разностными вы-

раженими.

где очередное значение температуры в точке является функцией предыдущего по времени и координатам. Применение метода конеч­ ных разностей при ручных расчетах было весьма затруднено. С по­ явлением ЭВМ этот метод нашел широкое распространение.

Метод конечных разностей обеспечивает большие преимущества при расчете нестационарных температурных полей, а также дает возможность учесть физические константы и коэффициенты тепло­ отдачи. Сущность метода сводится к замене производных через раз­

ностные отношения [42]. Представим, что нам дана функция t

=

f (х)

( р И С . 105). И з В е С Т Н О , ЧТО П р О И З В О Д Н у Ю ф у Н К Ц И Ю П р и

X =

Х(

можно

представить

 

 

 

 

 

 

 

 

 

t't = igat.

 

 

 

 

(І 68)

Выберем

на кривой

две соседние точки

А г _х ;

^_х )

и с

 

^•+ 1 ) таким

образом,

чтобы разности xt

— xt_± = xi+1

— xt

= h

были бы достаточно малы. Считая шаг h достаточно малым, мы имеем

13*

195

право приближенно принять равными и углы at = р\ = 7,., а сле­ довательно, рассматривать вместо касательной одну из секущих АВ или. ВС.

Тогда первая производная в точке і примет вид

или

л ~ to- v . —

^

^

^"-1

' ~ щ Y l _

АЕ

 

ft

или

 

 

 

 

 

 

Приближенное значение второй производной получаем, если заменим кривую на ломаную линию

/' ^_

_ j _

(

^1+1 ti

h — h-l

\

^1+1 %h -f- / і - !

 

,, K Q \

h ~

ft

\

ft

ft

) -

h?

[ i b J >

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

Пусть дан стержень длиной /. Распределение температуры для изолированного тонкого стержня описывается уравнением вида

Поскольку температура зависит от времени и координаты (от

двух переменных),

то обозначим

истинное

значение

температуры

в точке

через х

k,

где

і

— количество шагов

по

координате

х

до

рассматриваемой

точки,

a

k — количество шагов

по

координате

т

до рассматриваемой

точки. Обозначим: h — шаг

по

координате

х;

п — шаг по координате

т. Тогда

уравнение

(170) в

узле ik

заме­

няется

выражением

 

 

 

 

 

 

 

 

 

 

 

U, k+i-

ti, k

_j_

^ =

a ^ tt-i.

k - 2ti, k +

tt+1.

к

 

 

(171)

где £x и £ 2 — остаточные члены при шаге, не равном нулю.

Чем меньше шаг, тем меньше t,i и £ 2 и тем точнее замена диффе­ ренциального уравнения уравнением в конечных разностях. Преоб­

разуя уравнение

(171),

получаем

 

 

 

 

к ы = (1

Ж " )

* + I F

* +

k ^ h

~

( 1 7 2 )

Отбросив остаточный член h (а£2 t,i), вследствие его малости, найдем приближенное значение расчетной температуры в том же узле

к *+х = (1 + т )

к k + ^ " tt-i. k + tl+1. k ) .

(173)

По формуле (173) можно

определить температуру в

первом ряду

по х при т =

п, при известных краевых условиях т =

0.

Получен­

ные значения

температур в

первом ряду будут являться

краевыми

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

В случае решения нестационарных уравнений в частных произ­ водных параболического типа соотношение шагов, а также ошибка округления в численном решении определяет сходимость и устой­ чивость получаемых решений. Для некоторых соотношений шагов можно получить ряд упрощенных частных решений:

При

п = 1г213а

 

 

 

 

 

 

 

H.fc+i-

 

з

 

При /і =

h2/Qa

 

 

 

 

 

 

 

і

k+i

ti-vk + 4tj,k

+ U+i.k

 

 

4.

 

5

 

При n =

h2/\2a

 

 

 

 

 

 

 

і

_

h-i,k +

10^. k + fi+i, k

 

 

h, ft+i —

 

12

 

При

n =

h2/pa

 

 

 

 

 

 

 

j.

_ tj-v

k + (p — 2) ti, к + ti+i, к

 

 

h, k+i -

 

-

 

Формула

значительно

упрощается

при р = 2

 

Рассмотренный метод применяют при решении двух- и трехмер­

ных

задач. Однако при этом требуется значительное время (даже

на

ЭВМ). Например, для получения плоскостного стационарного

температурного поля охлаждаемой натрием по шести каналам ло­ патки требуется 3,5 ч рабочего времени ЭЦВМ М-220. Поэтому в на­ стоящее время разрабатываются новые разновидности метода конеч­ ных разностей, которые позволяют значительно сократить время расчета. К таким методам прежде всего следует отнести локальноодномерный метод переменных направлений.

Метод конечных разностей сложнее применить при рассмотрении тел более сложной формы, так как значительно усложняется про­ грамма расчета. В этом отношении наиболее целесообразен метод элементарных балансов [5], основанный на разбивке всего рассматри­ ваемого тела на ряд элементарных объемов, в пределах которых закон изменения температуры принимается линейным. Метод позволяет определить температуру в любых телах при любых граничных усло-

виях. Применение этого метода изложено в работе [80]. Однако при его использовании также необходимо значительное время для расчетов.

Метод регулярного режима [31] широко используют при решении ряда специальных задач, а в отдельных случаях — и при расчете температурных полей, особенно когда исследуется длительный переходный процесс с постоянными граничными условиями.

Следует отметить, что все приближенные методы решения диффе­ ренциальных уравнений требуют значительного времени. С исполь­ зованием ЭВМ это время хотя и сокращается, но для сложных узлов остается большим. К тому же программирование требует длительной, кропотливой работы специалистов. Поэтому решение задачи тепло­ проводности приближенными методами с применением ЭВМ целе­ сообразно в том случае, когда требуется выполнить многочислен­ ные расчеты.

§ 48. Методы аналогии

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

Механическую и акустическую аналогию применяют в ряде спе­ циальных задач. Гидродинамическая аналогия, основанная на тожде­ ственности (при математическом описании) функции тока и потен­ циала скорости, с одной стороны, и функции теплового потока и тем­ пературы, — с другой, позволяет решить задачи теплопроводности с помощью гидродинамического поля. Этот метод разработан до­ вольно глубоко, и с его помощью можно решать двух- и трехмерные нестационарные задачи и нелинейные уравнения параболического типа. Однако гидродинамические модели требуют больших экономи­ ческих затрат, очень громоздки, сложны в эксплуатации и создают определенные затруднения при решении задач с переменными теплофизическими характеристиками.

Вэтом отношении электроаналогия имеет большие преимущества

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

просты по своему устройству, могут быть выполнены с большой степенью точности и сравнительно дешевы в изготовлении.

Все существующие электрические модели можно разделить на две большие группы: 1) модели из сплошных сред и 2) модели из электри­ ческих сеток.

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

в ваннах

применяют слабые растворы электролитов, водопроводную

и

дистиллированную

воду.

 

 

 

С целью

увеличения

 

точно­

 

 

 

сти

измерений

и снижения

 

 

 

погрешности

на

электроды

 

 

 

подают ток большой частоты.

 

 

 

Модели

ванн

бывают объем­

 

 

 

ные и плоские в зависимости

 

 

 

от

решаемой

задачи.

Сама

 

 

 

ванна

обычно

коробчатого

 

 

типа

и изготовлена из изоля­

 

 

 

ционных материалов. Внутри

 

 

ванны

размещается

 

испы­

 

 

тываемая

модель, внутри ко­

 

 

 

торой

находится

электролит.

Рис.

106. Схема электролитической

ванны

При

необходимости

модели­

для определения температурного поля

охлаж­

рования различной

тепло­

даемой

лопатки.

 

проводности

можно заливать

 

 

 

в отдельные области электролит различной концентрации. К наруж­ ным и внутренним стенкам модели через специальные устройства подводятся граничные напряжения от индивидуальных делителей. Измерение потенциалов внутри модели производится специальным зондом.

Схема электролитической ванны для определения температур­

ного поля охлаждаемой

лопатки представлено на рис. 106.

Сопротивления Rx и

R2 определяют

по

уравнениям:

 

X

Р

_

X

Р .

 

« 1

FЭК

~

« 1

hi,

'

 

А .

Р

_

х

Р

 

 

а .

Рэп

 

И 2

 

 

где л—теплопроводность материала лопатки; коэффициент теплоотдачи от газа к лопатке;

коэффициент теплоотдачи от стенки канала лопатки

кохлаждающему агенту;

Рудельное сопротивление проводника;

Р

эк

эквивалентная поверхность электродов;

1

Л — толщина проводника;

 

 

шаг электрода.