Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Турбулентное смешение газовых струй

..pdf
Скачиваний:
10
Добавлен:
24.10.2023
Размер:
10.28 Mб
Скачать

160 Обобщение данных экспериментального исследования [гл. III

закономерности затухания относительной избыточной ско­ рости и температуры и искать решение задачи в виде за­

висимости координаты хс от основных параметров течения. При этом необходимо определить только зависимость

с1п(ж).

4. Для анализа закономерностей изменения концент­ рации вдоль струи при наличии спутного потока (т Ф 0 )

обратимся к условию сохранения потока массы примеси и используем интегральное соотношение (3.16). Обозначив

Uffl

U2

Аи ,

и и 2

= /01).

С

 

ф л

U1 — U2

тп’ Uyyi — U2

с т

 

 

 

(

).

получим

 

 

 

 

 

 

 

 

 

 

 

х

х

 

 

 

 

ст(&°)i+1 (1 — т) Аи°т^ fyrfdf] +

т ^ фгрйц

=

-£-'<?, (3-44)

 

 

 

о

о

 

 

гср

 

 

 

 

 

 

 

 

 

Г

 

1

 

 

 

 

 

-4 - Cm(b°)i+1

(1

т) Аит^ ф2/n ?dn +

 

 

 

 

 

 

 

О

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

/нур2г|Ч2г)

=

-

