Введение. Для снижения пылевыделения в рабочую зону от мест перегрузки порошкообразных материалов технологические емкости оснащаются системой обеспыливания (аспирации) [1, 2]. Эффективность и энергоемкость работы системы [3, 4] зависит от корректно подобранного, в соответствии с характеристиками аспирационного потока, состава оборудования [5, 6] и его эксплуатационного режима [7]. Пневмотранспортная загрузка материала характеризуется повышенной концентрацией пылевых частиц в емкости [8, 9]. Это объясняется значительным увеличением скорости движения и взаимодействием загружаемого потока и складированного (уложенного) материала [10, 11]. Поэтому при проектировании энергоэффективных систем обеспыливания необходимо учитывать изменяющиеся траектории пылевых частиц в полости технологических емкостей.Методология. Для точного учета поведения частиц в поле скоростей пневмотранспортного потока необходимо учитывать, как продольную, окружную составляющие скорости, так и поперечную. Несмотря на ее малость, приводит к изменению траектории по горизонтали. Потому для задания поля скоростей, в качестве исходных уравнений, приняты классические решения автомодельных движений в турбулентных струях [12], преобразованные в систему дифференциальных уравнений в обычных производных. С последующим их решением универсальным математическим пакетом Maple [13].Основная часть. Моделируем поведение частиц пневмотранспортной струе, поступающей в бункер через вертикальный канал диаметром . Для интегрирования исходных уравнений динамики твердых частиц [14], используем поле скоростей, описываемое системой уравнений [15]. Перейдем к безразмерной форме этих уравнений, используя в качестве характерного геометрического размера диаметр бункера .Тогда в соответствии с [12, 13] функции безразмерных составляющих скорости примут вид:                  (1)           (2)                        (3)где ;  и  – безразмерные координаты траектории пылевой частицы. Сгруппируем коэффициенты:            (4)где  – начальная скорость потока; А, С – постоянные автомодельности.С учетом этого систему (1)-(3) перепишем в следующем виде: ,                 (5)а система уравнений, характеризующая поведение частиц примут вид: . (6)где – критерий Фруда,  – критерий Стокса, учитывающие действие силы тяжести и аэродинамических сил соответственно.Результат. По полученным зависимостям составлена расчетная программа для универсального математического пакета Maple. В качестве исходных данных приняты конструктивные характеристики бункера и узла загрузки, физические характеристики потока (рис. 1. а).  а)                                                                                                б) Рис.1. а) листинг расчетной программы; б) графический результат расчета  Выводы. Приведенная методика расчета позволяет определять траекторию расчета пылевых частиц, необходимую для прогнозирования концентрации пылевых частиц в аспирационном потоке. На рисунке 1 б приведен результат графического расчета в среде Maple для твердых частиц различной крупности: 5, 10 и 15 мкм при движении в условиях не полной веерной струи, формируемой в узле загрузки. Траектории соответствуют рациональным настройкам загрузочного устройства, обеспечивающего отклонение потока от строго вертикального падения, при этом направление движения частиц рассматриваемых диаметров на некотором удалении от узла загрузки оседают в полости бункера.Источник финансирования. Грант Российского научного фонда (проект № 18-79-10025).