Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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 W— Ws 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 |
|
где под А понимается скоростное (w— ws) |
или температурное |
( Ts— Т) отставание, под А —соответственно |
ср^1 или «р ^1, а До |
соответствует А при ^ = 0. Отсюда ясен физический смысл пара
метров ф! |
или фг: |
величины |
1/фі и 1/ф2 численно равны [време |
||||||
ни, в течение которого отставания А |
уменьшаются |
в е раз. Не |
|||||||
которые авторы (например, [126]) называют |
параметры |
9Г1 |
|||||||
и с?2~1 временами релаксации. |
|
|
|
|
|
||||
Уравнения сохранения расхода газа и частиц имеют вид |
|||||||||
|
|
|
qwF = m = const, |
|
|
(1-6) |
|||
|
|
QswsF = ms— Wm = 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)