книги из ГПНТБ / Турбулентное смешение газовых струй
..pdf160 Обобщение данных экспериментального исследования [гл. 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) позволяют описать законо мерности развития течения в основном участке струи, если
известна связь характерного |
значения коэффициента |
турбулентной диффузии D° с основными определяющими |
|
параметрами течения т, п, 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
перемешивания не уменьшается до нуля при 1° —►0 , даже
если турбулентность во внешнем потоке отсутствует. Та кое неограниченное уменьшение интенсивности турбулент ного перемешивания могло бы реализоваться лишь при условии абсолютно равномерного распределения скорости на срезе соплового устройства, т. е. в том случае, когда
1° = 0 |
достигается при |
= 1 . |
ско |
||
Течение |
при неравномерном |
начальном профиле |
|||
рости и 1° = |
0 рассматривалось |
применительно к задаче |
|||
о следе |
за |
движущимися |
телами. Теоретически |
эта |
задача исследовалась в работах [15, 65]; в работе [6 6 ]
приведены результаты экспериментального изучения те чения за телом специальной формы при /° = 0. Теоре тический анализ относится либо к ламинарным течениям, либо к таким турбулентным течениям, в которых характе ристики турбулентного переноса определяются локаль ными градиентами скорости. Поэтому полученные в ука занных работах результаты не могут быть использованы для анализа рассматриваемого здесь случая течения, по скольку в соответствии с данными опытов было сделано предположение о постоянстве характерных значений па раметров переноса вдоль течения.
В случае нулевого избыточного импульса при турбу лентном течении определяющая роль исходных возмуще ний еще более очевидна, чем в рассмотренных уже слу чаях. Теоретический анализ показывает, что в таком те чении из-за быстрого уменьшения характерного дефекта скорости и слабого роста характерной ширины ее про филя [15, 65]происходит интенсивное (более интенсивное, чем в спутной осесимметричной струе) уменьшение харак терной турбулентной вязкости (диффузии), связанной с ло кальными градиентами. Поэтому естественно ожидать, что при нулевом значении избыточного импульса, так же как и в других случаях течения, в спутном потоке смеше ние будет определяться турбулентностью, сносимой пото ком (при постоянном значении коэффициента турбулент ного переноса).
Это подтверждается результатами измерения характе ристик турбулентности в следе за телом специальной фор мы при 1° = 0 , приведенными в работе [6 6 ], согласно ко
торым максимальное значение турбулентной вязкости практически постоянно вдоль всего течения.
I 4] Параметрическое описание течения в основном участке 167
Можно считать, что в реальных условиях турбулент ного течения при 1° = 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
видимому, достаточно полно отражают процессы, проис ходящие в потоке при 1° = 0. Оба итоговых соотношения (3.58) — (3.59) для характерных значений коэффициентов турбулентной диффузии при 1° = 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) |
с р . с |
с р . с |
|