Постоянная скорости диффузии единицы измерения. Коэффициенты диффузии
Лето! Солнце, тепло, в воздухе запах свежескошенной травы. Трава вон как далеко отсюда, а запах чувствуется. Это уже физика, подобное явление называется диффузией. И определяется оно как проникновение частиц соприкасающихся веществ друг в друга. Например, как сейчас, трава соприкасается с воздухом, и запах от скошенной травы распространяется по воздуху далеко от места косьбы. А характеризуется подобный процесс обычно такой величиной, как коэффициент диффузии.
Можно сказать, что всё происходящее вокруг нас во многом связано с процессами взаимопроникновения. Благодаря такому явлению мы даже живем. А что тут удивительного? Проникновение кислорода или питательных веществ в кровь - это и есть самая настоящая диффузия, одно вещество проникает в другое. Распространена диффузия в природе гораздо шире, чем нам кажется. Подобное явление не является чем-то экзотическим, а наоборот, достаточно широко представлено в окружающем мире. Даже то, что выдыхаемый человеком воздух не скапливается вокруг, а рассеивается в пространстве, тоже обусловлено диффузией.
Происходить взаимопроникновение может в самых разных телах - жидких, газообразных, твердых. Его причиной служит хаотическое движение молекул вещества. Кстати, такое явление, как диффузия считается в молекулярно-кинетической теории одним из факторов, подтверждающих ее. Так вот, возвращаясь к агрегатному состоянию тела: скорость проникновения вещества зависит от его и от такой величины, как коэффициент диффузии.
Что же это за таинственное понятие? Так называется количественная характеристика скорости переноса молекул одного вещества в какое-либо другое. Коэффициент диффузии, формула которого достаточно специфична, позволяет провести оценку количества вещества, прошедшего через единицу площади (квадратный метр) за единицу времени (секунду).
На практике установлено, что в газах взаимное проникновение происходит с наибольшей скоростью, а в твердых телах скорость проникновения минимальна. На коэффициент диффузии оказывают влияние температура тела и взаимная концентрация находящихся рядом веществ. С ростом температуры скорость взаимопроникновения увеличивается, с ростом концентрации вещества - тоже.
Таким образом, явление взаимопроникновения обусловлено так называемым градиентом концентрации или градиентом температуры. По типу диффузия разделяется на свободную и вынужденную. Вынужденная возникает под воздействием внешних сил. В зависимости от их вида вынужденную диффузию определяют как термо-, электро-, баро- и восходящую.
Диффузия широко используется в технике. Один из самых типичных примеров - диффузионная сварка. Суть подобной технологии проста: соединяют два различных тела (пусть это будут два металла), а затем подвергают их воздействию давления и температуры. Последняя по своей величине меньше, чем температура плавления веществ.
В итоге получается соединение двух разнородных материалов. Такая технология широко применяется в приборостроении и при изготовлении крупногабаритных деталей сложной формы, в опытном и мелкосерийном производстве. Подобная сварка может проводиться в разных условиях, вплоть до условий вакуума, все определяется конкретными требованиями, предъявляемыми к готовому изделию.
Другим, не менее распространенным, использованием явления диффузии является его применение для получения полупроводниковых структур. Одна из технологий создания р-n переходов основана на явлении взаимопроникновения. Под воздействием высоких температур, приближающихся к в кристаллической структуре получают области с нужной концентрацией примесей.
Рассмотрение понятия «коэффициент диффузии» позволило убедиться в чрезвычайно широком распространении явления взаимопроникновения веществ в природе, а также в разнообразных вариантах его использования в технике.
Как уже упоминалось, кинетическая теория газов позволяет просто объяснить факт медленности процесса диффузии, несмотря на большие значения скоростей тепловых движений молекул. Это обусловлено тем, что молекулы газа, чтобы попасть из одной точки в другую, вследствие столкновений вынуждены пройти длинный зигзагообразный путь, во много раз превосходящий расстояние по прямой между этими точками.
Кроме такого качественного объяснения, кинетическая теория позволяет и количественно оценить величину коэффициента диффузии и выразить его через молекулярные величины - длину свободного пробега молекул и скорости их тепловых движений.
Рассмотрим площадку 5 в сосуде с газовой смесью, перпендикулярную к оси X (рис. 51), вдоль которой поддерживается
постоянная разность концентраций (речь идет, следовательно, о стационарном процессе). Примем для определенности, что Из-за тепловых движений молекулы интересующего нас компонента будут переходить через площадку 5 как слева направо, так и справа налево. Ввиду существующей разности концентраций по обе стороны площадки возникнет некоторый диффузионный поток вдоль оси X, равный, очевидно, разности между числом молекул пересекающих площадки 5 в 1 с (перпендикулярно к ее плоскости) в направлении положительных значений X (вправо), и числом молекул пересекающих то же сечение и за то же время в противоположном направлении (влево):
Как определить число молекул, пересекающих площадки? Если бы все молекулы двигались с одинаковой скоростью направленной по оси X, то число молекул, переходящих в 1 с площадку в было бы равно где число молекул в единице объема.
В действительности существует распределение молекул по скоростям, но для грубой оценки мы примем, что у всех молекул одна и та же скорость, равная средней скорости Примем также, что тепловые скорости молекул равномерно распределены по трем взаимно перпендикулярным направлениям. Тогда из всех молекул единицы объема движется вдоль оси X и из них половина движется в положительном направлении оси X, т. е. по направлению к площадке, в то время как другая половина движется в противоположном направлении - от нее.
Следовательно, число молекул пересекающих площадки в 1 с слева направо, и число молекул пересекающих ту же площадку в противоположном направлении, выразятся соотношениями:
Здесь концентрации молекул с одной и с другой стороны от площадки. Относительно значений величин необходимо заметить, что они изменяются вдоль оси X вследствие столкновений молекул между собой. Поэтому к выделенной нами площадке молекулы подходят, имея те значения концентраций которые создались при последнем столкновении перед площадкой. Значит, мы должны считать, что это те числа молекул в единице объема, которые были на расстоянии (средняя длина свободного пробега) от площадки, по обе стороны от нее.
Диффузионный поток следовательно, определяется выражением
где разность концентраций между точками, отделенными друг от друга расстоянием в Разность эту нетрудно определить, если известно значение градиента концентрации - (будем полагать, что изменяется только по оси т. е. вместо можно писать Так как есть разность концентраций, приходящаяся на единицу длины, то на расстоянии она равна:
Эта формула справедлива, если X достаточно мало.
Таким образом, для диффузионного потока получаем выражение:
или, умножив обе части этого равенства на массу молекулы имеем:
Сравнивая его с уравнениями (40.2) и (40,3) закона Фика
находим интересующее нас выражение для коэффициента диффузии:
Из этого выражения видно, что коэффициент диффузии обратно пропорционален давлению газа (потому что и прямо пропорционален квадратному корню из температуры
При выводе формулы (42.1) не принималась во внимание диффузия второго компонента, крторая, разумеется, тоже происходит. И она не может не влиять на диффузию рассматриваемого компонента смеси.
Так, например, в смеси водорода и углекислого газа водород должен диффундировать значительно быстрее, чем углекислый газ, потому что при данной температуре средняя скорость тепловых движений молекул водорода почти в пять раз больше и, кроме того, длина свободного пробега молекул водорода тоже больше. Но это
значит, что объем водорода, переносимого в одном направлении, больше объема углекислого газа, переносимого в противоположном направлении. В таком случае в газе неизбежно возникает разность давлений, а значит, и поток газа в целом. Между тем мы определили диффузию как процесс, при котором давление газовой смеси во всех точках остается постоянным и газ, как целое, покоится. На самом же деле в процессе взаимной диффузии двух различных газов, одного в другой, превышение диффузионного потока одного из них над потоком другого уравновешивается течением всего газа по направлению к той области, где первоначально находились более быстро диффундирующие молекулы.
Этого обстоятельства мы, однако, не принимали во внимание при выводе формулы (42.1) для коэффициента диффузии, и эта формула справедлива в сущности только для диффузии молекул газа в среде того же газа. Такой процесс называется самодиффузией, а формула (42.1) выражает, следовательно, коэффициент самодиффузии.
С таким явлением мы имеем дело, например, когда газовая смесь состоит из двух различных изотопов одного и того же вещества, лишь незначительно отличающихся друг от друга своей массой, но не отличающихся никакими другими свойствами. Если один из изотопов радиоактивен, то такую самодиффузию легко наблюдать, так как за проникновением радиоактивных частиц можно следить по их излучению.
Коэффициенты диффузии
В выражениях для диффузионных газовых потоков коэффициенты молекулярной и турбулентной диффузии являются единственными параметрами, учитывающими свойства среды. Естественно, что эти величины имеют сложный характер, и их определение - одна из важных задач теории диффузионных процессов.
Коэффициенты молекулярной диффузии . Для газов со сходными молекулами (имеющими почти равные массы и эффективные сечения) Максвелл получил следующее выражение для коэффициента молекулярной диффузии:
где - длина свободного пробега молекул; 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.11) векторы j т и gradс коллинеарны*. Следовательно, согласно определению, направление вектора gradс является главным направлением тензора, а ось координат, соответствующая ему, - главной осью. Нахождение главных осей тензора коэффициентов диффузии для выработки - в ряде случаев задача неопределенная, так как для этого необходимо знать поверхности равных содержаний в потоке, т.е. поле содержаний, что обычно является конечной задачей исследований. Лишь в простых случаях диффузии главные направления могут быть определены достаточно просто. Например, при газовыделении с одной стенки gradс с некоторым приближением можно принять нормальным к этому боку и, следовательно, главные оси тензора будут направлены вдоль потока воздух и перпендикулярно к нему. В более сложных случаях главные оси тензора могут иметь и другие направления.
Следует отметить, что принятие тензора D т y симметричным для случая движения воздуха в горной выработке является также определенным допущением. Для неоднородного и неизотропного турбулентного потока, каким является вентиляционный поток в выработке, тензор коэффициентов турбулентной диффузии будет несимметричным. Ниже отмечается, что компоненты тензора коэффициентов турбулентной диффузии могут быть выражены через усредненное произведение (корреляцию) мгновенных значений пульсационной скорости и ni и пути перемешивания для содержания (здесь i ,j = х, у, z, и ni = и п ; u пу = v n u nу =w n). Для симметричного тензора должны соблюдаться равенства , что приводит к соотношениям . Однако для неизотропных вентиляционных потоков корреляция несимметрична относительно i и j , а это не отвечает приведенным равенствам. Несимметричность тензора коэффициентов турбулентной диффузии для шахтных вентиляционных потоков косвенно доказывается фактором различной интенсивности турбулентной диффузии в разных направлениях.
Отмеченные приближения, которые применяют при решении практических задач шахтной газовой динамики, в настоящее время не имеют оценки. Применительно к условиям диффузии в приземном слое атмосферы погрешности незначительны (в некоторых случаях они составляют 15-20 %). Однако степень анизотропности шахтных вентиляционных потоков значительно выше атмосферных, что может привести к необходимости учета факта несимметричности тензора диффузии.
Второй способопределения коэффициента турбулентной диффузии основан на использовании теории Прандтля о пути перемешивания, согласно которой компоненты потока газа можно определять как сумму трех слагаемых:
. (6.21)
Здесь, подобно тому, как это было принято в выражении (6.19) - , суммирование производится по двойному индексу (j ); i = j = х,у,z; ; L c - путь перемешивания для содержания.
Из выражения следует, что коэффициент турбулентной диффузии является тензором второго ранга
(6.22)
определяемым девятью компонентами -
Сопоставляя методы выражения коэффициента турбулентной диффузии по Буссинеску и Прандтлю, видим, что в первом случае коэффициент турбулентной диффузии остается неопределенным, во втором - определяется через характеристики турбулентного движения ().
В случае плоского потока () коэффициент турбулентной диффузии в поперечном к основному движению направлении определяется из выражения (6.21):
В случае изотропной турбулентности можно принять L cx = L су, что приводит к равенству
т.е. в этом частном случае коэффициент турбулентной диффузии является скаляром.
Если в уравнении (6.23) v n выразить по Прандтлю через путь перемешивания для импульса L, то для плоского потока получим выражение
, (6.24)
где среднее квадратичное значение v п
а 1 - коэффициент пропорциональности между u n и v п. Если принять, что
L / L с = а 2 = сопst, (6.25)
. (6.26)
Величина для случая диффузии газа является аналогом пути перемешивания для импульса по Карману (не тождественному прандтлевскому пути перемешивания).
Из уравнения видно, что, имея какие-либо гипотезы относительно величин l С (), можно, измеряя в потоке, определить коэффициент турбулентной диффузии. Наиболее простым допущением является отождествление l с путем перемешивания для импульса l ; во многих случаях такое приближение дает вполне удовлетворительные результаты.
Следующим шагом в этом направлении является принятие пропорциональности между l с и l ; значение коэффициента пропорциональности между ними зависит от свойств диффундирующего газа, разности содержаний газа в диффундирующем объеме и в среде. По имеющимся сведениям, этот коэффициент больше 1; для азота он равен ~ , для гелия ~ . Имеются попытки оценить l с через l и критерий Ричардсона, характеризующий затухание турбулентности под действием объемных (гравитационных) сил при диффузии активного газа.
Наконец, третий способ определения коэффициента турбулентной диффузии основан на представлении процесса диффузии как случайного движения жидких частиц, первоначально сконцентрированных в некоторой области. Бэтчелор показал, что и в этом случае коэффициент турбулентной диффузии является тензором второго ранга. Запись его (для случая однородной турбулентности), однако, имеет иной вид:
,
где у i , y j - лагранжевы координаты жидкой частицы, величины случайные, являющиеся функцией времени.
Представление коэффициента турбулентной диффузии в виде тензора имеет в основном теоретическое значение. В настоящее время практически ничего неизвестно о недиагональных компонентах этого тензора. Изученные в какой-либо степени компоненты тензора диффузии - это диагональные компоненты D т xx ,D т yу , D т zz , которые в дальнейшем и будут рассматриваться. Для простоты написания обозначим D т xx = D тх и т.д.
Необходимо отметить, что в общем случае коэффициент турбулентной диффузии является функцией координат. Это можно видеть, например, из уравнения (6.24), где величины dи/dу, v ´ n , l с для потоков в горных выработках являются функциями поперечных координат , а в некоторых случаях (изменение сечения по длине выработки, свободные струи) - и продольной координаты. Эти же величины являются и функциями скорости потока (точнее, числа Рейнольдса - Rе* потока), что говорит о существовании зависимости коэффициента турбулентной диффузии и от числа Rе.
Данные о коэффициентах турбулентной диффузии в горных выработках немногочисленны, что в значительной степени объясняется техническими трудностями их измерений. Имеющиеся сведения частично основываются на данных о коэффициенте турбулентного обмена для импульса и предположении о пропорциональности ему коэффициента турбулентной диффузии .
К.М. Тумаковой были установлены автомодельность поперечных составляющих относительного коэффициента турбулентной диффузии:
; - средняя скорость потока; α – коэффициент аэродинамического сопротивления; r - плотность потока; Н - высота выработки) по числу Рейнольдса, начиная от 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·10 4 до 3,72·10 4 относительные значения последних для штрекообразной выработки прямоугольного сечения, закрепленной рамной крепью из круглого леса с продольным калибром 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 , плотности воздуха r = = 1,22 кг/м 3 , диаметре выработки D = 2,5 м.
Рис. 6.1. Зависимость от у* = = у/Н (Н - высота выработки)
Рис. 6.2. Зависимость от z* = = z/В (В - ширина выработки)
Значения компоненты D т y ·10 3 (м 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)
Для диффузии углекислого газа в воздухе
где k : = 3,96·10 -4 м.
В формулах (6.29)-(6.33) использованы следующие обозначения:
Н- высота выработки, м;
Динамическая скорость, м/с;
u ср - средняя скорость воздушной струи, м/с;
α - коэффициент аэродинамического сопротивления, Н·с 2 /м 4 ;
r -плотность воздуха, кг/м 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 ; r - плотность воздуха, кг/м 3 .
К.Ю. Лайгна и Э.А. Поттер в своих последних работах* дают следующее выражение для среднего по поперечному сечению эффективного коэффициента диффузии:
, (6.39)
где , а число Рейнольдса рассчитывают по . Значение D э можно определить также по графикам, представленным на рис. 6.3.
Рис. 6.3. Графики к определению эффективного коэффициента турбулентной диффузии
Расчеты по приведенным формулам и графикам дают значения D э, порядка нескольких м 2 /с.
Дж. Тэйлор и К.Ю. Лайгна, исследуя влияние изогнутости канала на коэффициент диффузии, сделали вывод, что этот фактор может увеличивать D э до двух раз.
Для изучения газодинамических процессов при действии свободных струй В.Н. Воронин применил коэффициент турбулентной диффузии k т, определив его как отношение среднего содержания газа в поперечном сечении ядра постоянной массы свободной струи с я к среднему содержанию на ее границе с гр:
Значения 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/b 0 для плоской (б) струй (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(n 3/2 - 1,25n + 0,25)/(n 3/2 - 1),
а при п <3,12
а = 0,085 .
Подобное явление было установлено при распространении свободных ветровых струй в карьерах.
По В.Н. Воронину, коэффициент турбулентной диффузии струи, и начальном сечении которой уже имеется некоторое количество газа с содержанием с 0 (частично загазованная струя), определяется по формуле
. (6.46)
При этом принимается, что коэффициент турбулентной диффузии не зависит от турбулентной структуры газовоздушной среды вне свободной струи, т.е. газообмен между струей и средой определяется лишь течением в струе, а это, видимо, справедливо только для затопленных струй, распространяющихся в неограниченном пространстве. Поскольку через границу свободной струи происходит обмен турбулентными массами, то турбулентная структура струи должна зависеть от структуры движения и энергии привносимых в нее извне масс. При исследованиях затопленных струй, распространяющихся в ограниченных пространствах (карьеры, тупиковые выработки, камеры и т.п.), была установлена завиисимость их угла раскрытия [а следовательно, в соответствии с формулой (6.36) - и коэффициента структуры струи] от геометрии ограничивающих поверхностей, что должно быть связано с турбулентной структурой вторичных токов, заполняющих пространство между ограничивающими поверхностями и границей свободной струи*.
В общем случае структура вторичных токов должна зависеть от начального расхода воздуха в струе, и для точного описания газообменных процессов, связанных с распространением свободных струй, с помощью коэффициента турбулентной диффузии В.Н. Воронина необходимо определить зависимости его от диффузионных свойств внешней среды (например, от коэффициента диффузии D т.
Более строгим является исследование процессов газопереноса в свободных струях на основе ранее рассмотренных коэффициентов турбулентной диффузии и эффективных коэффициентов диффузии. Исследования для их установления применительно к струйным движениям в горных условиях были выполнены К.Ю. Лайгна, Э.А. Поттером и О.А. Суллакатко. Ими впервые получены выражения для расчета коэффициентов турбулентной диффузии ограниченной (степенной) струи , в частности, для эффективных коэффициентов продольной турбулентной диффузии:
круглая турбулентная струя, 30·10 3 < Rе < 730·10 3:
плоская турбулентная струя, 30·10 3 < Rе < 730·10:
(6.49)
Здесь - коэффициент стеснения струи; S - площадь поперечного сечения выработки; d - начальный диаметр струи; и - средняя начальная скорость струи; Н, В - соответственно высота и ширина выработки; b - начальная ширина струи.
ТЕМА №7. ИНТЕГРАЛЬНЫЕ ГАЗОДИНАМИЧЕСКИЕ ЭФФЕКТЫ В ШАХТАХ
Общие положения
Во многих случаях полезные результаты могут быть получены более простым интегральным методом.
Интегральный метод, или метод усредненных характеристик, основан на том факте, что локальные газодинамические эффекты в конечном итоге выступают в виде некоторых обобщенных, или интегральных, усредненных закономерностей, таких, например, как изменение среднего по сечению содержания газа на выходе из забоя, участка, шахты, изменение суммарного дебита газа из выработанного пространства в выработку и др. Предметом интегрального метода анализа является исследование газодинамических процессов в шахтах, описанных в терминах усредненных характеристик. В интегральном методе использованы такие усредненные характеристики процесса, как средняя скорость движения, среднее содержание и т.п. Усреднение может производиться по одному, двум, трем измерениям, а также во времени*. Использование усредненных характеристик не требует знания их полей, что существенно упрощает аппарат анализа, а также турбулентных характеристик потока, которые в этом методе обычно учитываются эмпирическими константами. Число параметров процесса сокращается. Все это делает интегральный метод достаточно простым и легко приводящим к конкретным результатам.
Однако замена локальных значений характеристик их усредненными значениями не всегда может пройти без заметного снижения достоверности получаемых результатов. В качестве примера можно привести расчет расхода воздуха по среднему содержанию метана для выработок, где имеются слоевые скопления газа: при достаточности полученного расхода в среднем он не всегда может обеспечить ликвидацию зон высокого содержания газа в пристеночных областях. В общем, правомерность и погрешность использования метода средних характеристик определяется их различиями, как в точке, так и во всей области движения, причем, чем больше области, где эти различия существенны, тем менее правомерно использование этого метода и тем больше возникающая при этом погрешность.
Основным соотношением интегрального метода является соотношение между содержанием газа с , объемным расходом газа J и расходом газовоздушной смеси Q* :
Если J и Q взаимонезависимы, изменение содержания прямо пропорционально расходу газа и обратно пропорционально расходу воздуха. В частности, монотонному изменению расхода воздуха в этом случае соответствует монотонное изменение содержания газа.
В ряде случаев расход газа в выработке зависит от расхода воздуха. При этом возможно появление так называемых переходных газодинамических процессов, при которых обратно пропорциональная зависимость между с и Q нарушается. В общем имеющий практическое значение характер зависимости с(Q) определяется соотношением (7.1) и зависимостью с (J ).
Основой метода усредненных характеристик является закон сохранения массы, который применяется к участку выработки конечной длины или к выработке в целом. Для выбранного участка определяются интегральные газовые потоки, поступающие в выработку и выходящие из нее, от всех действующих в выработке источников газовыделения. Алгебраическая сумма поступления газа в выработку и его выноса в виде этих потоков за некоторый промежуток времени определяет изменение газосодержания в объеме рассматриваемой выработки за этот же период.
Под интегральным газовым потоком от i -го источника понимается количество газа, поступающее в рассматриваемую выработку в единицу времени. Если i -й интегральный газовый поток в выработку объема V обозначить через J i то, согласно закону сохранения массы и в соответствии с вышесказанным, получим следующее наиболее общее дифференциальное уравнение переноса газа:
(7.3)
где п - число интегральных газовых потоков в выработке; ее - изменение среднего содержания газа в выработке за период времени
В стационарном случае
Интегральные газовые потоки могут поступать в выработку (выходить из нее) либо со струей воздуха, либо с ее твердых границ.
Интегральный газовый поток, вносимый (выносимый) в выработку вентиляционной струей,
(7.5)
где с ср - среднее содержание газа в поступающем (выходящем) в выработку воздухе; Q - расход воздуха на входе (выходе) в выработку.
Интегральные газовые потоки с твердых границ выработки могут иметь различное происхождение и рассчитываться разными способами. Так, интегральный газовый поток с обнаженной поверхности горных пород
(7.6)
где q г - абсолютное газовыделение с единицы обнаженной поверхности; S - площадь обнаженной поверхности.
Аналогично определяется интегральный газовый поток из граничащего с выработкой выработанного пространства.
Интегральный газовый поток из находящейся в выработке отбитой горной массы может быть определен как произведение абсолютного газовыделения q" г , отнесенного к единице массы отбитой горной породы, на общую отбитую массу М от:
(7.7)
Следует иметь в виду, что q г и q" г являются функциями времени.
Лабораторная работа №7
ИЗМЕРЕНИЕ КОЭФФИЦИЕНТА ДИФФУЗИИ ГАЗА
Диффузией называется явление проникновения двух или нескольких соприкасающихся веществ друг в друга. Процесс диффузии возникает в газе (так же, как и в любом другом веществе), если газ неоднороден по составу, т.е. если он состоит из двух или нескольких различных компонентов, концентрация которых изменяется от точки к точке. Процесс диффузии заключается в том, что каждый из компонентов смеси переходит из тех частей объема газа, где его концентрация больше, туда, где она меньше, т.е. в направлении уменьшения концентрации.
Всякий процесс, при котором параметры системы, участвующей в нем, с течением времени изменяются, называется нестационарным процессом, в отличие отстационарного процесса, при котором величины, характеризующие систему, не изменяются со временем. Диффузия, приводящая к выравниванию концентрации, т.е. к изменению разностей концентраций и самих концентраций компонентов, называется нестационарной диффузией. Можно себе представить и стационарную диффузию, когда тем или иным искусственным путем разность концентраций смеси поддерживается неизменной. Для этого нужно, например, в одну часть сосуда непрерывно добавлять данный компонент, а из другой части сосуда отбирать его в таком же количестве.
Целью настоящей работы является ознакомление с явлением диффузии газов, измерение коэффициента диффузии паров этилового спирта в воздухе при атмосферном давлении, а также ознакомление с экспериментальной методикой измерения коэффициента диффузии паров жидкостей.
1. Основы метода измерения коэффициента диффузии паров жидкостей
Согласно основному закону диффузии (закону Фика) плотность диффузионного потока какого-либо компонента (сорта молекул) пропорциональна градиенту концентрации этого компонента, взятому с обратным знаком:
Ix = − Ddn . | |
Смысл коэффициента диффузии D состоит в том, что он численно равен плотности диффузионного потока, т.е. количеству диффундирующего компонента, проходящего в единицу времени через единицу площади, перпендикулярной к направлению диффузии, при градиенте концентрации, равном единице. Для модели идеального газа коэффициент диффузии равен
π mkT | ||||||||||||
3 πnm πσ2 | ||||||||||||
При фиксированной температуре v | является постоянной, а | Следовательно, |
||||||||||
при постоянной температуре D ~ P 1 . С другой стороны, при фиксированном давлении
l ~ T , а v ~T . Следовательно, при постоянном давлении D ~ T3 / 2 . Эти заключения были тщательно проверены в экспериментах. СоотношениеDP = const соблюдается в довольно широком интервале давлений для не очень плотных газов с точностью до нескольких процентов. В зависимости от температурыD растет несколько быстрее, чем
T3 / 2 , что объясняется дополнительным уменьшениемσ с увеличениемT , приводящим к дополнительному увеличению l .
Коэффициент диффузии для многих газов в воздухе при нормальных условиях имеет
порядок D ~ 10 -5 м 2 /с , как это следует из (1.2), если учесть, что для них l ~ 10− 8 м,v ~ 500 м/с.
Если плотность потока диффундирующего компонента выражать его массой кг , протекающей через1 м 2 площади в единицу временис , то парциальная концентрация компонента в этом случае будет представлять собой парциальную плотность компонентаρ = nm . Уравнение (1.1) в этом случае примет вид
I x = −D | |||
При стационарной диффузии градиент концентраций | остается постоянным |
(неизменным во времени), поэтому постоянен и диффузионный поток. При нестационарной диффузии градиент концентраций изменяется (концентрации выравниваются). Соответственно этому изменяется со временем и диффузионный поток.
Пусть в сосуде с газом (воздухом) находится капля жидкости радиусом R , плотность паров этой жидкости вдали от капли (на бесконечности) -ρ′ ∞ ′ , на поверхности капли -ρ′ 0 ′ ,
причем ρ′ 0 ′ > ρ′ ∞ ′ (капля испаряется). Тогда убыль массы капли (M ) в единицу времени в соответствии с (1.3) может быть представлена в сферической системе координат в виде
− dM = −D | d ρ′′ | 4π r2 . | ||||||
После интегрирования этого выражения получим | ||||||||
− dM =4 πDR(ρ′0 ′ | −ρ′∞ ′) . | |||||||
Если ρ ′′ = 0 , что может быть достигнуто применением поглотителя паров вдали от |
||||||||
капли, то (1.5) можно привести к виду | ||||||||
4 π DR ρ ′′. | ||||||||
Из уравнения состояния идеального газа | P0 μ | |||||||
ρ′0 ′ = | ||||||||
где P 0 - давление паров жидкости на поверхности капли (давление насыщенного пара при температуреT 0 );T 0 - температура поверхности капли;μ - молярная масса жидкости;R 0 -
универсальная газовая постоянная.
Подставляя значение ρ 0 ′′ | из (1.7) в (1.6), получим | ||||||||||||||||||||||||||||||||||||||||||||||
D = − | R0 T0 | ||||||||||||||||||||||||||||||||||||||||||||||
4 πμP | |||||||||||||||||||||||||||||||||||||||||||||||
π R3 ρ , гдеρ - плотность жидкости, поэтому | 4 π R 2 ρ | ||||||||||||||||||||||||||||||||||||||||||||||
T ρ dR | |||||||||||||||||||||||||||||||||||||||||||||||
D = − | |||||||||||||||||||||||||||||||||||||||||||||||
2 Pμ | |||||||||||||||||||||||||||||||||||||||||||||||
Учтем, что | |||||||||||||||||||||||||||||||||||||||||||||||
4 πρ | 4 πρ | ||||||||||||||||||||||||||||||||||||||||||||||
3 2 / 3 | T ρ 1/ 3d (M ) 2/ 3 | ||||||||||||||||||||||||||||||||||||||||||||||
D = − | |||||||||||||||||||||||||||||||||||||||||||||||
Это и есть основное выражение для определения коэффициента диффузии в данной работе.
2. Описание экспериментальной установки
В данной работе использован относительный весовой метод измерения массы испаряющейся капли жидкости (см. рисунок).
Принципиальная схема экспериментальной установки:
1 - капля; 2 - подвес; 3 - пружина; 4 - стеклянный колпак; 5 - конденсатор; 6 - коромысло; 7 - стойка; 8 - поглотитель паров жидкости; 9 - тарелка; 10 - задающий LC-генератор; 11 - частотомер
Суть относительного весового метода, применяемого к настоящей работе, следующая. Капля 1 помещается шприцем на подвес 2, который крепится к плечу коромысла 6, припаянного в средней части к пружине 3, изготовленной из упругого материала (фосфористая бронза). Второе плечо коромысла заканчивается круглой пластинкой, служащей обкладкой конденсатора 5.
При изменении массы капли изменяется зазор между обкладками конденсатора. Это вызывает изменение емкости, которая, как известно, зависит от величины зазора. Конденсатор является составной частью задающего частоту LC - генератора 10, поэтому изменение его емкости влияет на генерируемую частоту. Это изменение частоты регистрируется частотомером. При этом множитель dM3/2 /dt в (1.10) будет пропорционален изменению генерируемой частоты: dM3/2 /dt df/dt.
Весы крепятся на стойке 7 и вся система помещается под стеклянный колпак 4 на тарелке 9. Для поддержания постоянного градиента концентрации при испарении капли служит поглотитель паров жидкости (силикагель) 8.
3. Методика проведения эксперимента
П о д г о т о в к а к о п ы т у
Ознакомиться с описанием и лабораторной установкой. Включить частотомер, дать ему прогреться в течение 15-20 мин.
Внимание! Питание частотомера осуществляется напряжением 220 В, соблюдайте осторожность при работе!
З а д а н и е
В настоящей работе необходимо измерить величину коэффициента диффузии паров этилового спирта (C2 H5 OH) в воздухе.
1. Аккуратно, не задевая весов, снять стеклянный колпак с тарелки. Навешать шприцем каплю этилового спирта, не дотрагиваясь иглой до подвеса. При навешивании
капли подставлять на тарелку под навес кювету для случайно оборвавшихся при навешивании капель.
2. Убрать кювету. Поставить колпак на тарелку.
3. В "ручном" режиме работы частотомера (так удобнее) снять зависимость изменения частоты f от времениt , т.е.f=f(t) . Измерения проводить через 20 с, контролируя время секундомером. Построить график зависимостиf=f(t) .
4. Смочить водой помещенный рядом со стойкой влажный термометр и по
установившимся показаниям определить температуру T 0. Пользуясь справочным материалом, прилагаемым к работе, найти значенияP 0 иρ .
ρ 1/ 3 | |||||||||||||||
D = − | |||||||||||||||
Скорость изменения частоты | (находится из | обработки | экспериментальной |
||||||||||||
зависимости f=f(t) методом наименьших квадратов;S - чувствительность весов (указана на экспериментальной установке).
6. Весь экспериментальный материал лучше представить в виде таблицы.
4. Контрольные вопросы
1. В чем состоит сущность явления диффузии в газах, твердых телах, жидкостях?
2. Дайте определение стационарным и нестационарным процессам диффузии.
3. В чем состоит физический смысл коэффициента диффузии? Может ли коэффициент диффузии быть отрицательной величиной?
4. Как зависит коэффициент диффузии газов от давления и температуры? Подумайте, как лучше всего с физической точки зрения организовать диффузионный процесс?
5. При испарении капли происходит понижение ее температуры. Почему? Что будет происходить при конденсации пара в капле?
Список литературы
1. Кикоин А.К., Кикоин И.К. Молекулярная физика. М.: Наука, 1976.
2. Матвеев А.Н. Молекулярная физика. М.: Высш. шк. 1981.
3. Варгафтик Н.Б . Справочник по теплофизическим свойствам газов и жидкостей.
В общем случае коэффициент диффузии зависит от следующих факторов:
температуры процесса;
давления в реакционной системе;
концентрации легирующих веществ;
кристаллического совершенства и наличия упругих напряжений в структуре.
2.14 Первый закон Фика
Если рассматривать диффузию как процесс переноса атомов в неравновесной системе, вызываемый только наличием градиента концентрации вещества, то выражение для потока вещества, проходящего в единицу времени через единицу площади сечения, перпендикулярного направлению перемещения вещества, выглядит следующим образом:
J = -D-gradC, (5.1)
где J - поток диффундирующего вещества; D - коэффициент пропорциональности, называемый коэффициентом диффузии и имеющий размерность см 2 /с; С - концентрация атомов вещества, см" 3 .
Знак "минус" в уравнении означает, что поток атомов (и, следовательно, сам процесс диффузии) направлен в сторону убывания концентрации примеси. Коэффициент диффузии - величина, экспоненциально зависящая от температуры
где Ea - энергия активации диффузии, эВ; D 0 - коэффициент диффузии при бесконечно большой температуре, см 2 /с; k постоянная Больцмана, 8,6*10 -5 эВ/К; Т- температура, К.
Уравнение (5.1) называют первым законом Фика, по имени швейцарского физика, получившего его, и говорит этот закон о том, что локальная скорость диффузии через сечение единичной площади в единицу времени пропорциональна градиенту концентрации рассоренного вещества, а направлена диффузия в сторону убывания концентрации примеси.
В технологии микроэлектроники используют плоскопараллельные пластины, и глубина проведения процессов диффузии много меньше толщины пластин. Это означает, что для теоретического анализа можно ограничиться одномерным случаем. Если концентрация примеси изменяется лишь в одном направлении, то выражение для потока в некоторой точке с координатой х в момент времени t будет иметь вид:
2.15 Второй закон Фика
С учетом закона сохранения вещества из первого закона Фика можно вывести выражение для второго закона Фика, согласно которому изменение концентрации вещества в некотором объеме должно определяться потоком внутрь этого объема
Уравнение (5.4) называют уравнением диффузии, и его решение для конкретных граничных условий описывает характер распределения диффундирующих частиц в различных точках среды в зависимости от времени.
2.16 Диффузия в полубесконечное тело из бесконечного источника.
Традиционно термодиффузионные процессы в технологии микроэлектроники проводят в две стадии. На первой стадии на поверхности подложки формируют тонкую легированную область. На второй стадии при отжиге структуры в атмосфере, не содержащей примесей, добиваются перераспределения введенных на первой стадии примесей в подложку. Первая стадия традиционного диффузионного процесса соответствует диффузии в полубесконечное тело из бесконечного источника. Вторая стадия диффузии соответствует диффузии в полубесконечное тело из ограниченного источника.
Полубесконечным телом любую монокристаллическую пластину можно считать при условии, что диффузия проводится на глубину, много меньшую толщины подложки. Бесконечным (или постоянным) источником можно считать такое состояние системы, когда количество примеси, уходящей с поверхности в объем полупроводника, равно количеству примеси, поступающей на поверхность. Это означает, что количество примеси вокруг пластины, например в газовой фазе, значительно больше, чем в ее объеме. Граничные условия в этом случае записываются как
С(х) =О при t = O;
С(х) =C 0 при t >О, х = О;
Решение второго закона Фико при данных граничных условиях описывается дополнительной функцией ошибок