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

книги из ГПНТБ / Клебанов, Ф. С. Аэродинамическое управление газовым режимом в шахтных вентиляционных сетях

.pdf
Скачиваний:
17
Добавлен:
19.10.2023
Размер:
4.74 Mб
Скачать

Глава III

ГАЗОВАЯ ДИНАМИКА ЭЛЕМЕНТОВ ВЕНТИЛЯЦИОННЫХ СЕТЕЙ

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

Причиной газовой нестаниионарности в выработках служат: I ) нестационарность источников газа и 2) изменение условий проветривания (аэродинамическая нестационарность).

Нестационарность источников газа определяется отдельным или сов­ местным изменением во времени интенсивности источников, их перемещением или изменением геометрии источников.

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

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

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

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

61

1. Обзор известных результатов

Необходимость исследования динамики распространения примесей в трубных потоках текучего возникла давно. Это было связано с вопро­ сами гемодинамики (распределение примеси в потоке крови), с опре­ делением скорости воды в' водопрЬводных магистралях и нефтепрово­ дах большой длины, с разработкой метода экспериментального опреде­ ления коэффициента молекулярной диффузии по изменениям распростра­ нения исследуемой примеси в трубе [ 2], Однако достаточно законченное

решение общих задач,

относящихся

к указанным областям,

было

получено сравнительно

недавно [2, 3 ].

Первоначально был рсследован

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

[2].

 

Если в поток жидкости или газа ввести какую-либо примесь,

раст­

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

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

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

Решение данной задачи приведено в работе [2 ] на основе уравне­

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

 

 

D( £ с

1

Эс

с)2с

dt

дс

(Ш,1)

Эг2

г

дг

дх*

i - i,2) дх

 

где D - коэффициент молекулярной диффузии, независимый по пред­ положению от концентрации; а - радиус трубы; uQ - скорость на оси трубы.

62

 

Обозначая

z =■г/а

 

 

 

 

 

и учитывая то, что

 

 

 

 

 

 

 

д2с

 

д2с

1

дс

 

 

 

 

 

 

'

 

'

+ ------- .

 

 

 

 

 

 

дх2

 

Эг2

г

дг

 

 

 

 

 

из формулы

(111,1) получаем

 

 

 

 

 

д2с

1 дс

а2 дс

а2и0

о. Эе

 

(Ш,2)

 

дъ2

z

3z

D

iJt

D

Эх

 

 

 

 

 

 

Граничное условие для равенства (111,2)

выражает тот факт,

что

стенка трубы непроницаема

 

 

 

 

 

 

 

 

 

 

 

 

при z = 1.

(Ш,3)

 

Решение уравнения (III,2)

дс

дается в виде

 

 

при ------ О

 

 

 

 

 

 

 

 

 

дх

 

 

 

 

 

- a t

 

1/2п-1/2 .

 

 

(Ш ,4)

 

с = е

 

I (aza

L)

),

 

 

где IQ(aza

/

D

)

- функция Бесселя первого рода нулевого поряд­

ка;

a

- постоянная времени.

 

 

 

 

 

Из граничного условия

(III,3) следует

 

 

 

 

 

1

/2

-1 /2

0.

 

 

 

(Ill,5)

 

