Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf

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

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

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

Добавлен: 27.06.2024

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

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

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

сматриваемую задачу излагаемыми

ниже методами;

почти все

эти допущения были введены Клигелем в работе {65].

 

1. Среда является двухскоростной

и двухтемпературной, т. е.

в каждой точке потока имеются две скорости (скорость газа и

скорость частиц) и две температуры

(температура

газа и тем­

пература частиц). При этом совокупность частиц считается неп­ рерывно распределенной по всему объему с условной плотно­ стью «газа» частиц равной произведению численной концент­ рации частиц в единице объема на массу одной частицы.

2.Давление создается только газом; влиянием частиц пренебрегается.

3.Течение — стационарное.

4.Массовый расход газа и массовый расход частиц вдоль по течению постоянны.

5.В любом поперечном сечении все параметры постоянны.

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

7.Система теплоизолирована; обмен теплом имеет место лишь между частицами и газом и осуществляется только путем конвекции.

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

9.Объемом, занимаемым частицами, можно пренебречь.

10.Вследствие высокой теплопроводности материала частип их температура по всему объему частиц постоянна.

11.Гравитационными и электрическими силами можно пре­ небречь.

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

Уравнение движения частицы можно представить в виде

лг3о, d w s С д Л Г2 Q(W— WS)2 dt

где г — радиус частицы; рв — плотность вещества жидкой или твердой частицы; w и ws — соответственно скорость газа и час­

тицы; Q— плотность газа;

cD— коэффициент

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

сферы.

 

 

При стоксовском режиме обтекания [84]

 

где число Рейнольдса потока газа, движущегося

относительно

частиц,

6 I w WS12r

 

рс

 

 

•П

 

здесь т] — коэффициент вязкости газа.

 

14


Число Маха при движении газа относительно частиц опреде­ ляется аналогично по формуле

I WWs I

М

а

I

В общем случае, когда энтальпия газа і і(р, Т), скорость звука

дд

дд / ді \ —1/_1_

ді

'dp'+ ~ d F \ d T j

