Введение. Проблемам изучения динамики пылегазовых потоков посвящены труды множества отечественных [1, 2, 3] и зарубежных авторов [4, 5]. Нами был предложен вероятностно-статистический подход к определению коэффициента аэродинамического сопротивления частиц в свободном падающем потоке материала [6]. Свободные потоки сыпучего материала встречаются при загрузках бункеров силосного типа, в которых, в отличие от укрытий [7], стенки бункера находятся на удалении от загрузочного устройства [8].В полифракционных материалах счетное количество мелких, пылевых частиц велико. Например, в гранулометрическом составе продукта мелкого дробления Михайловского горно-обогатительного комбината суммарная доля частиц трех младших фракций размером мельче 100 мкм составляет 12,6 %, однако оценочное количество таких частиц превышает 30 миллиардов, что в 50 раз больше количества частиц остальных фракций [9, 10]. Подобные частицы быстро достигают скорости витания "выпадая" из общего потока, но при значительных массовых долях не могут не оказывать влияния на аэродинамические свойства этого потока.Основная часть. Рассмотрим некоторые аспекты взаимодействия пылевых частиц с частицами многократно более крупными, т.е. поперечный дрейф мелких частиц, при опережении.Поперечный дрейф мелких частиц возможен за счет поперечного вектора аэродинамического сопротивления воздуха. Рассмотрим возможный перенос для случая двух частиц, падающих вдоль параллельных вертикальных траекторий. Крупная частица с диаметром  пусть падает по оси ОХ (рисунок 1) со скоростью равноускоренного движения ; ; ,   (1)где  – начальная скорость падения (при х = 0), м/с.Мелкая частица с диаметром  падает по криволинейной траектории, приближаясь или удаляясь от нее в зависимости от того, находится она выше или ниже перегоняющей ее крупной частицы. Уравнение динамики мелкой частицы зависит как от силы тяжести, так и от аэродинамической силы стоксовского типа. В векторной форме она имеет вид: ,      (2)где  – масса мелкой частицы, кг; – векторы, соответственно скорости воздуха и скорости частицы диаметром , м/с; – коэффициент динамической вязкости воздуха, Па*с.Проектируя это векторное уравнение на оси выбранной системы координат XOY, получим: ;              (3) ,                (4)где ; ; ,    (5) – проекции вектора  на оси OX и OY, м/с; – проекции вектора  на те же оси, м/с.  Рис. 1.  Схема дрейфа мелкой частицы диаметром  при ее падении параллельным курсом относительно вертикального перемещения равноускоренно падающей круглой (тяжелой) частицы диаметром   Будем полагать, что вертикальная составляющая вектора скорости воздуха равна нулю , а горизонтальная составляющая изменяется по закону линейного источника (стока), размещенного в точках оси абсцисс: ,             (6)где q – мощность (расход) источника или стока, м3/с; – скорость выдавливания ("всасывания") воздуха на боковой поверхности вертикального цилиндра диаметром , м/с.Положим, что величина этой скорости пропорциональна скорости падения тяжелой частицы  и изменяется экспоненциально относительно перемещаемого центра тяжести этой частицы в подвижной системе координат X1OY1: ,                    (7) Расчеты обтекания шара (рисунок 2) диаметром мм при м/с показали, что максимальное (в передней части падающей частицы) и минимальное (в задней части этой частицы) значение скорости  составляет: Расчетное значение скорости выдавливания и всасывания на боковой поверхности цилиндра образованной падающей шарообразной частицей диаметром , при этом описывается следующей функцией: ,       (8) где  – относительное расстояние между центрами частиц.Рассмотрим случай, когда мелкая частица в начальном положении (при t = 0) находится в точке  и проекции вектора скорости ее падения составляют: ,                      (9) ,                          (10) а тяжелая частица, падающая равноускоренно, находится при t = 0 в точке , т.е. догоняет мелкую. Составляющие вектора скорости тяжелой частицы: ;                      (11) .               s, м , м/с  Рис. 2.  Изменение скорости выдавливания (левый нижний квадрант) и всасывания (правый верхний квадрант) воздуха при  м и м/с  Сносящая скорость воздуха, обусловленная эффектом "выдавливания" или эффектом "всасывания" при обгоне тяжелой частицы в точке  в силу (6), (7) и (8) определяется следующим соотношением: ,            (12)а в точке : ,                    (13) ,                    (14) где в силу равноускоренного падения тяжелой частицы и равномерного падения мелкой частицы: ,         (15) .             (16)Таким образом уравнение динамики мелкой частицы с учетом принятых допущений можем записать в виде: ,    (17)где  – скорость седиментации частицы. ,                 (18) – плотность частицы, кг/м3, – динамическая вязкость воздуха (в расчетах принимается равной  Па*с).Поперечная скорость воздуха определялась по формуле .Здесь скорость вытеснения  определялась по кусочно-гладкой функции (8).Решение системы четырех уравнений (17) с учетом принятых начальных условий (9) – (11) осуществлялось численно в универсальной математической системе Maple (рисунок 3).Расчеты показали, что при пролете "крупной" частицы диаметром  м вблизи пылевой частицы, двигающейся с постоянной скоростью витания, последнюю относит вначале от "крупной" частицы, а затем в сторону теневого следа "крупной" частицы. Причем с увеличением отношения диаметров этих частиц эффект усиливается и ускоряется.   Рис. 3. Траектории частиц диаметром  при пролете вблизи частицы диаметром  м  Выводы. Изучена динамика пылевых частиц в потоке свободно падающего полифракционного сыпучего материала. Доказано наличие поперечных перемещений мелких, пылевых частиц, двигающихся со скоростью витания, при их опережении крупными частицами. Данный эффект необходимо учитывать при рассмотрении влияния фракционного состава перегружаемого материала на расходы эжектируемого им воздуха, и, как следствие, при определении объемов аспирации на предприятиях горноперерабатывающей отрасти.*Исследование выполнено при финансовой поддержке гранта РФФИ № 16-08-0007а и программы развития опорного университета на базе БГТУ им. В.Г. Шухова.