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

книги / Рудничная аэрология.-1

.pdf
Скачиваний:
6
Добавлен:
20.11.2023
Размер:
16.52 Mб
Скачать

месей (газы, пары, пыль) в воздушных потоках, а также распростра­ нение газов в горных породах. Газовая динамика имеет дело с диф­ фузией газообразных примесей.

Распространение газов в воздушном потоке происходит в резуль­ тате: 1) увлечения их движущимся воздухом, 2) молекулярного перемешивания газов с чистым воздухом, 3) турбулентного их пере­ мешивания. Соответственно различают конвективную, молекуляр­ ную и турбулентную диффузию.

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

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

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

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

n _ L* _ L

TL2 T *

т. e. диффузионный поток, выраженный в объемных единицах, имеет размерность скорости.

Соответственно указанным выше механизмам диффузии разли­ чают три вида диффузионных потоков: конвективный, молекуляр­ ный и турбулентный.

Конвективный перенос осуществляется движущимся воздухом. Следовательно, направление конвективного потока совпадает с на­ правлением осредненной скорости потока в точке.

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

щади этой плоскости будет Q = иЛ = и. Если концентрация диф­ фундирующего вещества в воздухе равна с, то вместе с количеством воздуха Q через площадку в единицу времени будет перенесено

cQ = си объемных единиц диффундирующего вещества. Следова­ тельно, конвективный диффузионный поток будет равен

7К=

си.

(XIV,1)

При трехмерном движении

компоненты

конвективного потока

по осям координат будут:

 

 

х = Ы \

(XIV,2)

h y = CV\

/■ *= *«*,

 

где и, V, w — компоненты скорости по осям Qx, Оу, Oz. Молекулярный и турбулентный диффузионные потоки направ­

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

Молекулярный диффузионный поток определяется градиентом концентрации с и эмпирическим коэффициентом пропорциональ­

ности JDm, называемым к о э ф ф и ц и е н т о м

м о л е к у л я р ­

н о й д и ф ф у з и и :

 

Ум = — grad с.

(XIV, 3)

Компоненты молекулярного потока по осям координат будут:

X ----

дс

9

 

дх

»

 

 

 

 

дс

 

(XIV,4)

1иу= - D * ду

дс

7MZ= - A * ~дГ

Формулы (XIV,3), (XIV,4) носят название п е р в о г о з а к о н а Ф и к а.

Формулы (XIV,4) предполагают, что коэффициент молекулярной диффузии не зависит от направления, т. е. DMх = D My = DMZ. Это следует из кинетической теории газов, согласно которой коэффициент молекулярной диффузии равен произведению свободной длины про­

бега молекул I на среднюю скорость их теплового движения ум, которые не зависят от направления:

D» = \ l v u.

(XIV.5)

Так как /м направлен в сторону уменьшения концентрации, градиент концентрации и ее производные по осям, берущиеся в урав­ нениях (XIV,3) и (XIV,4) в направлении диффузии, имеют отрица­ тельный знак. Поскольку диффузионный поток в направлении диф­ фузии есть величина положительная, правые части в уравнениях (XIV,3), (XIV,4) взяты со знаком минус.

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

U = —DTgrade.

(XIV,6)

Коэффициент пропорциональности DT носит название

к о э ф ­

ф и ц и е н т а т у р б у л е н т н о й

 

д и ф ф у з и и . Он

опреде­

ляется экспериментально.

 

 

 

 

Компоненты /т по осям координат будут:

 

7 т х

х

дс

 

 

дх

9

 

 

 

 

7т у =

D T у

дс

 

(XIV,7)

ду

9

U z = - D TZ

дс

 

 

 

 

dz

 

 

Как видим, в уравнениях (XIV,7) принято,что DTX=£DTy^ D TZ. Объяснение этому будет дано ниже. Знаки минус в уравнениях (XIV,6) и (XIV,7) имеют тот же смысл, что и в (XIV,3) и (XIV,4).

Обмен объемами газо-воздушной смеси при турбулентной диффу­ зии осуществляется под действием пульсационных скоростей. Рас­ сматривая перенос вещества под действием пульсационных скоростей как особый вид конвективной диффузии, можно написать следующие выражения для мгновенных турбулентных диффузионных потоков:

irx = cun; U у —«V.

2==

где с и ип% уп1 и>п— мгновенные значения концентрации и компо­ нент пульсационных скоростей в точке.

Осредняя эти потоки за достаточно большой промежуток времени, получим средние величины турбулентных потоков (черта означает знак осреднения по времени):

