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

В данной работе проводится сопоставление компьютерных реализаций численных алгоритмов решения уравнений математических моделей динамики газовзвесей с вязкой теплопроводной, невязкой теплопроводной и идеальной несущими средами. Математические модели разработаны в рамках континуальной методики моделирования динамики многофазных сред. В исследовании моделировался часто встречающейся в горной промышленности процесс взаимодействия ударной волны, движущейся из однородного газа в газовзвесь. Актуальность исследования данного течения неоднородных сред связана с экранированием аэрозольными завесами промышленных взрывов. При моделировании для вязкой среды задавались однородные граничные условия Дирихле, для невязкой среды однородные граничные условия Неймана. Уравнения математической модели интегрировались конечно-разностным методом Мак–Кормака. Для преодоления численных осцилляций применялась нелинейная схема коррекции сеточных функций. Программа, реализующая континуальную методику динамики многофазных сред, состояла из блока задания граничных условий, блока, реализующего численное решение, блока учета межфазного взаимодействия. В результате сопоставления численных расчетов математических моделей динамики газовзвеси с идеальной, невязкой теплопроводной и вязкой теплопроводной несущими средами было выявлено, что в процессе движения газовзвеси наибольшее влияние на интенсивность межфазного обмена импульсом оказывает учет вязкости несущей среды газовзвеси.