{ Q ~

др

- 1/2

Если же энтальпия газа і = с ѵТ,

то скорость звука а = у/~ уДТ, где х =

= С р / ( с р R ) показатель адиабаты

газа.

Вводя функцию сопротивления

f D(M ,R e)= -^-p- , формулу

(1.1) можно преобразовать к виду

 

 

dt

= Т

(М’ К6) ~Т~ (™ - ws)-

(1 • 2)

2

 

r2QB

 

Отсюда

 

 

 

 

 

ws —— = Vl(w —ws),

(1.3)

dx

где

4

D r2ga

Конвективный обмен теплом между частицами и газом опи­ сывается соотношением

 

“7" лг36всв ~ ~ =

“ 4лг2а (Ts—T),

(1.4)

где

св — теплоемкость вещества жидкой или твердой

час­

 

тицы;

 

 

 

 

Т и Ts— соответственно температура газа и частицы;

 

 

а — коэффициент теплоотдачи от частиц к газу.

 

 

Введем числа Нуссельта и Прандтля:

 

 

Nu

2га

Рг

 

 

 

 

к

где М — коэффициент теплопроводности газа. Тогда из соотношения (1.4) получим

* > * ^ = Ъ ( Т - Т 3),

(1.5)

dx

 

 

 

где

 

 

 

3_ _Nu_

V I

NuyiCp

 

2 Рг

е всвг2

ЗРг f DcB

 

15


Отметим, что размерность параметров <рі и <рг выражена в с-1. При постоянных значениях q>4, фг, w и Т из соотношений

(1.3) и (1.5) следует

Д —Д0

 

где под А понимается скоростное (wws)

или температурное

( TsТ) отставание, под А —соответственно

ср^1 или «р ^1, а До

соответствует А при ^ = 0. Отсюда ясен физический смысл пара­

метров ф!

или фг:

величины

1/фі и 12 численно равны [време­

ни, в течение которого отставания А

уменьшаются

в е раз. Не­

которые авторы (например, [126]) называют

параметры

9Г1

и с?2~1 временами релаксации.

 

 

 

 

 

Уравнения сохранения расхода газа и частиц имеют вид

 

 

 

qwF = m = const,

 

 

(1-6)

 

 

QswsF = msWm = const,

 

 

(1.7)

где F — площадь проходного сечения;

 

 

 

W — отношение расхода частиц к расходу газа.

 

 

В соответствии с допущением 4 W = const.

 

 

в мо­

Применим теорему импульсов

к среде, находящейся

мент t в элементарном объеме Fwdt, заключенном

между сече­

ниями F и F + dF. Изменение количества движения среды,

нахо­

дящейся в данном объеме,

отнесенное к периоду

времени dt,

равно d(QW+ QswsFw)-,

кроме того, необходимо учесть изменение

импульса

частиц,

проходящих

за

время dt

через границы:

•—ü?[tos/rQs(ffif—щ5)].

Приравнивая

общее изменение

количества

движения

изменению

внешних сил

давления: —d(pF) + pdF, и

пользуясь

уравнениями (1.6) и (1.7), придем

к

следующему

соотношению

QW dw JrQswsdws~\-dp = 0.

 

 

(1.8)

 

 

 

 

Аналогично может быть получено

и уравнение

сохранения

энергии

 

 

 

 

 

 

 

 

 

г + Qi

F (w

+d(pwF),

где е и es— соответственно внутренняя энергия газа и частиц. Произведя тождественные преобразования и интегрирование,

получим

^ + и ф + ^ ) + у = соп81.

(1-9)

16



Полагая, что для газа выполняется уравнение Клапейрона

P = qRT

(1.10)

и соотношение е — сѵТ, R = cpсѵ, а для частицы es — cBTs, полу­ чим

+ + W (св7\ + = Е0= const. (1.11)

Дифференциальные уравнения (1.3), (1.5) и (1.8), а также конечные соотношения (1.6), (1.7), (1.10) и (1.11) определяют семь неизвестных (р, q, qs, w, ws, T и Ts) в любом сечении сопла заданной формы.

Для удобства вычислений преобразуем полученную систему. Исключим Q, Qs, р и dT из уравнений — расхода, Клапейрона и энергии и подставим в уравнение импульсов (1.8). После пре­ образований получим

d ln w __

1 d ln F I

W

M2

(*— 1)

Ws_

d w , I

dx

 

 

w

—Г- +

M2 — 1 dx ^ М2 — 1

 

wdx

 

+ -

 

dTs

 

 

( 1. 12)

 

 

Tdx ) '

 

 

где показатель адиабаты

 

 

 

 

 

 

w2

 

 

 

 

a M2

 

 

 

cV

%RT

 

 

 

 

 

 

 

В правой части дифференциального уравнения (1.12) в соот­ ветствии с уравнениями (1.3) и (1.5) отсутствуют неизвестные производные. Его численное интегрирование можно производить, например, методом Рунге—Кутта.

Отрицательный член, стоящий в правой части этого уравне­ ния в фигурных скобках, обусловливает сдвиг критического се­ чения, т. е. сечения, где М = 1, в расширяющуюся часть сопла, что характерно для неизэнтропических течений. Как следует из уравнения (1.12), сдвиг увеличивается с ростом содержания час­ тиц в газе, а также с ростом производных dws/dx и dTs/dx.

Соотношение (1.12) можно представить и в несколько иной форме. Выразим w через М и Т и исключим dT с помощью уравнения энергии. Окончательно получим

 

 

( М 2 — l ) 'r f M 2

_ _

 

 

 

2М2 Л +

 

МА

 

 

 

dF

■1 1 + хМ2

W

,

I

1 \

Wwdw.s

(1.13)

F

X 2 + (х — 1) М2 RT

(■cBdTs-{-wsdws)

RT

 

 

 

 

Уравнения (1.3), (1.5), (1.11) и (1.13) позволяют определить неизвестные М (х), Т(х), ws(x) и Ts(x), так как скорость w(x)