Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 171
Скачиваний: 4
Ш ^ |
Л'+Я « ■ ^•= а |
(4.2) |
|
||
х |
5 |
|
Здесь т — некоторый фиксированный объем, ограниченный |
по |
|
верхностью S, а п — единичная внешняя нормаль к этому |
объ |
|
ему. Напомним также, |
что входящие в эти формулы плотности |
Q и Qs, представляют собой отношения соответственно масс газа и частиц к объему, в которой они заключены, и не являются по этому собственно плотностями газа и вещества частиц.
Уравнение движения двухфазной среды можно |
записать в |
|
виде |
|
|
X |
X |
|
- f f |
[pn + QaV s [(Vs- V ) n ] } d S , |
(4.3) |
•Js*J |
|
|
где Ps и P—соответственно векторы внешних массовых сил, дей ствующих на частицы и на газ; знак минус перед поверхностным
интегралом означает, что нормаль п — внешняя; второе слагае мое в этом интеграле отражает приобретение (или потерю) коли чества движения вследствие захвата (или потери) объемом т новых частиц. Следует отметить, что в отличие от уравнений (4.1) и (4.2) уравнение (4.3) записано для объема т, связанного с газом и, следовательно, подвижного, поэтому внесение произ водной по времени под знак тройного интеграла, стоящего слева, недопустимо. В то же время При переходе от подвижного к фик сированному контуру вид слагаемых, стоящих справа, остается тем же.
Примёним к |
левой части формулы (4. 3) |
известное соотно |
||
шение (см., например, работу [84], с. 96) |
|
|
||
т |
т0 |
5 |
|
, 4 4 ) |
|
|
|||
где А ■ произвольный вектор; |
|
|
|
|
г — объем, скорость перемещения |
границы которого V= |
|||
= V ( S ) , |
|
|
у X в дальнейшем |
|
to— фиксированный объем (индекс «О» |
||||
опускается). |
|
|
|
интег- |
Приводя подобные члены под знаком поверхностного |
||||
рала, получим вместо формулы (4.3) |
|
|
|
|
7 Г \ \ \ |
= ^ |
( е ^ + о Л ) * - |
|
|
|
J j [/ш + еН(Ня) + е*ИД1/^)]с/5. |
(4. 5) |
127
Интегральное уравнение (4.5) удобно тем, что оно записано применительно к фиксированному объему т и может быть ис пользовано для непосредственных расчетов. Кроме того, оно име ет ясный физический смысл; изменение по времени количества движения в фиксированном объеме т определяется воздействием внешних массовых сил, работой сил давления на границе объема и потоком импульса газа и частиц в этот объем.
На основании неизменности объема т полная производная d/dt, стоящая перед тройным интегралом, может быть внесена под его знак и заменена частной производной.
Уравнение энергии в интегральной форме выводят аналогич но. При интегрировании по подвижному объему т, связанному с газом, получаем
|
|
|
|
|
|
(4. 6) |
где |
е и es — соответственно |
внутренняя энергия |
газа |
|||
|
частиц; |
|
|
|
получать газ |
|
|
gQdx и QsQsdx — внешнее тепло, которое могут |
|||||
|
и частицы, находящиеся в объеме dr; |
|
||||
|
отношения q / qp и qs/qb соответствуют доле |
поверхности dS, |
||||
занятой газом или частицами (здесь |
рг — собственно плотность |
|||||
газа, а рв— собственно плотность вещества частиц). |
знаком по |
|||||
|
Таким образом, первое слагаемое, |
стоящее |
под |
|||
верхностного интеграла, характеризует работу |
сил давления на |
|||||
границе объема dt, а второе — энергию, вносимую |
или уноси |
|||||
мую частицами, попадающими в |
объем dt или уходящими |
из |
||||
него. |
|
|
|
|
|
|
|
Применим к первому слагаемому |
(4.6) формулу (4.4). Полу |
||||
чим для фиксированного объема t |
интегральное уравнение |
|
||||
|
d_ |
|
|
dt — |
|
|
|
dt |
|
|
|
|
-Ш [Q(-Q+^)+Q*(Qe+^y,)]rft+
t
f lf S = 0 .
(4. 7)
128
Уравнение (4.7), как и (4.5), также легко интерпретируется физически: изменение энергии, заключенной в фиксированном объеме т, равно подведенной энергии внешних источников, ра боте массовых сил, сил давления (на границе), а также энергии газа и частиц, проносимой через границу. Знак полной производ ной, стоящей перед первым интегралом уравнения (4.7), может быть внесен под интеграл и заменен частной производной d/dt.
Уравнения движения и энергии частиц можно представить в виде
dVs
dt f ~\-Ps — Qв 4P’
^ = g + Q s , dt
где полные производные взяты вдоль траекторий частиц, f и q — соответственно сила, обусловленная вязкостью, и тепловой по ток, отнесенные к единице массы частиц и вызванные взаимодей ствием частиц с газом. Последнее слагаемое в формуле для ус корения частиц соответствует воздействию градиента давления на поверхность, занятую частицами.
В эквивалентной форме эти уравнения получат вид
( ? , Ѵ ) ? 4 + - ^ |
+ — Ѵ / > - / - Я , = 0 , |
( 4 . 8 ) |
dt |
QB |
|
( V s V ) e , + - ^ - q - Q , = 0 . |
(4.9) |
dt
Величины f и q могут быть представлены в виде f = i y - V , ) K
q = (T— Ts)<?2,
где <рі и <р2— функции от р, Т, Ts и (V— Vs).
Плотности собственно газа и частиц, а также их внутренние энергии связаны с давлением и температурой формулами
аг = Ѳ г ( л Л ; е=е{р,Т)\ |
) |
|
QB= const; |
es^=es{Ts). |
j |
Вместо внутренней энергии газа удобнее |
пользоваться его эн |
|
тальпией |
|
|
i = e{p,T) + p / Q r = i(p,T).
В подобластях непрерывности параметров соотношения (4.1), (4.2), (4.5) и (4.7), записанные в интегральной форме, могут быть заменены дифференциальными уравнениями.
5 |
3739 |
129 |
Воспользуемся формулой Гаусса—Остроградского
Я™ М Я « '
S т
и произвольностью объема т; из выражений (4.1) и (4.2) в под областях непрерывности параметров получим
^ - + ѵ ( е іо = о |
и % - + v (Civ g = o . |
(4 .ii) |
dt |
dt |
|
Эти соотношения часто записывают в виде
■^ + div (рІ/) —0 и - ^ - + div(e,K,) = 0, |
(4.12) |
причем следует отметить, что для осесимметричного течения
div V == -д и |
, |
д ѵ |
, V |
дх |
4- |
dtf |
У |
а для плоского — последнее слагаемое отсутствует. Аналогич ным образом из уравнения (4.5) следует
~ (QV+ е Л ) - Q P - QSPS+ Ѵр+ VQl/2+ V Q ^ = о .
Заменяя здесь Ps по формуле |
(4.8), используя |
уравнения |
|
(4.11), тождество |
|
|
|
— |
+ |
|
(4.13) |
Q |
Qr |
QQb |
|
и выполняя преобразования, получим |
|
||
^ ѵ ) ^ + 4 г + — ѵ ^ + — 7 - ? = 0 - |
(4 Л 4 ) |
||
dt |
рг |
6 |
|
Подобным образом в подобластях непрерывности параметров можно получить и уравнение энергии в дифференциальной фор ме.
Из уравнения (4.6) получаем
_д_ |
<4 |
|
+*.) - ö C Q + P i/) - e ,( Q ,+ /> ^ ) + |
|||
dt |
|
|||||
+ ѵ A l f ѵ + іг-ѵ* |
|
Vi |
= 0 . |
|||
+ * , { V - V ) |
||||||
Qr |
Qb |
|
|
|
|
|
Используя равенства |
(4.8) — (4.11), (4.13), а также |
следую |
||||
щее из равенств |
(4. 11) |
и |
(4. 13) |
тождество |
|
|
|
V Q ^ |
I |
6 |
d Q r |
I V Q s ^ _ _ q |
|
|
6г |
|
|
dt |
бв |
|
130
и выполняя преобразования, получим уравнение энергии в срав нительно простом виде
---- Цѵ'Ѵ/7 + ^ -) + ѵѴ = 0, |
(4.15) |
|
dt Qr V |
dt 1 |
|
где ^ = - ^ [ ( К , - У ) / + ^]-(3 . e
§ 4.2. Классификация сильных разрывов
Уравнения сохранения расхода, количества движения и энергии, записан ные в интегральной форме (4.1) — (4.4), позволяют получить соотношения на сильных разрывах.
Рассмотрим только плоские разрывы, так как в окрестности любой неосо бой точки произвольной поверхности разрыва достаточно малый элемент ее всегда можно считать плоским.
Выберем систему координат, движущуюся с постоянной скоростью, совпа дающей в данный момент со скоростью поверхности разрыва и расположен ной так, чтобы одна из ее осей была нормальна к плоскости разрыва1. Выде лим элементарный цилиндр с образующей, параллельной этой оси, и основа ниями, лежащими по обе стороны от плоскости разрыва на малых расстояни ях от нее. При переходе в объемном интеграле к пределу, при сближении по верхностей, находящихся по обе стороны от разрыва, учитывая непрерывность
производных от всех параметров по времени |
и по осям, лежащим |
в плоско- |
||
«■> |
—> |
*+ |
q, Q и Qs, получим |
|
сти разрыва, а также непрерывность /, |
Р, |
Р,, |
|
|
[QK„]=o И |
[ е ^ да] = о , |
(4.16) |
где индекс «п» означает, что составляющая рассматриваемого вектора (в дан ном случае — скорости) нормальна к плоскости разрыва; квадратные скобки заменяют разность стоящих в них параметров за и перед разрывом.
Условиям (4.16) эквивалентны соотношения
ОУп — 3 и |
s n ~ Js> |
(4.17) |
причем j и j s непрерывны при переходе через разрыв. |
и (4.4) сво |
|
Интегрирование по объему цилиндра |
в соотношениях (4.3) |
дится к интегрированию по направлению, нормальному к плоскости разрыва. Из уравнения (4.5) следует
и ] + |
7 - М |
= |
° : |
(4.18) |
|
Ув |
|
|
|
U [ ^ ] |
= о. |
|
|
|
Далее, из выражения (4.3), получаем |
|
|
|
|
У [ V + у V ,] + [р] « + [ б ,? , (V ^ - ѵ л) ] = 0 , |
|
|||
откуда, использовав соотношения (4.17) и (4.18), получим |
|
|||
3 ІУл] + is ІУот] + [р\ |
= |
0; y[V t] = 0 . |
(4.19) |
|
Из соотношения для теплообмена между частицами и газом |
(4.6) следует |
|||
М > і ] = 0 . |
|
(4.20) |
||
1 Если линии разрыва стационарны, то система координат |
неподвижна. |
|||
5* |
|
|
|
131 |