Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.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 — 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