
книги из ГПНТБ / Клебанов, Ф. С. Аэродинамическое управление газовым режимом в шахтных вентиляционных сетях
.pdfГлава 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) |
|
4Кт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