Ij (аа

^ D

1/*)=

 

 

 

 

Корень уравнения

(III,5),

соответствующий самому нижнему зна­

чению параметра а,

есть

 

 

 

 

 

 

aaj

D

 

= 3,8,

 

 

 

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а2

 

 

 

 

(Ш,6)

 

 

 

 

3,82D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если в какой-либо момент времени примесь вводится на длине тру­

бы

L;

то время, необходимое для того, чтобы сделать заметным изме­

нение концентрации в результате конвекции,

имеет порядок L/uq,

и,

следовательно,

в соответствии с формулой (111,6)

определяющим усло-

63

вием для возможности учета в механизме переноса примеси только молекулярной диффузии является неравенство

(Ш.7)

14,4 D

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

<Э2с

дс

 

 

К —

 

 

(III,8)

 

=

<9xj

dt

 

где К -

 

эффективный коэффициент диффузии, определяемый по формуле

а

2

2

 

 

 

и~.

 

 

К =

 

О

 

(Ш,9)

 

 

 

1920

 

 

Координата

х^

 

 

 

uol

 

 

Это означает, что рассматривается перенос примеси через плоскость которая движется с постоянной скоростью, равной средней скорости

ламинарного потока в трубе.

 

 

коор­

Этим самым в рассмотрение вводится подвижная система

динат.

 

 

 

 

случаев распростра­

Известны решения уравнения (111,2) дця двух

нения примеси:

 

 

 

 

 

1)

примесь массы М в момент времени t = 0

концентрируется в

точке

х = О,

 

 

 

 

 

2) примесь с равномерной концентрацией с0 вводится в трубу с

постоянной скоростью в точке

х = 0 начиная с момента времени

t= О

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

с - V ) .

 

 

 

для 1-го случая

 

 

Решение уравнения (Ш,8)

 

 

 

1 м _2 -3/2

К

-1/2

х ,

 

(Ш Л О )

с =—Ма ^тг

ехр ( ------- ).

 

 

2

 

 

4Kt

 

 

64

Для 2-го случая

;r0 '

i l l

~ “

f ( W

) " р" 1,1 >0;

( " U l )

где

 

z

 

о

 

 

 

 

 

 

 

,

2

г

е

-z2 j'

 

 

erf = ------/

 

dz.

 

 

В 1-м случае длина Ц , которая содержит 90% примеси, определя­ ется из выражения

 

1

L, К

-1 /2

—1/2

erf ( —

' t

) = 0,9,

 

4

1

 

 

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

ц >

4,65К 1/211/2-

("1.12)

Если

с тах -

максимальная концентрация в точке

то концентрация на концах интервала Ц = 0,3126 стах-

Во 2-м случае длина L2 перемещающейся зоны примеси, в которой

концентрация изменяется от 0,9 сс

до 0,1 cQ, находится из равенства

.гГ

( i

L2K -1/2r ,/2

) = 0,8;

 

 

А

А.

 

 

откуда

 

 

 

 

ч

= 3,62K1/2t1/2

 

(III,13)

Результаты теоретического и экспериментального исследования про­ цесса распределения примеси в турбулентном потоке в трубах описаны в работе [ 3 ] и коротко сводятся к следующему.

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

. К = 10,1 av*,

(111,14)

ИЛИ

 

К = 7,14avyl/2,

(ШД5)

5 820.

65

где u - радиус трубы; v - средняя скорость текучего в трубе; v* - динамическая скорость; у - коэффициент Шези-Дарси

Г 'Т

V ,/ ~ •

Формулы (Ш,14) и (111,15) применимы как для гладких, так и для шероховатых труб.

2.Распределение концентрации примеси, первоначально сконцентри

ванной в точке х = О (при t = 0),

описывается кривой Гаусса распре­

деления вероятностей

 

с - A t" 1/ 2 exp [ -(x -vt)2

(111,16)

4Kt J

 

где A - постоянная, зависящая от общего количества диффундирующей примеси.

Выражение (111,16) является решением дифференциального уравнения

в форме (111,8) для начальных условий, соответствующих импульсу при­ меси.

3. Максимальная концентрация импульса примеси переносится со средней скоростью потока, т.е. центр распределения примеси в каждый данный момент находится в точке

х= vt.

Отклонение от симметрии выражения (III,16) незначительно при уоловии

х40v*

— » ----- . a v

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

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

5. Числовое значение эффективного коэффициента диффузии может быть определено из эксперимента по формуле

2

Х 1/2

= In 2,

(111,17)

т1/2

 

 

66

где 2tj/'2 -

период времени, в течение которого концентрация приме­

си в данном

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

центрации для этого пункта.

Последнюю формулу можно представить также в виде

V3T2

К = •___Ш_

'

(HI.18)

4х.1п 2

 

где х - длина трубы от точки ввода примеси до точки наблюдения.

6. Расчет распределения концентрации в случае, когда одно текучее следует за другим, может также производиться на основе решения диф­ ференциального уравнения (III,8) с использованием эффективного коэф­ фициента диффузии, определяемого по формуле (111,14).

Решением указанного уравнения является выражение (III,16), где

Xj = х —vt.

7. Длина L вдоль трубы, на которой концентрация примеси изменя­ ется в пределах 0,99 < с < 0,01, определяется выражением

о

V*

(111,19)

L2

= 437 аХ — ,

 

v

 

т.е. T~Y1/2L —X

Экспериментальная проверка основных положений теории Тэйлора применительно к условиям шахтных выработок была предпринята Ходкинсоном и Личем [4 ]. Их выводы сводятся к следующему.

Максимальная концентрация примеси (газ, витающая пыль) распро­ страняется со средней скоростью воздуха в полном согласии с теори­ ей Тэйлора.

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

Длина импульса примеси - величина постоянная для данного рассто­ яния от источника при любых скоростях воздуха.

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

L = 0,42 х.

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

L = 5,4 х1/2

Это совпадает с одним из выводов теории Тэйлора (III, 19).

67

2.Уравнение процесса распространения

нвыноса газообразной примеси потоком воздуха

в Длинной прямой выработке

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

1.В длинной прямой выработке движется воздух со средней ско

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

ном потоке устанавливается на уровне

С2 (первоначально считаем, что

концентрация была равна Cj = 0,

т.е,

воздух не содержал примеси

газа).

 

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

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

16, а)

 

при

х < 0

с(х,0) = <

(III,20)

при х > 0.

с2=0

Ф )

г * 2 г=а

с(х)

> 2

\т>о

* 2

* /

 

М

И Ш

-

тс(х)

f=Q

с/х)

Уг

_______ Е ____

 

Ш /Ш М-

с/х)

 

сг=0

г = в

С1

' Ф )

Т>0

 

XT*

Рис.16. Графическое изображение ? различных краевых усло-

Ш Ш И М ВИЙ ИЛЯ УРавнения Рао“ пространения газа в длин-

ной прямой выработке

4*2

Ф )

Г*2

68

В действительности в рассматриваемой схеме при х <0 концентра­ ция газа в воздушном потоке равна нулю. Однако целесообразно при­ нять при х<0 концентрацию равной с2, так как очевидно, что данный случай можно интерпретировать также в ином виде, а именно:-не как начало действия в момент т= 0 источника газа, сосредоточенного в пункте х= 0, а как вытеснение из выработки чистого воздуха пото­ ком газовоздушной смеси с концентрацией

Если в выработке длительное время действовал точечный источник

газа,

а затем в момент т= 0

его действие прекратилось, то возника­

ет ситуация обратная той, что

отображается условиями (111,20), т.е.

чистый воздух начинает постепенно вытеснять загазованный воздух

(рис. 16,6).'

 

 

 

Краевые условия при этом принимают вид

 

I" 0

при

х < 0

 

с (х, 0) = ^

при

х > 0.

(111,21)

 

I Cj

2.

В общем случае, когда начальный уровень концентрации гааа в

выработке не равен нулю, т.е. по выработке движется частично зага­ зованный воздух, вместо формул (111,20) и (111,21) имеем (рис. 16, в, г)

С2 при

х < 0

с(х,0)

( 111, 22)

Cj при

х > 0. ■

3. В выработке в точке х=0 действует точечный источник газа постоянной интенсивности. Первоначальная объемная скорость воздуха Qj в момент т= 0 резко (неплавно) изменяется до величины Q2-

Краевые условия для данной ситуации определяются также выраже­ нием (111,22).

4. Интенсивность точечного источника газа резко (неплавно) изме­ няется от G] до G2 при постоянстве объемной скорости воздуха в вы­ работке. Данный случай также имеет краевые условия в виде формулы

(111,22).

Уравнение, которое описывает в подвижной системе координат про-, цесс рассеивания и выноса примеси в выработке под совместным влия­ нием молекулярной и турбулентной диффузии и конвективного переноса, для краевых условий (111,20) имеет вид

Вс

д2 с

(111,23)

--------К -----

at

ах

 

Х1 = х —vt; с* = — ,

с2

где v - скорость движения воздуха в выработке.

69

Из уравнения (Ш,23) следует формула для определения концентрации газа, рассеиваемого и выносимого потоком воздуха в длинной прямой выработке:

1

х—vt

(HI,24)

c^x .t) i [ l

_Ф ( --------- )];

 

2^Kt

 

 

X—vt

(111,25)

c(x,t) = — [ 1 —Ф (

2

 

 

Из формулы (III,25) следует, что граница раздела загазованного и чистого врздуха перемещается со средней скоростью потока v и явля­ ется центром зоны перемешивания. Загазованный воздух проникает в глубь чистого воздуха на расстояние от центра зоны перемешивания, равное половине всей длины зоны. Концентрация газа впереди центра зоны перемешивания изменяется от 0 до с2/2, а позади зоны от с2/2 до с2- С увеличением времени движения потока длина зоны переме­ шивания возрастает, но пределы изменения концентрации газа в зоне остаются, очевидно, теми же, т.е. С2 > с >0.

Из выражения (111,25) можно получить формулу для распределения концентрации газа вдоль выработки для некоторого фиксированного мо­ мента времени. Если же зафиксировать како&-пибо пункт вдоль выра­ ботки, то это же выражение дает формулу для изменения концентрации газа во времени в данном пункте.

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

c(x,t)

'1

X—vt

 

(111,26)

[ 1

+Ф (

)].

 

 

 

2 / к 7

 

 

Для начальных условий (Ш,22) при с2 > Cj (рис. 16, в)

имеем

 

 

п _ ® (

■ i c ^ - n

+ c ,.

(111.27)

 

2

 

2 ^

 

 

Если же

С2 <

(рис. 16,г), то

 

 

C( x ,t) - .^ L 1 5

[1 + ф ( - ^ L ) ] +c2.

(Ill,28)

Из выражений (111,27)

и (111,28) видно, что граница раздела воздуха

с концентрацией газа с2

и воздуха

с концентрацией газа

Cj по мере

70

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