h x = ~ u n; jTy = cvn) îrz = cwn-

(XIV,8)

Коэффициент турбулентного обмена для импульса и коэффициент турбулентной диффузии. В § 38 было введено понятие коэффициента турбулентного обмена для импульса е, компоненты которого опре­ деляются выражениями (VI,30) и (VI,27),

TT« = tTM = pe2- J - ,

(XIV,9)

где тТХ2 = тТ2х — касательное

турбулентное

напряжение в пло­

скости хОу,

направленное

вдоль оси Ох;

р — плотность потока.

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

Здесь и далее ось Ох направлена по потоку, Оу — перпендику­ лярно боковым стенкам, Oz — перпендикулярно кровле и почве выработки.

Касательное напряжение в плоскости является результатом

переноса через

нее импульса

 

 

 

т! / = » - * „

(XIV,10)

где P t —- поток

импульса

через единицу площади

плоскости, пер*

пендикулярной

оси Oi.

 

Первый индекс при т указывает, что плоскость, в которой дейст­ вует касательное напряжение, перпендикулярна оси Oi, второй — что данная компонента напряжения направлена вдоль оси О/; Oi,

О)

— общие обозначения координатных осей, в которых i и j могут

принимать

значения х, у, z (например, при i = х Oi = Ох; при

j =

z Оу =

Oz и т. д.).

 

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

пульса в единицу времени через площадку в трехмерном потоке, перпендикулярную оси Oi и движущуюся вдоль нее со средней

местной скоростью

щ. Имеем

 

 

 

P i ^ J i - u J i ,

(XIV,И)

где J t — компонента импульса в

направлении оси Oi.

 

В уравнении (XIV,И) первый члеп выражает общий поток

импульса, второй

член — поток

импульса, вызываемый

только

осредненной скоростью; последний следует вычесть, поскольку пло­

щадка движется вдоль оси со скоростью ut\ разность этих двух слагаемых определяет только турбулептный перенос импульса.

Если импульс отнести к единице объема,

то

J i = РЩ

(XIV, 12)

и его поток будет

 

Я, = (ри?-й?).

(XIV, 13)

Полная скорость щ равна сумме осредненной и пульсационной

составляющих:

 

Щ = щ + иП1.

(XIV, 14)

Подставляя уравнение (XIV,14) в (XIV,13), произведя вычисле­

ния

и имея

в

виду, что ип1 =

0, получаем

 

 

 

 

Pi = püJi= рип tunl.

(XIV,15)

В

главе

VI

было показано,

что

 

 

 

 

=

 

(XIV,16)

где lj — путь перемешивания в направлении оси Oj.

Заменив в уравнении (XIV,15) одну величину иш ее значением

из выражения (XIV, 16), получим

 

p i = *m h isr-

(XIV’17)

Здесь правая часть обозначает сумму трех членов, получающихся при i = const и j = х, у, z.

Сопоставляя уравнения (XIV,17) с (XIV,10) и (XIV,9), видим,

что в общем случае

 

 

е =

= в//.

(XIV,18)

Символ e/у обозначает сумму девяти компонент коэффициента обмена, получающихся при i = х, у, z и / = х, у, z, т. е. коэффициент турбулентного обмена для импульса является тензором второго ранга.

В случае плоско-параллельного потока, движущегося в направле­ нии оси Ох, из выражения (XIV, 17) для потока импульса Рг через плоскость xOz [этот поток импульса определяет т2х и, следовательно,

профиль скорости u(z)] имеем

= z, ^ =

 

=

 

 

P

I --- г

du

i

T du .

---г

du \

---г

du

 

 

z — p y u nzl x

dx

+

u n J'y~ ô ^ ~ ir

u nzl z ' J f J

— Pu n J 'z - Q ^ -

В

этом

частном

случае

коэффициент

обмена в

направлении

оси Oz

 

 

 

 

е// = е22 =

е2

 

 

(XIV,19)

 

 

 

 

 

 

 

можно

считать скалярной

величиной.

 

 

 

Из выражений (XIV,17) и (XIV,12) следует, что турбулентный поток импульса можно выразить через его производную и коэффи­ циент турбулентного обмена:

P ‘ = 8'/ 4 f * (XIV,20)

Сопоставляя выражение (XIV,20) с выражениями для турбу­ лентного диффузионного потока примеси (XIV,7), видим, что они тождественны: поток субстанции определяется через производную некоторой величины, характеризующую ее количество или J) и некоторый коэффициент (DT или е^). Это сходство не случайно: оно определяется тождественностью турбулентного механизма пере­ носа любой субстанции — обменом объемами среды под действием пульсационных скоростей. Поэтому коэффициент турбулентной диффузии может быть выражен аналогично е/у, а именно как осредненное во времени произведение некоторой диффузионной скорости s и некоторого пути перемешивания для концентрации 1С:

DT= —slC9 (XIV,21)

где s — мгновенная скорость пульсационного движения элементар­ ного объема газо-воздушной смеси в турбулентном потоке; — расстояние, на которое перемещается этот объем до потери

им своей индивидуальности.

Вобщем эти величины отличны от пульсационных скоростей ил

ипути перемешивания /, определяющих перенос импульса в тур­

булентном потоке. Однако в большинстве случаев принимают s = = ип, /с = /, т. е.

DT = B. (XIV,22)

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

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

ветственно в s и

1С. Такие газы

называются д и н а м и ч е с к и

а к т и в н ы м и .

С уменьшением

раэности плотностей газа и воз­

духа уменьшается различие между s и ип, и I.

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

Уравнение диффузии. Общее уравнение турбулентной диффузии

имеет вид

,

 

,

 

. dCw

 

дС

дСи

dCv

(XIV,23)

dt

*

дх

'

ду

dz

 

где

 

 

 

 

 

 

Сс -J- с;

и= и + ип;

(XIV,24)

V = V + Và

W= W+ wn;

С, и, и, w — мгновенные значения концентрации и компонент ско­ ростей;

с, и, и, w — их осредненные по

времени значения;

с, иП1 уп, wn — их пульсационные

составляющие.

Подставляя в уравнение (XIV,23) соотношения (XIV,24) и про­ изводя осреднение по времени, оогласно правилам осреднения Рей­

нольдса (см. § 38),

получаем

 

 

дс , - дс

- дс , дс

0 /n 0с \ i д ( п дс \ . д ( п

дс \

dt + иТх

+ V ду

Wdz"" дх

дх ] + ду \Р*У ду ) + dz

dz ) ‘

(XIV.25)

Сравнивая уравнение (XIV,25) с выражениями (XIV,2) и (XIV,7) видим, что три последних слагаемых левой части уравнения (XIV,25) представляют собой сумму производных от конвективных диффу­ зионных потоков по осям координат, а слагаемые правой части — то же, от турбулентных диффузионных потоков. Следовательно, уравнение (XIV,25) выражает з а к о н с о х р а н е н и я м а с с ы п р и д и ф ф у з и и : суммарное количество газа, вносимое в неко­ торый произвольный объем и выносимое из него конвективным и тур­ булентным диффузионным потоком, равно приращению количества

газа в этом объеме. Это приращение оценивается членом -J^-; при

— = 0 (стационарная диффузия) количество вносимого газа равно

количеству

выносимого

газа.

 

 

 

Полное уравнение диффузии получим, если к правой части

формулы (XIV,25)

добавим

производные от молекулярных диффу­

зионных потоков

(XIV,4):

 

 

 

 

дс

. — дс

, Z

дс

.

дс

