Решение: масса газа, переносимая в результате диффузии, определяется законом Фика:
[ m=-Dcdot frac{Delta rho }{Delta x}cdot Scdot t, ]
где D – коэффициент диффузии; ∆ρ/∆x – градиент плотности, т.е. изменение плотности, приходящееся на 1 м толщины слоя почвы; S – площадь поверхности слоя; t – длительность диффузии.
Выразим искомый коэффициент диффузии и рассчитаем его. При этом нужно учесть следующее: значение градиента плотности – отрицательное число, что соответствует сущности процесса диффузии- зависимость плотности от расстояния в направлении движения диффундирующей массы выражается убывающей функцией, градиент которой — отрицательная величина.
[ D=-frac{m}{frac{Delta rho }{Delta x}cdot Scdot t},text{ }D=-frac{14,5}{-1,45cdot 100cdot 24cdot 3600}=1,16cdot {{10}^{-6}}. ]
Ответ: 0,0116 см2/с.
В выражениях для диффузионных газовых
потоков коэффициенты молекулярной
и турбулентной диффузии являются
единственными параметрами, учитывающими
свойства среды. Естественно, что эти
величины имеют сложный характер, и их
определение — одна из важных задач
теории диффузионных процессов.
Коэффициенты молекулярной
диффузии. Для
газов со сходными молекулами (имеющими
почти равные массы и эффективные сечения)
Максвелл получил следующее выражение
для коэффициента молекулярной
диффузии:
,
где
—
длина свободного пробега молекул;vм
— скорость их теплового
движения; черта означает среднее значение
величины. При нормальных условиях
имеет порядок 10-5см,
vм
= 10-4
÷10-5
см/с.
В силу статистической
однородности молекулярного движения
величины
и
,
а, следовательно, и коэффициент
молекулярной диффузии не зависят от
направления. Коэффициент молекулярной
диффузии слабо зависит от содержания
диффундирующего газа. С увеличением
температуры он возрастает пропорциональноТ1+a,
где Т
— абсолютная температура
среды, а — коэффициент,
изменяющийся от 0,5 до 1. С увеличением
давления коэффициент уменьшается в
обратно пропорциональной зависимости.
Выше отмечалось, что в
шахтных условиях молекулярная диффузия
имеет подчиненное значение в процессе
переноса газов. Кроме того, изменения
содержания газов, температуры и давления
воздуха в активно вентилируемых горных
выработках относительно невелики.
Поэтому при решении задач газопереноса
в шахтах можно принимать Dм
= const.
Следует иметь в виду, что коэффициент
молекулярной диффузии газа в среду
равен коэффициенту молекулярной диффузии
среды в этот газ. Средние значения
коэффициентов молекулярной диффузии
некоторых газов приведены ниже.
Газ
Температура, °С Коэффициент
диффузии, см2/с
Аммиак в воздухе 0
0,217
Водород в воздухе —
0,634
Метан в воздухе —
0,196
Оксид углерода в воздухе
— 0,129-0,138
Углекислый газ в воздухе 0
0,142
Коэффициенты турбулентной
диффузии. В
теории турбулентности коэффициент
турбулентной (или вихревой) диффузии
вводится как некоторый коэффициент
пропорциональности. При этом для его
выражения используют три принципиально
различных подхода.
В первом способе
коэффициент
турбулентной диффузии
определяют, следуя Буссинеску, как
коэффициент пропорциональности
между потоком газа и градиентом содержания
в соответствии с формулой (6.11) — jт
= — Dт·grad
c.
Известно, что произведение
вектора, каким является в формуле (6.11)
градиент содержания, на некоторую
величину [в выражении (6.11) ею является
коэффициент турбулентной диффузии Dт]
может дать вектор [в формуле (6.11) это
вектор потока газа] лишь в случае, если
эта величина является скаляром или
тензором. Коэффициент турбулентной
диффузии не может быть скаляром в связи
с тем, что в случае равенства производных
от содержания по направлениям компоненты
газовых потоков по этим направлениям
также были бы равны, что в условиях
существенно неоднородного и неизотропного
турбулентного воздушного потока в
выработках невозможно вследствие
различия компонент пульсационных
скоростей.
Таким образом, остается предположить,
что коэффициент турбулентной диффузии
в горной выработке — тензор. Можно
показать, что в условиях неоднородной
и неизотропной турбулентности коэффициент
турбулентной диффузии — тензор второго
ранга. Тогда компоненты газового потока
будут иметь следующее выражение:
(6.17)
Здесь
(6.18)
(i,j
= х,у,z)
является тензором
коэффициентов турбулентной диффузии
второго ранга с компонентами Dтxx,
Dтxу,
…, Dтzz.
Выражение (6.17) может быть записано в
свернутом виде
.
(6.19)
где правая часть представляет
собой сумму трех значений
,
получающихся, если
фиксировать i,
а j
придавать последовательно
значения х, у, z
(суммирование по
двойному индексу).
Выражение (6.19) обычно
упрощают, принимая, что оси Ох,
Оу, Оz,
являются главными
осями тензора. Если тензор
симметричный, то
и, следовательно, коэффициент турбулентной
диффузии определяется только диагональными
компонентамиDтxx,Dтyу,Dтzz.
Для однородной и изотропной турбулентности
имеет место сферическая симметрия
газовых потоков. Следовательно,
.6.20)
В этом частном случае
может рассматриваться как скаляр.
В выражении (6.11) векторы jт
и gradс
коллинеарны.
Следовательно, согласно определению,
направление вектора gradс
является главным направлением
тензора, а ось координат, соответствующая
ему, — главной осью. Нахождение главных
осей тензора коэффициентов диффузии
для выработки — в ряде случаев задача
неопределенная, так как для этого
необходимо знать поверхности равных
содержаний в потоке, т.е. поле содержаний,
что обычно является конечной задачей
исследований. Лишь в простых случаях
диффузии главные направления могут
быть определены достаточно просто.
Например, при газовыделении с одной
стенки gradс
с некоторым приближением
можно принять нормальным к этому боку
и, следовательно, главные оси тензора
будут направлены вдоль потока воздух
и перпендикулярно к нему. В более сложных
случаях главные оси тензора могут
иметь и другие направления.
Следует отметить, что
принятие тензора Dтy
симметричным для случая движения воздуха
в горной выработке является также
определенным допущением. Для
неоднородного и неизотропного
турбулентного потока, каким является
вентиляционный поток в выработке,
тензор коэффициентов турбулентной
диффузии будет несимметричным. Ниже
отмечается, что компоненты тензора
коэффициентов турбулентной диффузии
могут быть выражены через усредненное
произведение (корреляцию) мгновенных
значений пульсационной скорости иni
и пути перемешивания
для содержания
(здесьi,j
= х, у, z,
иni
= ип;
uпу
=vn
unу
=wn).
Для симметричного тензора должны
соблюдаться равенства
,
что приводит к соотношениям.
Однако для неизотропных вентиляционных
потоков корреляциянесимметрична относительноi
и j,
а это не отвечает приведенным равенствам.
Несимметричность тензора коэффициентов
турбулентной диффузии для шахтных
вентиляционных потоков косвенно
доказывается фактором различной
интенсивности турбулентной диффузии
в разных направлениях.
Отмеченные приближения, которые применяют
при решении практических задач шахтной
газовой динамики, в настоящее время не
имеют оценки. Применительно к условиям
диффузии в приземном слое атмосферы
погрешности незначительны (в некоторых
случаях они составляют 15-20 %). Однако
степень анизотропности шахтных
вентиляционных потоков значительно
выше атмосферных, что может привести
к необходимости учета факта несимметричности
тензора диффузии.
Второй способ
определения
коэффициента турбулентной диффузии
основан на использовании теории Прандтля
о пути перемешивания, согласно которой
компоненты потока газа можно определять
как сумму трех слагаемых:
.
(6.21)
Здесь, подобно тому, как это
было принято в выражении (6.19) —
,
суммирование производится по двойному
индексу (j);
i
= j
= х,у,z;
;
Lc
— путь перемешивания для содержания.
Из выражения
следует, что коэффициент турбулентной
диффузии является тензором второго
ранга
(6.22)
определяемым девятью
компонентами —
Сопоставляя методы выражения
коэффициента турбулентной диффузии по
Буссинеску и Прандтлю, видим, что в
первом случае коэффициент турбулентной
диффузии остается неопределенным, во
втором — определяется через характеристики
турбулентного движения ().
В случае плоского потока
()
коэффициент турбулентной диффузии
в поперечном к основному движению
направлении определяется из выражения
(6.21):
,
.
В случае изотропной
турбулентности можно принять Lcx
= Lсу,
что приводит к равенству
,
(6.23)
т.е. в этом частном случае коэффициент
турбулентной диффузии является скаляром.
Если в уравнении (6.23) vn
выразить по Прандтлю через путь
перемешивания для импульса L,
то для плоского потока
получим выражение
,
(6.24)
где
среднее квадратичное значениеvп
;
а1
— коэффициент
пропорциональности между un
и vп.
Если принять, что
L
/ Lс
= а2
= сопst,
(6.25)
то
.
(6.26)
Величина
для случая диффузии газа является
аналогом пути перемешивания для импульса
по Карману [21] (не тождественному
прандтлевскому пути перемешивания).
Из уравнения
видно, что, имея какие-либо гипотезы
относительно величинlС(),
можно, измеряя
в потоке, определить коэффициент
турбулентной диффузии. Наиболее простым
допущением является отождествлениеlс
путем перемешивания
для импульса l;
во многих случаях такое приближение
дает вполне удовлетворительные
результаты.
Следующим шагом в этом
направлении является принятие
пропорциональности между lс
и l;
значение коэффициента пропорциональности
между ними зависит от свойств
диффундирующего газа, разности содержаний
газа в диффундирующем объеме и в среде.
По имеющимся сведениям, этот коэффициент
больше 1; для азота он равен ~,
для гелия ~.
Имеются попытки оценитьlс
через l
и критерий Ричардсона, характеризующий
затухание турбулентности под действием
объемных (гравитационных) сил при
диффузии активного газа.
Наконец, третий
способ определения коэффициента
турбулентной диффузии
основан на представлении процесса
диффузии как случайного движения жидких
частиц, первоначально сконцентрированных
в некоторой области. Бэтчелор показал,
что и в этом случае коэффициент
турбулентной диффузии является тензором
второго ранга. Запись его (для случая
однородной турбулентности), однако,
имеет иной вид:
,
где уi,
yj
— лагранжевы координаты
жидкой частицы, величины случайные,
являющиеся функцией времени.
Представление коэффициента
турбулентной диффузии в виде тензора
имеет в основном теоретическое значение.
В настоящее время практически ничего
неизвестно о недиагональных компонентах
этого тензора. Изученные в какой-либо
степени компоненты тензора диффузии —
это диагональные компоненты Dтxx,Dтyу,Dтzz,
которые в дальнейшем
и будут рассматриваться. Для простоты
написания обозначим Dтxx
= Dтх
и т.д.
Необходимо отметить, что в
общем случае коэффициент турбулентной
диффузии является функцией координат.
Это можно видеть, например, из уравнения
(6.24), где величины dи/dу,
v´n,
lс
для потоков в
горных выработках являются функциями
поперечных координат [23], а в некоторых
случаях (изменение сечения по длине
выработки, свободные струи) — и продольной
координаты. Эти же величины являются и
функциями скорости потока (точнее, числа
Рейнольдса — Rе
потока), что говорит о существовании
зависимости коэффициента турбулентной
диффузии и от числа Rе.
Данные о коэффициентах турбулентной
диффузии в горных выработках
немногочисленны, что в значительной
степени объясняется техническими
трудностями их измерений. Имеющиеся
сведения частично основываются на
данных о коэффициенте турбулентного
обмена для импульса и предположении о
пропорциональности ему коэффициента
турбулентной диффузии [24].
К.М. Тумаковой были установлены
автомодельность поперечных составляющих
относительного коэффициента турбулентной
диффузии:
;
— средняя скорость потока; α
– коэффициент аэродинамического
сопротивления;
— плотность потока; Н
— высота выработки)
по числу Рейнольдса, начиная от Rе
= 13600, а также равенство вертикальной и
горизонтальной поперечных составляющих
коэффициента диффузии.
Их значения в ядре
потока равнялись 0,02, а на расстоянии
0,13H
и 0,8H
от кровли — 0,03.
В ряде случаев хорошие результаты
получаются, если использовать средние
по высоте (ширине) выработки значения
коэффициентов турбулентной диффузии.
Коэффициент турбулентной диффузии
может быть рассчитан по характеристике
рассеивания газа. Для случая однородной
и изотропной турбулентности в
равномерном потоке воздуха (без градиента
скорости) распределение содержания
газа в газовом факеле за источником
газовыделения описывается гауссовой
кривой ошибок:
,
(6.27)
где с —
содержание газа в точке с координатами
х,у;z
— расстояние от
источника вниз по потоку; у
— расстояние от точки,
соответствующей максимальному
содержанию газа сmах
в плоскости х =
соnst,
измеряемое в направлении, перпендикулярном
направлению движения воздуха; и
— скорость потока
воздуха.
Если в формуле (6.27) с
выразить как часть
сmах,
то из нее можно определить Dт.
Например, полагая с
= сmах/2,
получим
Dт=,(6.28)
где
—расстояние от оси
газового факела до точки в его поперечном
сечении, в которой с
= сmах/2.
Все величины в уравнении
(6.27) поддаются прямому измерению: и
и х
измеряют непосредственно
на месте эксперимента,
—по графику зависимости
с(у), построенному
на основании измерения содержания на
расстоянии от источника, равном х.
Поскольку выражение (6.27) справедливо
лишь для однородной и изотропной
турбулентности, то в силу равенства
(6.20) по нему определяют диагональные
члены тензора коэффициентов диффузии,
не зависящие от координат.
Известно, что турбулентность
шахтных вентиляционных потоков
неизотропна; ее можно считать однородной
лишь в направлении основного течения
(при неизменных форме сечения,
шероховатости стен и расходе воздуха).
Поэтому для шахтных условий выражение
(65.28) дает, во-первых, неточные значения
Dт
и, во-вторых, лишь
некоторые средние значения поперечной
компоненты тензора Dтy.
Погрешности будут возрастать по мере
приближения источника газа от оси потока
к стенке, так как при этом источник
попадает в области все большего градиента
скорости, т.е. все большей анизотропии
турбулентности.
Учитывая экспериментальное
подтверждение аналогии Рейнольдса для
процессов переноса импульса и пассивной
примеси при течениях в пристеночной
области, коэффициенты диффузии шахтных
вентиляционных потоков при диффузии
пассивных газов в первом приближении
можно принимать равными коэффициенту
турбулентного обмена для импульса. Для
чисел Рейнольдса от 1,25·104
до 3,72·104
относительные значения последних для
штрекообразной выработки прямоугольного
сечения, закрепленной рамной крепью
из круглого леса с продольным калибром
7,5, относительной шероховатостью в
направлении вертикальной оси 8,9,
горизонтальной (перпендикулярной
основному движению) 8,4 приведены на
графиках рис. 6.1 и 6.2, где y
— координата, перпендикулярная бокам
выработки, z
— кровле и почве. Пересчет относительных
значений турбулентного обмена импульса
в абсолютные производится по формуле
ε =
ε*v*D
, где D
— характерный линейный
размер потока (например, диаметр).
Приведенные на графиках данные
соответствуют средним по сечению
абсолютным значениям коэффициентов
турбулентного обмена для импульса εу
и εz,
порядка 5·10-3
м2/с
при средней скорости воздуха в выработке
uср
=1 м/с, коэффициенте трения α
= 15·10-3
Н·с2/м4,
плотности воздуха
= = 1,22 кг/м3,
диаметре выработки D
= 2,5 м.
Рис. 6.1. Зависимость
от у* = = у/Н (Н — высота
выработки)
Рис. 6.2. Зависимость
от z*
= = z/В
(В — ширина
выработки)
Значения компоненты Dтy·103
(м2/с),
полученные для некоторых видов выработок,
приведены ниже:
модель штрекообразной выработки,
площадь поперечного сечения13,4×14,2 см,
средняя скорость воздушной струи 0,25
м/с……………………………………1,1
квершлаг, закрепленный
анкерами, площадь поперечного сечения
24,5 м2,
скорость воздушной струи 0,5-1,2
м/с………………………………………………..2,4÷4,1
то же, площадь поперечного
сечения 23 м2,
скорость воздушной струи 1,1 м/с……..6,8
квершлаг без крепления,
сечение сводчатое, площадь поперечного
сечения 11,8 м2,
скорость воздушной струи 1,7
м/с……………………………………………….5,1
то же, площадь поперечного
сечения 7,5 м2,
скорость воздушной струи 0,8 м/с……1,8
штрек без крепления,
сечение сводчатое, площадь поперечного
сечения 10 м2,
скорость воздушной струи 0,27
м/с………………………………………………………0,8
Для расчета продольных Dтx
и поперечных Dтy
компонент коэффициента турбулентной
диффузии метана в воздухе можно
использовать приведенные ниже формулы.
Для штрекообразных выработок
;
(6.29)
,
(6.30)
где
,
а число Рейнольдса не рассчитывается
по
.
Для элемента S
(м2)
поперечного сечения штрекообразной
выработки при средней скорости по
площади элемента и’ср
(м/с):
.
(6.31)
Для круглых гладких и шероховатых труб
,
(6.32)
где R
— радиус трубы.
Для широкого прямого канала
.
(6.33)
Для диффузии углекислого газа в воздухе
,
(6.34)
где k:
= 3,96·10-4
м.
В формулах (6.29)-(6.33) использованы
следующие обозначения:
Н- высота выработки, м;
—
динамическая скорость, м/с;
uср
— средняя скорость воздушной струи, м/с;
α —
коэффициент аэродинамического
сопротивления, Н·с2/м4;
-плотность
воздуха, кг/м3;
ν — кинематический коэффициент
вязкости, м2/с;
S
— площадь поперечного
сечения выработки, м2.
По этим формулам для некоторых
средних условий (uср
= 1 м/с; Н= 2,5
м,
=
0,1 м/с; R
= 1 м) значения компонент
Dтx,
Dтy
составляют порядка 10-3
м /с.
Коэффициент турбулентной
диффузии D
характеризует рассеивание газа в
потоке за счет работы турбулентных
пульсаций. В ряде случаев на перемещения
диффундирующего газа налагаются более
сильные движения, вызываемые наличием
сдвига (градиента) скорости потока.
Именно к таким потокам — «потокам со
сдвигом» — относятся шахтные
вентиляционные потоки.
В 1951 г. В.Н. Воронин показал,
что при движении газового облака по
выработке его продольная деформация
определяется профилем скоростей. В
1953 г. Дж. Тэйлор опубликовал решение
задачи продольной турбулентной
диффузии примеси от мгновенного источника
в круглой трубе. Им было показано, что
продольное рассеивание примеси,
вызываемое градиентом скорости,
существенно больше, чем рассеивание,
вызываемое турбулентными пульсациями
скорости. Дж. Тэйлор предложил оценивать
суммарный эффект продольного рассеивания
примеси относительно плоскости,
движущейся со средней скоростью
потока, коэффициентом, который получил
название эффективного коэффициента
диффузии DЭ:
DЭ
= Dг
= Dтх,
(6.35)
где Dг
— коэффициент диффузии, вызываемый
наличием градиента скорости (коэффициент
градиентной диффузии, или коэффициент
дисперсии, по Дж.Тэйлору виртуальный,
эффективный, действительный коэффициент
диффузии); Dтх
— коэффициент
продольной турбулентной диффузии.
Газовый поток, вызываемый градиентной
диффузией, определяется выражением
,
(6.36)
где
,
—отклонения
соответственно содержания и скорости
потока в точке от их средних по поперечному
сечению значений;
;
;
с, и
— соответственно
усредненное по времени содержание и
скорость потока в точке;
—
их средние по поперечному сечению
значения.
Коэффициент Dэ
может быть определен в точке или быть
усредненным.
Исследования С.П. Грекова и А.Е. Калюсского
позволили получить следующее выражение
для эффективного коэффициента диффузии
штрекообразной выработки:
;
(6.37)
по И.Ф. Ярембашу
.
(6.38)
Здесь v
— кинематический
коэффициент вязкости воздуха, м2/с;
uср,
— средняя скорость воздушного потока,
м/с; D
— диаметр выработки,
м; α
— коэффициент аэродинамического
сопротивления выработки, Н·с2/м4;
S
— площадь поперечного сечения выработки,
м2;
— плотность воздуха, кг/м3.
К.Ю. Лайгна и Э.А. Поттер в
своих последних работах
дают следующее выражение для среднего
по поперечному сечению эффективного
коэффициента диффузии:
,
(6.39)
где,
а число Рейнольдса рассчитывают по
.
Значение Dэ
можно определить также по графикам,
представленным на рис. 6.3.
Рис. 6.3. Графики к определению
эффективного коэффициента турбулентной
диффузии
Расчеты по приведенным
формулам и графикам дают значения Dэ,
порядка нескольких м2/с.
Дж. Тэйлор и К.Ю. Лайгна,
исследуя влияние изогнутости канала
на коэффициент диффузии, сделали вывод,
что этот фактор может увеличивать
Dэ
до двух раз.
Для изучения газодинамических
процессов при действии свободных
струй В.Н. Воронин применил коэффициент
турбулентной диффузии kт,
определив его как отношение среднего
содержания газа в поперечном сечении
ядра постоянной массы свободной струи
ся
к среднему содержанию
на ее границе сгр:
.
(6.40)
Значения kт
зависят от условий распространения
свободной струи и изменяются от 0,3 до
0,9.
Основываясь на подобии полей скоростей
и содержаний в ядре постоянной массы
свободной струи, В.Н. Воронин получил
следующие выражения для коэффициента
турбулентной диффузии чистых (не
содержащих газа в начальном сечении)
свободных струй:
для основного участка круглой струи
kт
= 1÷1,84А;
(6.41)
для основного участка плоской струи
kт
= 1÷1,44А´;.
(6.42)
В приведенных формулах
(6.43)
,
(6.44)
где RЯ
— радиус ядра постоянной
массы; и — скорость
в точке с координатами
х,
у; и0
— осевая скорость; φя
— относительная координата границы ядра
постоянной массы;
;
(6.45)
а —
коэффициент структуры свободной струи,
зависящий от начальных турбулентности
и профиля скорости (по В.Н. Воронину, для
круглой струи а =
0,044÷0,053, для плоской а
= 0,09÷0,12).
Значения коэффициентов турбулентной
диффузии, рассчитанные по приведенным
формулам, даны на рис. 6.4.
Приведенные выражения
справедливы для свободных струй в
неограниченном пространстве. В условиях
горных выработок свободные струи часто
распространяются в ограниченных объемах,
при этом воздухообмен между струей и
окружающим воздухом определяется не
только структурой струи, но и структурой
воздушных потоков в окружающей ее среде,
которая в свою очередь зависит от
геометрии ограничивающих поверхностей
и их шероховатости и в общем случае
отлична от таковой в неограниченных
объемах. В результате коэффициенты
турбулентной диффузии струй в ограниченных
пространствах отличаются от таковых в
неограниченных пространствах. Впервые
это было отмечено Ю.М. Первовым, который
предложил учитывать его соответствующим
изменением коэффициента структуры
а.
а
Рис. 6.4. Зависимость kт
от
для круглой (а) и аl/b0
для плоской (б) струй
(l
— длина струи, S
— площадь ее начального сечения)
С учетом п
— отношения ширины
камеры к ширине выработки, подводящей
воздух, согласно Ю.М. Первову:
для струи, выходящей из квадратного
гладкого отверстия, при
n
>2,33а
= 0,077(n-0,5)(n
+ 1).
При п <
2,33 коэффициент
структуры не зависит от степени
ограничения и равняется 0,42;
для струи, выходящей из
круглой гладкой трубы, при п
> 2,33
а = 0,062
(п -0,5)
При n<2,33
а =
0,034;
для плоской струи при п
> 3,12
а = 0,2(n3/2
— 1,25n
+ 0,25)/(n3/2
— 1),
а при п
<3,12
а= 0,085 .
Подобное явление было установлено при
распространении свободных ветровых
струй в карьерах.
По В.Н. Воронину, коэффициент
турбулентной диффузии струи, и начальном
сечении которой уже имеется некоторое
количество газа с содержанием с0
(частично загазованная
струя), определяется по формуле
.
(6.46)
При этом принимается, что
коэффициент турбулентной диффузии не
зависит от турбулентной структуры
газовоздушной среды вне свободной
струи, т.е. газообмен между струей и
средой определяется лишь течением в
струе, а это, видимо, справедливо только
для затопленных струй, распространяющихся
в неограниченном пространстве.
Поскольку через границу свободной струи
происходит обмен турбулентными массами,
то турбулентная структура струи должна
зависеть от структуры движения и энергии
привносимых в нее извне масс. При
исследованиях затопленных струй,
распространяющихся в ограниченных
пространствах (карьеры, тупиковые
выработки, камеры и т.п.), была установлена
завиисимость их угла раскрытия [а
следовательно, в соответствии с формулой
(6.36) —
и коэффициента
структуры струи] от геометрии ограничивающих
поверхностей, что должно быть связано
с турбулентной структурой вторичных
токов, заполняющих пространство между
ограничивающими поверхностями и границей
свободной струи.
В общем случае структура
вторичных токов должна зависеть от
начального расхода воздуха в струе, и
для точного описания газообменных
процессов, связанных с распространением
свободных струй, с помощью коэффициента
турбулентной диффузии В.Н. Воронина
необходимо определить зависимости его
от диффузионных свойств внешней
среды (например, от коэффициента диффузии
Dт.
Более строгим является исследование
процессов газопереноса в свободных
струях на основе ранее рассмотренных
коэффициентов турбулентной диффузии
и эффективных коэффициентов диффузии.
Исследования для их установления
применительно к струйным движениям
в горных условиях были выполнены К.Ю.
Лайгна, Э.А. Поттером и О.А. Суллакатко.
Ими впервые получены выражения для
расчета коэффициентов турбулентной
диффузии ограниченной (степенной)
струи [12], в частности, для эффективных
коэффициентов продольной турбулентной
диффузии:
круглая турбулентная струя,
30·103
< Rе
< 730·103:
плоская турбулентная струя,
30·103
< Rе
< 730·10:
(6.49)
Здесь
— коэффициент стеснения
струи; S
— площадь поперечного сечения выработки;
d
— начальный диаметр струи; и
—средняя начальная
скорость струи; Н, В —
соответственно высота
и ширина выработки; b
— начальная ширина
струи.
Примеры расчета кинематического коэффициента диффузии в газах.
Примеры расчета кинематического коэффициента диффузии в газах. Пример X1 * 2.Двухкомпонентный бензол-коэффициент диффузии диоксида углерода при температуре n1, 2 определяет Г-273,2 + 45 = 318,2 ÅK. Решение. Физико-химические константы(согласно таблице) Молекулярная масса: Бензола M] = 78, и Углекислый газ м3 = 44.01 См. таблицу-фракции для расчета молярного объема при температуре кипения Le Ba. 11-2 (единица измерения cmg {mol) Углерод 14.8. Водород 3.7 6-членное кольцо-15,0 Углекислый газ 34,0 Температура кипения (°к): Бензол T’TSP. 1 = 353,3 Диоксид углерода. 。 7 ′ Кип. 2-194.7 Важные константы: Бензол 7cr. Я-561.7<sup class=»reg»>®</sup>К; РКР= 47.7 Диоксид углерода. 。 ГКЛ-а-304.3°к; РКР> 2 −73.0 в Расчетные физико-химические константы молярного объема (СМГ {моль)) Бензол^ » ute<, = 6.14.8 + b-3.7-15 = 96 Диоксид углерода. 。 ^ МКП <2-34″ ° 0,0748 см2; секунды / 78.11 +44.01 1/2 я 78.11-44.01) 0.0043 (318.2) e / g 1(961 /<sup class=»reg»>®</sup>+ 34 | /<sup class=»reg»>®</sup>)2 2.
Углекислый газ-определяет кинематический коэффициент диффузии воздушной системы при температуре 1300°К и давлении 01.2.
Людмила Фирмаль
- Р = 0,96. 0.00837 (318.2)>; * 1 (96 | / s + 341?3) 2(318,2 + 368) 78.11 + 44.01 В С 78.11-44.01 / = 0,0678 см7; секунды Сі 2 = 1.47.0.96 / ■ 353.3 194.7 = 368 ^ 2.74•10 «<•0.771.823-47.7、/ с * 73,/ С. 561.75 / / а•304.35 / ’ 2ч、 И, 2 1 ^ x (~^ gtag) ’^ 0 ′ 0664 см Согласно таблице система с бензолом Х1 3 кинематическая диффузионная способность-углекислый газ 01.2 = 0.0715 см2 {сек. Ошибка расчета (%): Метод гириланда+4.6 Закон Арнольда— 5.2 Закон слатери-6,9 Пример X1-3. Рассчитайте кинематический коэффициент диффузии системы азот-монооксид углерода 0|, *с использованием методов: Хиршфельдера, Бирда, Кертисса (/=15°С; Р = 1 Ат). Решение.
Постоянная мощность (вкладка Мб ): Азот Е,/ А = 91,5°к, о = 3,681 а
Окись углерода…. Е2 / л = 110.0°к,<Р2 = 3.590 с Молекулярная масса M
| = M2 = 28 110.0 = 1 100.4°к Р = / ^ р = ^ 91.5 ′ ™т 273 + 15 ЛВ»
Ширина 2.87 Таблица интерполяционных данных. UP-1, выглядит так: 1 * =
2.87;и (|-1)* = 0.960 поправочный коэффициент (вкладка UI-2) −1.0035
3,681 + 3,590 0|, 2 2-3, 635 По формуле (Х1-25): 1.0035 = 0.197 см ^ / с
/ 28 + 28 | / а / 2•28•28 / 3. Метод Слаттери: 7 ^ = 318.2. п /
561.7-304.3 = 0,770 0,002628 * 2883 / г и»-«-1(3.635)*•0.960 O * найдено
экспериментально. 2 = 0,211 см2 / С.
- Погрешность расчета составляет-6,6%. Пример X1-4 [13]. Используя методы Хиршфельдера, Бирда и Кертисса, определите коэффициент диффузии системы углекислый газ (I) углекислый газ(2) при температуре f = 80°C и давлении p = 2 атм(используя номограмму Куонга).Решение. Согласно таблице Mb, вы можете увидеть следующее: Углекислый газ (МГ = 44). 。 .a, 3.996 A; E | / 6«190°K окись углерода (M2 = 28)a2 = 3.590 A; e2 / A = 110°K Мощность постоянная смеси: ля., ^ 3’996 + 3’590 = 3.79 A + =(190-110) n-5 = 144°K T = 273.2 + 80 = 353.2 ’K; G * = = 2.45 6), s 144 Согласно таблице, значение t * = 2.45. VII*! (Интерполяция) соответствует 1.007. / Л ^ + МПО. Д._ / 44 + 28 0.5 _ 2 2 I2M1M2)2•44•28.
Выполните действия в номограмме(рисунок XM): 1) соединяя прямые p = 2 и 7 = 353.2, определяем точки на оси P. 2) соединяя точку на оси φ с 0 = 3.79, определяется точка на оси оси. 3) точка на оси P соединяется с P-1.007 для определения точки на оси 5. 4) путем соединения точки на оси 5 с M = 2.42, 0 = 0.103 см2 / с определяется. Согласно формуле (X1-25), 0 = 0,103 см *(сек. Пример X1-5. 20°С и 0 в 1|, g-0,151 см21 секунды(из таблицы X1-5).
Применяют формулу (Х1-39), исключающую отношение поправочного коэффициента. Решение. Постоянная мощность : Воздух= 97,0 Диоксид углерода.
Людмила Фирмаль
- 。 Е2 / −190.0 = В 97190 = 135.4 К = 2,162; = 9.58; = 1.017%1)* = 0.7491 КТХ 293.2 с 2 135.4 KT7 1300 С | г 135.4 VII I в соответствии с таблицей: После соответствующих замен, вы получаете: 0.300-ы = 0.151 ^ = 1.96 cmsec XI * 5 экспериментальное значение согласно таблице 0 | зд> 0-1. 97 см2 / С. Погрешность расчета составляет −0,5%. Пример X1-6.Рассчитайте кинематический коэффициент диффузии системы монооксид углерода-азот при температуре f = 50°C и давлении p = 105 Ат. Система согласно таблице х1-5 коэффициент диффузии монооксида углерода-азота при температуре 15°с, давлении 1 0,211 см2!Считанные секунды. ^ 15°С Решение. Важные константы: Азот ГКР., = 126°к; /?cr. 1 = 33.5 АГ Окись углерода…. ТКР. 2 = 134,5°к; РКР. 2 * 34.6 в Псевдокритические параметры (*1 = A ’ 2™05): Т’ПС. cr.«= 0.5 −126 + 0.5•134.5 = 130.2° К РПС. cr.= 0.5•33.5 4 0.5•34.6 = 34.0 на. Да. О 50°с = О0 15 С №.Г ’ 0.211. ^ 273.2 + 50 * / * 273.2 4 15 ^ 0,236 см21 сек.
Коэффициент диффузии при 50°C и I-уравнение (X1-37): Заданный параметр: 323.2 105. * pr-130 2-Lf 34» ^ 1050.、 1г (зп) П (почтовый°) −0.95. 0,95; ОП = 0.95 * 0.236 105 = 0,00214 см2 (секунда Согласно рисунку рисунка X1-3: Расчетные значения не были проверены экспериментально. Соотношение рассчитанных значений выглядит следующим образом: ИП 0.00214-5°= 0.236 = 0.009 В аналогичном случае (см. Коэффициент самодиффузии азота под давлением 104,5 Ат) Er / d0 * * 0,0084.Поэтому можно предположить, что результат расчета близок к фактическому значению ОП. Пример X1-7.
Смотрите также:
Методические указания по гидравлике
Возможно эти страницы вам будут полезны:
- Полуэмпирические формулы для расчета кинематического коэффициента диффузии.
- Диффузия в многокомпонентных газовых системах.
- Выбор метода расчета кинематического коэффициента диффузии в газах.
- Критерий Шмидта
Коэффициент диффузии жидкостей, газов и паров (Таблица)
Диффузия газов и паров
Коэффициент самодиффузии газов, D’ вычисленный теоретически, оказывается приближенно равным величине: D’= 1/3 u * λ, где u и λ обозначают среднюю скорость и среднюю длину свободного пробега молекул газа. Если же имеет место взаимная диффузия двух различных газов, то коэффициент диффузии D определяется выражением
где D’ и D» — коэффициенты самодиффузии первого и второго газов, а n1 и n2 — их концентрации. Эта формула выведена в предположении, что диффузия происходит стационарно, т. е. что суммарная концентрация обоих газов остается постоянной (n1 + n2 = const). При этом условии коэффициент диффузии первого газа во второй (D12) оказывается равным коэффициенту диффузии второго газа в первый (D$i), т. е. что D12 = Dh=D.
В справочной таблице диффузия предполагается, что газ (пар) диффундирует в пространство, занятое другим газом (паром); оба тела находятся при постоянных условиях давления и температуры (последняя указана в таблице). Значения коэффициентов диффузии К даны в см2/сек.
Газы |
t (°C) |
К |
Аргон — гелий |
15 |
0,70 |
Водород — азот |
12,5 |
0,73 |
— воздух |
0 |
0,64 |
— кислород |
0 |
0,69 |
— кислород |
14 |
0,77 |
— метан |
0 |
0,62 |
— окись углерода |
0 |
0,64 |
— углекислота |
0 |
0,54 |
— этилен |
0 |
0,48 |
Воздух — водяной пар |
0 |
0,22 |
— водяной пар |
15 |
0,26 |
— кислород |
0 |
0,18 |
— сероуглерод |
0 |
0,09 |
Кислород — азот |
0 |
0,18 |
— азот |
12,5 |
0,20 |
Окись углерода — этилен |
0 |
0,11 |
Углекислота — водород |
0 |
0,54 |
— водород |
18 |
0,60 |
— водяной пар |
18 |
0,15 |
Бензол — водород |
0 |
0,29 |
— воздух |
0 |
0,14 |
— кислород |
0 |
0,18 |
— метан |
0 |
0,15 |
— окись углерода |
0 |
0,13 |
Водяной пар — водород |
0 |
0,69 |
— углекислота |
0 |
0,13 |
Кислота уксусная — водород |
0 |
0,40 |
— воздух |
0 |
0,11 |
— углекислота |
0 |
0,07 |
Сероуглерод — водород |
0 |
0,36 |
— воздух |
0 |
0,09 |
— углекислота |
0 |
0,06 |
Спирт амиловый — водород |
0 |
0,23 |
— воздух |
0 |
0,06 |
— углекислота |
0 |
0,04 |
Спирт метиловый — водород |
0 |
0,50 |
— воздух |
0 |
0,13 |
— кислород |
0 |
0,18 |
— метан |
0 |
0,15 |
— окись углерода |
0 |
0,13 |
— углекислота |
0 |
0,08 |
Спирт этиловый — водород |
0 |
0,37 |
— воздух |
0 |
0,10 |
— углекислота |
0 |
0,07 |
Диффузия водных растворов, жидкостей
Раствор диффундирует в чистую воду. Значения коэффициентов диффузии выражены в см2/сутки. Концентрация выражена для электролитов в грамм-эквивалентах на литр раствора, для неэлектролитов — в молях. Температура t °С указана в таблице.
Вещество |
t (°С) |
Конц. |
D |
Аммиак |
12 |
1.0 |
1,42 |
4 |
3,55 |
1,06 |
|
3,75 |
0,45 |
||
Барий хлористый |
8 |
0,33 |
0,65 |
2.4 |
0,66 |
||
Бром |
12 |
0,1 |
0,8 |
Глицерин |
10,14 |
0,125 |
0,356 |
0,875 |
0,342 |
||
1,75 |
0,300 |
||
Кадмий сернокислый |
19,0 |
2,0 |
0,246 |
Калий азотнокислый |
17,6 |
0,02 |
1,28 |
0,3 |
1,26 |
||
1,4 |
1,10 |
||
Калий бромистый |
10 |
1,0 |
1,13 |
Калий сернокислый |
19,6 |
0,02 |
1,01 |
0,28 |
0,86 |
||
0,95 |
0,79 |
||
Калий углекислый |
10 |
3,0 |
0,60 |
Калий хлористый |
17,5 |
0,02 |
1,36 |
0,9 |
1,52 |
||
18,0 |
1,0 |
1,330 |
|
3,6 |
1,338 |
||
Кальций хлористый |
9 |
1.5 |
0,72 |
10 |
0,27 |
0,68 |
|
2,0 |
0,68 |
||
Кислота азотная |
12 |
0,55 |
1,91 |
19,5 |
0,10 |
2,07 |
|
0,90 |
2,26 |
||
3,90 |
2,46 |
||
Кислота серная |
12 |
1,0 |
1,12 |
18 |
0,35 |
1,32 |
|
0,85 |
1,34 |
||
2,85 |
1,60 |
||
4,85 |
1,90 |
||
Кислота соляная |
12 |
0,50 |
2,07 |
1,0 |
2,09 |
||
19,2 |
0,10 |
2,21 |
|
0,90 |
2,63 |
||
3,20 |
3,89 |
||
Кобальт хлористый |
10 |
1.5 |
0,46 |
Литий бромистый |
10 |
2,3 |
0,80 |
4,4 |
0,90 |
||
Литий хлористый |
18 |
0,10 |
0,951 |
1,0 |
0,920 |
||
4,2 |
0,956 |
||
Медь сернокислая |
17 |
0,10 |
0,39 |
0,50 |
0,29 |
||
1,95 |
0,23 |
||
Медь хлористая |
10 |
1.5 |
0,43 |
Натрий азотнокислый |
10,5 |
3,0 |
0,76 |
5,0 |
5,0 |
0,83 |
|
13 |
0,6 |
0,90 |
|
6,0 |
0,77 |
||
Натрий бромистый |
10 |
2,9 |
0,86 |
Натрий сернокислый |
10 |
1,4 |
0,66 |
10,4 |
1,29 |
0,49 |
|
Натрий углекислый |
10 |
2,4 |
0,39 |
Натрий хлористый |
12 |
1,0 |
0,90 |
15 |
0,02 |
0,94 |
|
0,1 |
0,94 |
||
0,9 |
0,97 |
||
3,9 |
1,02 |
||
Сахар тростниковый |
12 |
0,5 |
0,28 |
1,0 |
0,25 |
||
18,5 |
0,30 |
0,31 |
|
0,97 |
0,24 |
||
1,97 |
0,43 |
||
Свинец азотнокислый |
12 |
0,22 |
0,71 |
0,82 |
0,66 |
||
Серебро азотнокислое |
12 |
0,02 |
1,03 |
0,10 |
0,98 |
||
0,90 |
0,88 |
||
3,9 |
0,53 |
||
Спирт этиловый |
11 |
0,05 |
0,73 |
0,25 |
0,69 |
||
0,75 |
0,62 |
||
3,75 |
0,45 |
||
Цинк сернокислый |
19,5 |
0,025 |
0,50 |
0,55 |
0,36 |
||
2,95 |
0,33 |
_______________
Источник информации: КРАТКИЙ ФИЗИКО-ТЕХНИЧЕСКИЙ СПРАВОЧНИК/ Том 1, — М.: 1960.
Поделитесь ссылкой с друзьями: