Книжная полка Сохранить
Размер шрифта:
А
А
А
|  Шрифт:
Arial
Times
|  Интервал:
Стандартный
Средний
Большой
|  Цвет сайта:
Ц
Ц
Ц
Ц
Ц

Изучение траекторий факелов вредных веществ, распространяющихся в атмосфере промышленных предприятий

Покупка
Основная коллекция
Артикул: 619896.01.99
Самсонов, В. Т. Изучение траекторий факелов вредных веществ, распространяющихся в атмосфере промышленных предприятий [Электронный ресурс] / В. Т. Самсонов // Научно-технический прогресс и охрана труда. Сборник научных работ институтов охраны труда ВЦСПС. - Москва : ПРОФИЗДАТ, 1983. - с. 60-69. - Текст : электронный. - URL: https://znanium.com/catalog/product/467546 (дата обращения: 29.03.2024). – Режим доступа: по подписке.
Фрагмент текстового слоя документа размещен для индексирующих роботов. Для полноценной работы с документом, пожалуйста, перейдите в ридер.
ВЦСПС

ВСЕСОЮЗНЫЙ ЦЕНТРАЛЬНЫЙ ОРДЕНА «ЗНАК ПОЧЕТА» НАУЧНО
ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ОХРАНЫ ТРУДА

НАУЧНО-ТЕХНИЧЕСКИЙ ПРОГРЕСС 

И ОХРАНА ТРУДА

СБОРНИК НАУЧНЫХ РАБОТ ИНСТИТУТОВ ОХРАНЫ

ТРУДА ВЦСПС

МОСКВА ∙ ПРОФИЗДАТ ∙ 1983

УДК 628.51:614.72

Канд. техн. наук В. Т. САМСОНОВ

(ВЦНИИОТ ВЦСПС, Москва)

ИЗУЧЕНИЕ ТРАЕКТОРИЙ ФАКЕЛОВ ВРЕДНЫХ ВЕЩЕСТВ, 

РАСПРОСТРАНЯЮЩИХСЯ В АТМОСФЕРЕ ПРОМЫШЛЕННЫХ

ПРЕДПРИЯТИЙ

Организованные     производственпература и плотность выбросов ко
ные выбросы в атмосферу вредных
леблются  в широких  пределах     и

веществ на промышленных площадмогут значительно отличаться     от

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

правленные вертикально вверх. Темпри действии на них горизонтально

60

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

Искривление факела под дейст
вием ветра теоретически заканчивается на бесконечной высоте. Однако в практических расчетах принимается, что осевая линия становится параллельной направлению 
ветра, если угол между ними не превышает 10°. Расстояние по вертикали от устья трубы до уровня, 
на котором осевая линия становится условно горизонтальной, считают начальным возвышением факела и в расчетах загрязнения атмосферы промышленными выбросами 
прибавляют к геометрической высоте трубы.

В ряде случаев благодаря искус
ственному увеличению начального 
возвышения 
удается 
значительно 

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

Однако нахождение оптимальных 

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

методик расчета траекторий факелов.

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

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

Для невозмущенных застройкой 

воздушных потоков наиболее обоснованными представляются зависимости, приведенные в работах 
[1, 2]. В частности, формула    для

вычисления 
гидродинамического 

подъема
(м) имеет вид [1]:

(1)

где
— коэффициент; 
— скорость 

газовоздушной смеси в устье трубы, м/с; 
— диаметр устья тру
бы, 
м; 
—скорость 
ветра, 

м/с;
—плотность выброса и ат
мосферного воздуха, кг/м3.

Как показывают опыты, влияние 

зданий на траекторию факелов 
очень велико.

Так, в работе [3] показано, что 

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

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

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

Гидродинамический подъем без 

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

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

от горизонтального направления, 
изгибается.

Если в выбросе содержится 

пыль, то вследствие оседания частиц осевая линия факела с удалением от устья источника непрерывно понижается.

Таким образом, под воздействием 

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

Изучение траекторий факелов в 

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

подъем отсутствовали).

В опытах использовали линейный 

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

В опытах измеряли профили тем
ператур в нескольких сечениях факела. По построенным на схематическом чертеже профилям избыточных температур проводили температурную осевую линию, которая 
совпадала с одной из линий тока. 
Установлено, что всплывания факела при скорости ветра более 
4 м/с не возникает. Опыты проводили при нескольких значениях 
толщины вытеснения пограничного 
слоя набегающего потока.

Путем обработки опытных дан
ных получена формула для вычисления координат осевой линии факела (или линии тока), искривленной вследствие искажения скоростного поля ветра зданиями и сооружениями:

,    (2)

где 
—приращение 
верти
кальной координаты осевой линии 
факела над уровнем эффективной 
высоты

62

трубы, м; 
—эффективная!

высота трубы, равная сумме геометрической высоты и начального 
возвышения 
факела, 
м; 
—

текущее значение толщины вытеснения, м;
—толщина вытеснения 

в сечении потока, проходящем через 
ось трубы, м.

Сопоставление вычисленных
по 

формуле (2) координат осевой линии факела с результатами опытов 
показало удовлетворительное их совпадение. Толщину вытеснения
в 

заданных сечениях потока вычисляли    по    номограмме   рис. 1   (здесь

— начальное значение
в сече
нии потока, проходящем через подветренную кромку крыши, при вычислении
над 
подветренной 
цир
куляционной зоной (вычисляется по 
кривой 3);
— общее обозначение 

длины циркуляционных зон единой, 
подветренной и наветренной; — расстояние от срывной кромки крыши 
здания, м).

Таким образом, формула (2) опи
сывает траекторию факелов, образующихся от линейных и точечных 
источников в случаях, когда конвективно-динамический подъем отсутствует.

На втором этапе изучали траек
торию изотермической струи, развивающейся 
в 
деформированном 

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

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

Источник выбросов представлял 

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

Были выполнены измерения про
филей 
продольной 
составляющей 

скорости воздушного потока, обтекающего модели зданий.

Опыты показали, что измеренные 

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

Рис. 1. Номограмма для вычисления параметра
1 — подветренная зона;

2 — единая зона; 3 — наветренная зона

составляющей 
скорости 
потока 

ветра в слое, в пределах которого 
происходит искривление факела.

Следовательно, чтобы расчетные 

значения гидродинамического подъема факела соответствовали опытным данным, необходимо в формулу (1) ввести поправку на скорость 
невозмущенного воздушного потока.

Для определения этой поправки 

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

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

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

что
в сечении,    проходящем

через источник). Высота обратного 
потока
может быть рассчитана 

по 
графикам 
рис. 
2 
(здесь 

и
—длина

циркуляционной зоны соответственно единой, наветренной и подветренной).

Похожая картина формирования 

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

64

Различие между зоной обратного потока и ветрозащищенной зоной заключается в направлении движения 
воздушных масс: в первой осредненное 
движение 
осуществляется 

навстречу потоку ветра, во второй —
по ветру.

На рис. 2, б, в ветрозащищенные 

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

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

(м):


(3)

где 
— поправка на величину

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

Анализ опытных данных позво
лил заключить, что расчетная скорость равна скорости ветра на высоте 
= 0,4
. Это соотноше
ние получено следующим образом. 
Подставляя в формулу (3) измеренные в опытах значения 
, 

и 
, вычисляли поправку на 

скорость 
и   по   ней   нахо
дили   высоту

. 
Установлено, 
что 
= 

(0,40— 0,42)
.

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

скорость 
невозмущенного 
потока 

ветра):

,
(4)

где
— плотность застройки, рав
ная отношению суммарной длины 
циркуляционных зон на участке от 
наветренной границы застройки   до

Рис. 2. Границы циркуляционных зон (1), обратных потоков (2)
и  ветрозащи
щенных зон (3):

а - единая  зона;  б — наветренная зона;  в — подветренная  зона

рассматриваемого сечения потока 
к полной длине площадки в тех же 
пределах.

Если скорость ветра задана на 

высоте флюгера ( =10 м), то скорость 
невозмущенного 
ветра

(м/с) можно вычислить по формуле:
где
— обозначение выраже
ния в фигурных скобках формулы 
(4)  для  высоты,    соответствующей
значению параметра
=5 в рас
сматриваемом

5     Зак.   158

сечении потока, и для высоты =10 
м у наветренной границы предприятия.

В зависимости (3) поправка 

является неизвестной величи
ной. В связи с тем, что она зависит 
от определяемой величины подъема 

, решить уравнение (3) в явном 

виде не представляется возможным.

Предлагается следующий спо
соб вычисления гидродинамического подъема 
. Вычисляют гидро
динамический подъем факела, рассчитываемый по скорости невозмущенного (равномерного) потока

65

ветра
(м) по формуле:

(5)

Затем путем попыток находят гидродинамический подъем 
из ус
л о в и я
. Д л я  

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

произведение 
было бы 

больше
и 
меньше

(здесь
—
произвольно при
нимаемое 
значение 
гидродинами
ческого подъема, м). Поправка 
вычисляется для высоты 
= 