i [ ( D „ + D j £ ) +

 

-тг +

r v л—

1

dz

 

dt

1

дх

1

ду

 

 

 

+ ^ г [ ( 0 „ + а д | - ] + £ [ < 0 „ + а д & ] .

(XIV,26)

При 4 - = 0 имеем случай с т а ц и о н а р н о й д и ф ф у з и и * dt

Поскольку в уравнение диффузии входит скорость потока, ре­ шать его необходимо с уравнением движения (см. главу VI). Кроме того, для решения необходимо знать зависимость коэффициентов диффузии от скорости и координат, а также начальные и граничные условия (т. е. концентрацию газа при t = 0 и на стенках выработки).

Уравнение неразрывности при диффузии. В § 35 было показано, что при постоянной плотности р уравнение неразрывности (VI,4)

dp

, dpи

,

dpi;

 

dpw

f

dt

' dx

'

dy

'

dz

(XIV,27)

превращается в

du

.

du .

dw

_~

 

 

 

 

dx

+ "d7“T""“âr'~u*

(XIV,28)

При диффузии в общем случае плотность среды переменна. Тем не менее уравнение (XIV,28) справедливо и в этом случае. Действи­ тельно, так как

р = р0-Д р С ,

(XIV,29)

где Ро — плотность воздуха;

 

Ар — разность плотностей воздуха и диффундирующего газа, то, подставив выражение (XIV,29) в (XIV,27) и произведя дифферен­ цирование с учетом уравнений (XIV,24), получим

A dC а / dCa I dCi;

I дСю \ ,

du . dv . dw \ _

- АР - д Г - ^ { - д Г + - д Г + —

)+РоЬг+17+1пг)=0-

Сопоставляя полученное соотношение с (XIV,23), получаем

ди

|

dv

,

dw

(XIV,30)

dx

'

dy

"*

dz

T . e., несмотря на изменение плотности потока воздуха вследствие изменения в нем концентрации диффундирующего газа, с м е с ь в о з д у х а и г а з а в е д е т с е б я к а к н е с ж и м а е м а я ж и д к о с т ь .

Особенности диффузии активных газов. Рассмотрим поток воз­ духа в наклонной выработке, плотность которого вследствие про­ исходящих в нем процессов диффузии изменяется по оси Oz (рис. 123), увеличиваясь от кровли к почве. Выделим в потоке некоторый эле­ ментарный объем со, находящийся в данный момент в слое 1 плот­ ностью PiПусть в следующее мгновение он переместится под дейст­ вием пульсационной скорости в слой 2 плотностью р2, причем р2 < <CPv Допустим, что расстояние между слоями 1 и 2 равно длине прандтлевского пути перемешивания /, т. е. объем со в слое 2 сохра­ няет еще свою индивидуальность и, следовательно, его плотность равна рх. Тогда вследствие разности плотностей смеси в объеме со и в окружающих его объемах появится направленная вертикально вниз (по действию силы тяжести) выталкивающая сила, которая

будет стремиться вернуть объем со в исходное положение. Эта объем­ ная сила, отнесенная к единице объема, составляет

/ = £ (P 1 -P 2)=gd p, (XIV,31)

где g — ускорение силы тяжести.

Очевидно, подобная сила будет действовать на объем со на всем его пути перемещения из слоя 1 в слой 2.

Проекция силы / на ось Oz (/2), действуя против пульсационной скорости, перемещающей объем со, как бы уменьшает ее на некоторую величину Vf.

Рис* 123. Схема действия объемной силы при диффузии активного газа

Величина vf не зависит от направления пульсационной скорости vn: при перемещении объема со вниз под действием ип на столь же малое расстояние /, что и вверх, сила /2 также будет препятствовать его движению, величина ее не изменится, но направлена она будет уже вверх, т. е. по-прежнему против пульсационной скорости ип.

Все приведенные выше рассуждения касались случая, когда плотность потока возрастала от кровли к почве. Очевидно, что при уменьшении плотности от кровли к почве те же объемные силы будут способствовать перемещению объема со под действием пульсационных сил, поскольку их направление будет совпадать с направлением пульсационной скорости ию осуществляющей перемешивание в на­ правлении оси Oz.

Поэтому в общем случае истинная («диффузионная») скорость, под действием которой происходит диффузия активного газа, будет равна алгебраической сумме пульсационной скорости vn в направле­ нии оси Oz и условной скорости vf\

S = vn± Vf. (XIV,32)

229

Знак минус в уравнении (XIV,32) берется при увеличении плот­ ности от кровли к почве, знак плюс — при ее уменьшении в том же направлении.

Выражение (XIV,32) справедливо и в случае, если вместо мгно­ венных скоростей использовать среднеквадратические:

S9 = v'n±v\.

(XIV,33)

При vf > О (плотность потока уменьшается от кровли к почве) объемные силы способствуют перемешиванию газа в направлении, перпендикулярном основному движению, при < 0 (плотность потока увеличивается от кровли к почве) они препятствуют пере­ мешиванию.

Уменьшение интенсивности перемешивания при Vf < 0 является результатом снижения поперечных пульсационных скоростей или — при постоянной осредненной скорости — уменьшения интенсив­ ности турбулентности. Очевидно, что если отрицательная скорость vf имеет ту же абсолютную величину, что и пульсационная уп, т о , согласно выражению (XIV,32), s = 0 и интенсивность турбулент­ ности

5

8 = -=-

и

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

Затухание

турбулентности

характеризуется ч и с л о м Р и ­

ч а р д с о н а

g

dp

 

 

 

 

R i = ^ i - c o s p ,

(XIV,34)

где — угол

наклона выработки

к горизонту.

объемных сил

Число Ri

пропорционально

отношению работы

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

Можно показать, чтс

Штрихами здесь обозначены среднеквадратические значения ско­ ростей.

Как отмечалось ранее, при = vu происходит затухание турбу­ лентности. Этому соответствует критическое число Ричардсона

Шк»= (* £ У

(XIV’35)