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

К параметрам, характеризующим интенсивность гидродинамического воздействия, следует отнести интенсивность динамического излучения:

1,                                                    (1)

где  Е – энергия импульса давления, генерируемого в канал статора, Дж; t1 – время совмещения канала ротора с каналом статора, с;  Smax=ah ; a, h – ширина и высота прямоугольного канала, м.

Энергия импульса давления определяется по формуле:

2,                              (2)

где c – скорость звука в жидкости, м/с; ρ – плотность жидксти, кг/м3.

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

3,                                                (3)

где P - давление на входе в канал статора, Па; P0 - давление насыщенных паров обрабатываемой жидкости, Па; V – скорость потока жидкости на входе в канал статора, м/с2.

Число импульсной кавитации показывает степень развитости кавитационных явлений в канале статора аппарата. При С<1 кавитация считается развитой. Чем меньше С, тем интенсивнее кавитация.

Из формул (2)-(3) следует, что интенсивность гидродинамического и кавитационного воздействия зависит в основном от скорости и ускорения потока жидкости на входе в канал статора, а также от давления в канале статора при неизменных геометрических параметрах РИА (радиус ротора, ширина и высота каналов в роторе и статоре) и физических свойствах жидкости (плотность, вязкость и скорость звука).

Движение потока жидкости в каналах ротора и статора описывает трёхмерное уравнение Навье–Стокса совместно с уравнением неразрывности для несжимаемой жидкости:

4  (4)

Здесь Vx, Vy и Vz – проекции скорости на оси координат x, y и z.

Сложность гидромеханических процессов не позволяет использовать уравнение Навье–Стокса для нахождения кинематических и динамических характеристик потока жидкости (Vx=Vx(x, y, z, t), Vy=Vy(x, y, z, t), Vz=Vz(x, y, z, t), P=P(x, y, z, t)) в каналах ротора и статора РИА аналитически, так как описание нестационарных гидромеханических процессов в РИА связано с использованием нелинейных дифференциальных уравнений в частных производных с переменными начальными и граничными условиями, многие из которых обычно неизвестны. Особую трудность вызывает то, что одна из границ (поверхность стенки ротора и выходной периметр канала ротора) подвижна.

В настоящее время для компьютерного моделирования турбулентных течений, реализующихся в большинстве программных приложений, наиболее широко распространенной методикой является применение осредненных по Рейнольдсу уравнений Навье–Стокса, дополняемых той или иной моделью турбулентности. Имеется обширная литература по предсказательным возможностям RANS-моделей турбулентности применительно к двумерным течениям, для которых в настоящее время относительно легко обеспечить сходимость решения по сетке. Однако, для трехмерных турбулентных течений, и особенно для тех, которые характеризуются сложной геометрией и многовихревой структурой, получение сеточно-независящих численных решений представляет сложную задачу.

Для решения системы дифференциальных уравнений (4) использовался программный продукт на основе конечно-элементного анализа – ANSYS CFX компании ANSYS Inc. Комплекс ANSYS CFX – это универсальная программная система конечно-элементного анализа, работающая на основе геометрического ядра Parasolid . Модуль ANSYS CFX использует вычислительный блок по давлению, что позволяет применять эти продукты для решения широкого круга инженерных задач. ANSYS CFX предназначен для моделирования ламинарного и турбулентного течения ньютоновских и неньютоновских сред, как однофазных, так и многофазных потоков, на дозвуковых, так и на сверхзвуковых режимах. ANSYS CFX полностью встроен в расчетную среду ANSYS Workbench графическую оболочку, в которую интегрированы все расчетные модули и инженерные инструменты компании ANSYS . 

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

Расчет РИА в ANSYS CFX состоит из двух этапов:

  1. Создание 3D модели для расчета (проектирование и моделирование области пространства, в которой происходят исследуемые процессы);
  2. Просчет 3D модели (задание точности расчета, граничных и начальных условий данной модели).

На первом этапе исследуемая модель может быть создана любым редактором векторной графики:  Autodesk Inventor, Autodesk AutoCAD, Kompas, SolidWorks, T-Flex CAD и др.

Основные требования к создаваемой 3D модели:

  • 1) модель должна обладать некоторым внутренним объемом, то есть представлять собой некоторую твердотельную область пространства с ненулевым объемом, но не объект образованный некоторой поверхностью (оболочкой) с нулевой толщиной.
  • 2) модель должна быть в соответствующем формате, адаптивном для среды ANSYS, то есть файл модели, созданный в векторном графическом редакторе должен иметь соответствующее расширение, в данном случае *.sat. При необходимости файлы с другими расширениями можно конвертировать в вышеуказанный формат при помощи программы CADFix.  На рис. 3 показана 3D модель РИА, интегрированная в ANSYS.

5

Рис. 3. 3D модель внутренних объемов РИА в ANSYS.

Используя данную методику, можно определить численное значение параметров обрабатываемой среды (скорости и давления) в любой точке пространства ротора и статора.

В программном продукте ANSYS CFX был произведен расчет скорости и давления жидкости в области пространства, ограниченной трехмерной моделью внутреннего объема РИА. Задача решалась для следующих начальных условий:

t = 25, °С – начальная температура обрабатываемой жидкости (вода);

Pр = 4·105 Па –давление на входе в канал ротора.

Pс = 4·104 Па –давление на выходе из канала статора.

На рис. 4 представлены примеры результата расчета полей скорости и давления в каналах ротора и статора РИА. 

Для расчета интенсивности излучения необходимо рассчитать средние значения скорости и давления потока по плоскости на входе в канал статора в момент времени, когда скорость потока жидкости максимальна. Это позволяет сделать специальный модуль программы. По средним значениям давления и скорости потока на входе в канал статора определяем  величины интенсивности излучения и числа гидродинамической кавитации по формулам (1)-(3)

6

Рис. 4. Поле давления и скорости в канале ротора и канале статора при их частичном совмещении.

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

 Далее - Компьютерное моделирование гидродинамики потоков жидкости статического гидродинамического кавитатора >>