0,4
. Возвышение
находят пу
тем графической интерполяции, для 
чего 
строят 
график 
в 
ко
о р д и н а т а х  
,
Н а  

э т о м  
графике 
по 
значе
нию
= 
=
находят 

гидродинамический 
п о д ъ е м  

.  Д а л е е  в ы ч и с л я ют 

расчетную скорость ветра
на высоте 

=0,4 
по    формуле

Определив расчетную скорость 

ветра, можем рассчитать координаты траектории конвективно-динамического факела по формуле:

„
(6)

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

м/с2; 
и
— абсолютная тем
пература соответственно среды и выброса, °С.

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

(7)

где
—приращение ординаты 

осевой линии пылевого факела, м;

— скорость седиментации сред
ней по удельной поверхности частицы 
пыли, 
м/с 
[4];
—

расстояние от устья источника 
вдоль ветра, м; и.ф—расчетная скорость ветра на уровне эффективной 
высоты трубы, м/с.

Вертикальные координаты тра
ектории факела, вычисляемые по 
формуле (6), следует сложить с приращением 
координат
, 
опре
деляемым по формуле (2), а для 

66

пылевых факелов — еще и по формуле (7). Окончательная формула для расчета траектории факела 
имеет вид:

(8)

При выводе формул принято, что 

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

Если устье источника находится 

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

Динамическая 
струя 
выброса 

«пробивает» обратный поток, если 

при
.  Прове
рить превышение гидродинамического подъема факела над высотой 
обратного потока можно по формуле (3) при
=0,2 , Если ока
жется, что 
, то рас
чет траектории факела следует выполнять по изложенной методике.

Сопоставление расчетных траек
торий факелов с экспериментальными данными показало удовлетворительное их соответствие.

На основе результатов исследо
ваний разработана инженерная методика расчета координат осевой 
линии факела. Основные положения 
этой методики рассмотрены на конкретном примере. Условия примера 
соответствуют 
условиям 
экс
перимента, выполненного автором в 
аэродинамической трубе ВЦНИИОТ. 
Благодаря этому оказывается возможным сравнить результаты расчета с опытными данными.

Пример.
Источник нулевой высоты с диа
метром устья 4   мм  установлен  в

центре крыши модели узкого здания размерами 50X50X400 мм. Здание обтекается ветром в поперечном 
направлении. Толщина вытеснения
пограничного слоя, набегающего потока ветра равна 22 мм. Скорость 
выброса 
=8 м/с, скорость    ветра

= 3,35 м/с.
Решение. Прежде всего необходи
мо определить размеры циркуляционных зон, найти координаты их 
границ и линий нулевой скорости. 
Здание является узким. На основании характерных размеров здания— длины =8H и ширины b —
= H, с учетом относительной толщины вытеснения
= 0,44 по фор
мулам работы [5] находим длину
единой циркуляционной зоны. Она 
равна 6,58H (здесь H — высота здания, мм).

Далее определяем толщину вы
теснения в сечении, проходящем через источник, необходимую для вычисления скорости ветра на расчетной высоте. На номограмме 
рис. 
1 
по 
значе
нию
=0,50/6,58 = 0,076 нахо
дим
=0,0235*6,58*50 = 7,73 мм.

Рассчитываем возвышение факе
ла: 
,=2,2*4*8/3,35 = 21 мм. 

Затем определяем действительное 
возвышение факела, для чего примем несколько пробных значений
с целью построения интер
поляционного 
графика. 
При
мем
=30 мм. Поправку на ско
рость
на вы
соте     
= 0,430=12 мм       при
= (12,00+7,73)/7,73 = 2,552 и 
=1   вычисляем  по  форму
ле (4):  
=0,711.     Произведе
ние

равно 30*0,711=21,34 мм. 

Положим
=25 мм. Поправка 

равна     0,631;  
=

= 15,8 мм. Примем  
=40 мм. 

к

Рис.    3.        Интерполяционный 
график для   нахождения   гидродинамического   подъема   фа
кела (к Примеру)

Поправка
равна 0,852. Произ
ведение
= 34,1 мм. На ин
терполяционном графике рис. 3 по 
значению 
=21    мм    находим

= 29,95 мм. Расчетную ско
рость ветра 
на      высоте  

=0,4Х

Х29,95 мм = 11,98 мм вычисляем по 
соотношению
. Поправ
ка
, найденная по формуле (4)

при =1 и 
=2,55, соответ
ствует 0,711. Расчетная скорость
равна 3,35*0,711 = 2,38 м/с.