(ЬУ~ D el,

где (D — D/RuД.

Для больших удалений можно в квадратных скобках

пренебречь слагаемым, содержащим Аи°т. Тогда из двух указанных выражений получим следующие соотношения:

dcm ------c l o

dx° — для осесимметричного течения,

dcm------ clD

dx° — для плоского течения.

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

§ 4] Параметрическое описание течения п основном участке ПИ

соотношения (3.44) при Аи„ -> 0:

стЪ ^ pi Q

1 с р

Отсюда получаем

d —-----О ^ - d x ° — для осесимметричного течения.

с т

 

fJl

1

--- /

Р,

d ~ ~ D '

т dx° — для плоского течения.

т

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

 

(3.45)

для плоского течения

 

С-т

(3.4G)

 

т п ср V D°xa

где D° — Dlu^R = D/m.

В соотношениях (3.43) и (3.44), а также в последующих выкладках используется понятие средней плотности рср. Нужно отметить, что значение средней плотности в разных соотношениях ив левой и правой частях этих соотношений не одинаковы. Тем не менее для получения качественных зависимостей,характеризующих в основном приближенную пропорциональность одних величин другим, представ­ ляется возможным произвести соответствующие дей­ ствия с параметром рср как с единой алгебраической ве­ личиной.

Соотношения (3.45) и (3.46) позволяют описать законо­ мерности развития течения в основном участке струи, если

известна связь характерного

значения коэффициента

турбулентной диффузии с основными определяющими

параметрами течения т, п, 6 °, 6 °

и D°v а также если опре­

делено значение пср. Поскольку в опытах исследовалось осесимметричное течение, в дальнейшем на нем будет сосредоточено основное внимание.

Удовлетворительного описания параметров струи в спутном потоке можно достигнуть, если рассмотреть три

6 Г. Н. Абрамович и др.

1G2 Обобщение данпых экспериментального исследования [гл. III

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

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

Di ~ (\ I \ /(>)*-• ,

(3.47)

в плоском течении — комплекс

Д 1 ~ |/„ |/р н .

(3.48)

Очевидно, что для больших удалений характерными значениями плотности и скорости являются их значения в спутном потоке, т. е. перемешивание, связанное с сило­ вым воздействием струи на поток, описывается коэффици­ ентами турбулентной диффузии, которые определяются соотношениями (3.47) — (3.48).

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

ной диффузии D -г, а также турбулентность, генерируемая в начальной части течения в пограничных слоях, которая характеризуется коэффициентом диффузии Ds и играет доминирующую роль при Г = 0. Можно предположить, что суммарный эффект переноса примеси определяется су­ перпозицией указанных процессов переноса:

D = klD 1 + k2D 2 + k3D 3,

(3.49)

где ки к2 и к3 — некоторые коэффициенты, зависящие от расстояния и основных параметров течения. В этом случае соотношение (3.45) можно представить соответственно для

§ 4] Параметрическое описание течения в основном участке 1G3

осесимметричного и плоского случаев в виде

, (3.50)

" c p l'1 ’ 2I 7° Г 2+ k 2 nc\nD2 m + k3ncv3mD3

а соотношение (3.4(>) — в виде

-Т- .(3.51)

Km'*-117П\+к2пп1Р2тЮ1+ kmnl PamtV*)

Коэффициент к'а зависит от параметра т, поскольку третье слагаемое рассматриваемой суммы (3.49) описы­ вает характеристики турбулентного переноса, обусловлен­ ные исходной неравномерностью профилей скорости в по­ граничных слоях. При m ~ 0 u m 1 исходная неравно­ мерность с достаточной точностью учитывается величиной избыточного импульса поэтому

при т ~ 0

и т J?;> т%к'а = 0,

(3.52)

где т* соответствует

I° (т%) — 0 .

 

При фиксированных значениях п, 6 Ь 62

и Daмаксималь­

ное значение коэффициент к3 принимает

при т == т*.

Из условий (3.52) можно определить для течения в за­

топленной струе и течения за источником, не возму­ щающим поток, коэффициенты к[ и кг’.

В случае затопленной струи = 0) величина /° ж 1

(с точностью до величин и б р), поэтому для осесиммет­ ричной струи имеем

для плоской струи

с —■V ср 1

Сопоставляя полученные соотношения с зависимостями

6*

104 Обобщение данных экспериментального исследования [гл. III

(3.40) и

(3.41), определим значения коэффициентов к[ и

параметра »cpi:

осесимметричной струи),

= 11 (для

к{ =

9,5 (для

плоской

струи),

ср1 = п.

 

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

ьк

\сиуЧу = ^ сиуЧу.

оо

Всоответствии с условиями задачи и = const. При этом

1

2г Ф°)г+] Ф(Л) Л^Л

о

Профили концентрации на больших удалениях извест­ ны, поэтому имеем

5 ф (л)

^ 0,45, ^ ер (л) л йц ~ 0,129.

о

о

Ширину профиля концентрации можно связать с продоль­ ной координатой и коэффициентом турбулентной диффу­ зии, используя определение коэффициента диффузии:

 

 

U

dY1

(3.53)

 

 

D = li in

dt

 

где Y 2 — дисперсия примеси.

 

 

При

больших

удалениях и обычном условии переноса

х = ut

имеем

У* _ 2Рх

 

 

 

 

 

 

Для струйных

профилей величина дисперсии связана

с шириной профиля соотношением

 

У2 ж 0,14 Ъ\

§ 4] Параметрическое описание течения в основном участке 165

В соответствии с этим получаем соответственно для осе­

симметричного

и плоского

течений

 

 

1

И

______ 1

(3.54)

Ст ~

3,68.0°ж°

т ~~ У 2,891Гх°

 

 

Остается установить связь между значением макси­ мальной концентрации ст и параметрами т и п в равномер­ ном потоке.

На больших расстояниях от источника, где ширина профиля концентрации в соответствии с соотношением (3.53) определяется значением коэффициента турбулентной диффузии и не зависит от размеров источника, характер­ ное значение концентрации в данном сечении зависит от соотношения потоков массы, втекающих в пределы струи:

_

G\ __

Gi___p m i____ 1_

m

C H -G 2~

G‘i p2u2 m u'

Следовательно, при т ф 1 и п ф 1

в соотношениях (3.54)

может быть внесена поправка:

 

 

1

С-т

i

(3.55)

°т ~ 3,68mnD'’x°

пгп У 2,89JL>*x“

 

 

Из сопоставления соотношений (3.50), (3.51) и (3.55) имеем

к2~ к2п~ 34, пср2 = п.

6 . Рассмотрим случай нулевого избыточного импульса

Г — 0. Как уже указывалось, этот случай соответствует условию т ^ 1 , т. е. такому превышению скорости в ис­

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

они весьма существенны и именно они определяют сме­ щение. Это значит, что интенсивность турбулентного

166 Обобщение данных экспериментального исследования [гл. III

перемешивания не уменьшается до нуля при —►0 , даже

если турбулентность во внешнем потоке отсутствует. Та­ кое неограниченное уменьшение интенсивности турбулент­ ного перемешивания могло бы реализоваться лишь при условии абсолютно равномерного распределения скорости на срезе соплового устройства, т. е. в том случае, когда

= 0

достигается при

= 1 .

ско­

Течение

при неравномерном

начальном профиле

рости и =

0 рассматривалось

применительно к задаче

о следе

за

движущимися

телами. Теоретически

эта

задача исследовалась в работах [15, 65]; в работе [6 6 ]

приведены результаты экспериментального изучения те­ чения за телом специальной формы при /° = 0. Теоре­ тический анализ относится либо к ламинарным течениям, либо к таким турбулентным течениям, в которых характе­ ристики турбулентного переноса определяются локаль­ ными градиентами скорости. Поэтому полученные в ука­ занных работах результаты не могут быть использованы для анализа рассматриваемого здесь случая течения, по­ скольку в соответствии с данными опытов было сделано предположение о постоянстве характерных значений па­ раметров переноса вдоль течения.

В случае нулевого избыточного импульса при турбу­ лентном течении определяющая роль исходных возмуще­ ний еще более очевидна, чем в рассмотренных уже слу­ чаях. Теоретический анализ показывает, что в таком те­ чении из-за быстрого уменьшения характерного дефекта скорости и слабого роста характерной ширины ее про­ филя [15, 65]происходит интенсивное (более интенсивное, чем в спутной осесимметричной струе) уменьшение харак­ терной турбулентной вязкости (диффузии), связанной с ло­ кальными градиентами. Поэтому естественно ожидать, что при нулевом значении избыточного импульса, так же как и в других случаях течения, в спутном потоке смеше­ ние будет определяться турбулентностью, сносимой пото­ ком (при постоянном значении коэффициента турбулент­ ного переноса).

Это подтверждается результатами измерения характе­ ристик турбулентности в следе за телом специальной фор­ мы при = 0 , приведенными в работе [6 6 ], согласно ко­

торым максимальное значение турбулентной вязкости практически постоянно вдоль всего течения.

I 4] Параметрическое описание течения в основном участке 167

Можно считать, что в реальных условиях турбулент­ ного течения при = 0 характеристики турбулентного

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

Рассмотрим осесимметричное течение. В этом случае для выражения величины относительного избыточного импульса через основные параметры течения можно ис­ пользовать при малой толщине пограничных слоев следую­ щее приближенное соотношение:

оо

Г2 ^ рп щ) у dy • (pxu\R2)~x

о

~ 1 — то — 2бх — 2б1р (1 — то) — 2то2п62. (3.56)

Пренебрегая величиной 2б1р по сравнению с едини­ цей, получим

1 — то* ~ 2 (6 t + то^г62).

(3.57)

Соотношение (3.57) показывает, что влияние наруж­ ного пограничного слоя на смещение режима течения с ну­ левым избыточным импульсом от значения то, равного единице, в значительной мере зависит от относительной плотности наружного потока п. Чем больше п, тем сильнее проявляется влияние наружного пограничного слоя.

Соотношения (3.56) и (3.57) могут быть использованы для определения характеристик турбулентного переноса при то то* (7° = 0). Если значение то* близко к единице, то характеристики турбулентности будут близки к харак­ теристикам турбулентности в следе за пластиной (в кото­ рую может быть предположительно развернут контур

168 Обобщение данных экспериментального исследования [гл. II

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

Коэффициент турбулентной диффузии определяется согласно соотношению (3.48) силовым воздействием такой пластины на поток в целом:

Тогда в соответствии с обозначениями, принятыми в соот­ ношениях (3.45) и (3.57), при т 1 имеем

6j + тЛ(6°

ср

Если значение параметра т заметно отличается от единицы, то характеристики турбулентности будут опре­ деляться сдвиговым течением. Используя соотношение (3.20) для зоны смешения, получим

Dz

1+ п (1 — т /

«1

2 1 -р тп

Протяженность области, где в основном порождается турбулентность при 7° = 0, пропорциональна длине на­ чального участка. Ниже по течению порождение турбу­ лентности в этом случае относительно невелико и величина коэффициента турбулентной диффузии соответствует его значению в конце начального участка струй. Протяжен­ ность последнего при значении иг*, заметно отличающем­ ся от единицт.1, может быть найдена из соотношения (3.34)

П ( 1 + тп)

X

Из двух последних соотношений имеем

о1 т

(3.59)

Явления, связанные с порождением турбулентности в сле­ де за пограничными слоями, а также в зоне смешения, по-

§ 4] Параметрическое описание течения в основном участке 169

видимому, достаточно полно отражают процессы, проис­ ходящие в потоке при = 0. Оба итоговых соотношения (3.58) — (3.59) для характерных значений коэффициентов турбулентной диффузии при = 0 показывают, что ком­ плекс ncpD°, который согласно (3.45) определяет вели­ чину ст , в обоих случаях слабо зависит от отношения плотностей п и практически зависит только от разности (1 — та*). При этом суммарный эффект, вызываемый по­ рождением турбулентности в начальном участке струи, дает значение коэффициента турбулентной диффузии, про­ порциональное разности (1 — та*):

ncpD°a — 1 — та, —

+ талиба.

(3.60)

Это соответствует значениям

пср = пср. с Для

следа и

геср = (Vn/( 1 + /

п)?‘ для

слоя смешения.

 

Аналогичным

образом

можно

проанализировать слу­

чай Г — 0 для плоского течения, для которого связь от­ носительного избыточного импульса с определяющими параметрами течения дается уже точным соотношением, аналогичным (3.56):

Г = $57Г (I T “

m) dy0 +

( - £ - l ) dlf =

 

 

= (1 — та) (1 — 6iP) — 6i— таг/гб2.

(3.61)

Пренебрегая

 

о

 

величиной б1р по сравнению с единицей,

нолучаем

 

 

 

 

1 — та* ~

+ пСпЬ.,.

(3.61')

Так же как и в осесимметричном течении (3.57), влия­ ние наружного пограничного слоя на смещение режима течения с нулевым избыточным импульсом от значения та = 1 к значениям то < 1 зависит от относительной плот­ ности наружного потока п. Однако влияние пограничных слоев на это смещение в плоском случае слабее, чем в осе­ симметричном.

Если значение та* близко к единице, то по аналогии с (3.58) получаем

DK

б1 + т°. ”62

(3.62)

с р . с

с р . с

 

Соседние файлы в папке книги из ГПНТБ