Файл: Жаров Г.Г. Судовые высокотемпературные газотурбинные установки.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 230

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

Временные краевые условия сводятся к заданию скалярной функ­ ции, которая определяет распределение температуры в теле в неко­ торый момент времени

t = f (х; у; г).

Пространственные (граничные) краевые условия могут быть за­ даны четырьмя способами:

— температурой поверхности тела в любой момент времени (гра­ ничные условия первого рода);

•— тепловыми потоками на поверхности тела в любой момент времени (граничные условия второго рода);

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

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

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

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

§45. Решение дифференциального уравнения теплопроводности классическими методами

Выведенное в § 44 уравнение теплопровод­ ности (137) является однородным линейным дифференциальным уравнением второго порядка с частными производными. Такое диф­ ференциальное уравнение имеет бесчисленное множество частных решений. Классический метод решения состоит в том, что находятся



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

 

л

 

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].

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

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

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

от

изображения

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

к

изображению

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

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