Рассмотрим искривление оси фа
кела на высоте 30 мм от линии нулевых скоростей. На высоте
ось струи уже изогнута, и ее вертикальный импульс незначителен. 
Здесь на осевую линию начинает 
влиять направление воздушного потока, его повышение или понижение под влиянием препятствий.

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

Определим координаты линии ну
левой скорости единой зоны. Для 
этого воспользуемся графиком, приведенным на рис. 2, а.

На схематическом чертеже рис. 4 

наносим контур здания (вид на 
стену, параллельную направлению 
ветра). Откладываем длину единой
5*

циркуляционной зоны от наветренной стены вдоль направления ветра. Затем координаты точек, лежащих на кривых 1 и 2 рис. 2,а, переносим на рис. 4. Для этого проводим вспомогательную линию, соединяющую наветренную кромку 
крыши   с   концом   циркуляционной

67

x H
0,5
1
2
3
4
5
6
6,5

 ,мм 8,46 9,9 13,7 18,7 25,9 35,3 49,3 57,6

Рис. 4. Схематический чертеж к Примеру:

/ — здание; 2 — источник выброса; 3 — граница обратного потока;  4 — гра
ница единой  зоны;   5 — рассчитанная  траектория  факела

зоны (линия ОС на рис. 4). Эта 
линия является осью абсцисс (начало координат в точке О).

При переносе границ с рис. 2, а 

на рис. 4 необходимо умножить ординаты на длину 
единой зо
ны, а абсциссы — на
По 

рис. 4
путем измерения определяем ординаты
границы обратного по
тока от уровня земли.

Затем вычисляем координаты осе
вой линии факела в предположении, что динамический подъем отсутствует. Вычисления выполняем 
для высоты
=30 мм по фор
муле: 
=30(30+
)/(30 
+ 

7,73) = 0,795(30+
).

Получены следующие результа
ты:

Здесь
— превышение ординат 

осевой линии факела над уровнем 

= 30 мм.

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

линии факела с учетом только гид

родинамического подъема по формуле:
(здесь 
— коор
дината, отсчитываемая от оси источника). 
Получены 
следующие 

значения координат:

мм
0
25
75
125
175
225
275
300

мм
0
25
36,5
42,75
47,9
52
55,6
57,1

мм
10
39,6
55,18
60,75
60,8
57,9
52,4
45,4

Здесь ордината
отсчитывается от 

уровня устья трубы.

На рис. 4 нанесена осевая линия 

факела, построенная по результатам вычислений. Сопоставление вычисленных координат этой линии с 
опытными данными показало их 
удовлетворительное соответствие.

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

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

Вы вод ы

1. В настоящее время отсутствует 

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

x H
0,5
1
2
3
4
5       6
6,5

мм 8,46

,46

9,9
,9

13,7
18,77 25,9 35,3 

49,3

49,3
57,6

60 63,5

,5

65
60
50
35
15
0

0
4,6
,6

9,08
8
2,9 ‒ 4 , 1  ‒13,2 ‒21,7

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

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

3. 
Положительные
результаты 

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

расположение источника по отношению к зданиям и др.). Учет положения траектории факела в пространстве в расчетах загрязнения 
атмосферы 
позволит 
обоснованно 

принимать соответствующие инженерные решения по обеспечению 
чистоты приточного воздуха на промышленных предприятиях.

Л И Т Е Р А Т У Р А

1. И в а н о в
Ю.   В.    Экспериментальное   исследование   струй, развивающихся   

в потоке. — В кн.:  Теория и расчет вентиляционных струй.    Сб. 
трудов.    Л.,    

1965. (ВЦНИИОТ ВЦСПС).

2. Э п ш т е й н   А. М. О форме оси турбулентной струи в неограниченном го
ризонтальном поперечном потоке. — «Инженерно-физический журнал», 1965, т. IX, 
№ 4.

3. Н и к и т и н В. С, М а к с и м к и н а Н.Г.    Унифицированный метод расчета

загрязнения    атмосферы    выбросами    низких    источников— «Водоснабжение   и   санитарная техника», 1980, .№ 4.

4. С а м с о н о в В.  Т.  Изучение на   моделях выделения и распространения пыли

при перегрузке сыпучих материалов. — В кн.: Совершенствование техники безопасности и производственной санитарии. Сб. науч. работ институтов охраны труда ВЦСПС. М., Профиздат, 1980.

5. Методические указания по расчету концентраций вредных веществ в атмо
сфере предприятий   тракторного   и   сельскохозяйственного   машиностроения, застраиваемых зданиями разной высоты и формы. М, 1981.   (ВЦНИИОТ ВЦСПС).