Примеры разрушительного воздействия кавитации

Водослив ГЭС (диаметр 15.2 м)
[Brennen C.E. Hydrodynamics of Pumps,Oxford University Press.1994. 293 p].

Лопасти гидротурбины

Вращающиеся детали гидронасосов

Гильза цилиндра двигателя

Применение

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

Примеры применения кавитации

Обработка водомазутной смеси
[Геллер C.В. Приготовление водомазутных эмульсий
посредством волновой диспергации
//Новости теплоснабжения.- 2010 г., №04 (116) 2010 г. ]

Эмульгирование

Очистка деталей сложной формы (фильтров)

1. Разработана методика численного исследования ударного воздействия пузырьков на тело, состоящая из трех этапов [1, 6, 11]. На первом этапе методом граничных элементов рассчитывается форма пузырька и кумулятивной струи, поля скоростей и давлений в жидкости. На втором этапе в предположении абсолютной твердости тела и без учета вязкости жидкости (это соответствует достаточно большим пузырькам и телам, по свойствам близким к металлическим) рассчитывается ударно-волновая динамика в струе и окружающем газе. Результатом второго этапа является дискретная динамическая нагрузка на поверхности тела в ходе удара струи. Далее эта нагрузка аппроксимируется сплайнами и аналитическими зависимостями. На третьем этапе полученная аппроксимация нагрузки используется для расчета деформаций поверхности тела и полей напряжения в приповерхностном слое.

Этап 1. Расчет сжатия пузырька у стенки.
Этап 2. Расчет ударной нагрузки на стенку.
Этап 3. Расчет воздействия ударной нагрузки на тело.

Рис.1. Методика расчета ударного воздействия пузырьков на тело.

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

Рис.2. Сжатие пузырька у стенки в условиях экспериментов [Philipp A., Lauterborn W. Cavitation erosion by single laser-produced bubbles // J. Fluid Mech. 1998. V. 361. P. 75–116] до начала тороидальной фазы движения (a) и один из моментов сжатия в условиях экспериментов [Tong R.P., Schiffers W.P., et. al. The role of `splashing' in the collapse of a laser-generated cavity near a rigid boundary // J. Fluid Mech. 1999. V. 380. P. 339-361 ] на тороидальной фазе движения (b). Верхний ряд – эксперимент, нижний – расчет по разработанной методике.

Рис.3. Расширение-сжатие пузырька в воде вблизи стенки в условиях расчетов [Tong R.P., Schiffers W.P., et. al. The role of `splashing' in the collapse of a laser-generated cavity near a rigid boundary // J. Fluid Mech. 1999. V. 380. P. 339-361]. Формы пузырьков в последовательные моменты времени, поля давления и скорости жидкости в окрестности пузырька в один из моментов сжатия: a), c) – результаты по созданной методике, b), d) – результаты [Tong R.P., Schiffers W.P., et. al. The role of `splashing' in the collapse of a laser-generated cavity near a rigid boundary // J. Fluid Mech. 1999. V. 380. P. 339-361].

3. Изучено влияние на динамику пузырька начального расстояния до стенки и отклонения формы пузырька от сферической. Установлено, что при малой начальной несферичности у пузырька, примыкающего к стенке или удаленного от нее на расстояние, не превышающее половины среднего по объему начального радиуса пузырька, образуется кумулятивная струя, направленная к стенке. При этом с увеличением начального расстояния до стенки у сплюснутого пузырька возрастает скорость соударения конца струи с противоположной ей частью поверхности пузырька, а у вытянутого - давление внутри пузырька [8, 9].

Рис.4. Поля скорости и давления в жидкости и формы пузырьков при их сжатии: верхний ряд - в момент удара кумулятивной струи по стенке, нижний ряд - в один из моментов развитой тороидальной фазы движения для изначально слегка сплюснутого е=0.9 (а), сферического е=1 (b) и слегка вытянутого е=1.04 (с) пузырьков; е – отношение полуосей сфероида.

4. Методика расчета сжатия кавитационного пузырька вблизи твердой стенки обобщена на случай его распада на отдельные фрагменты [5].

Рис.5. Сжатие пузырька у стенки в двух случаях его распада на отдельные фрагменты: a) – пузырек касается стенки, b) – находится от нее на некотором удалении. Верхний ряд - формы пузырьков в пять моментов времени, нижний ряд - поля давления (бар) в моменты 4, 5.

5. Разработаны алгоритмы расчета ударно-волновых процессов в сжимаемых средах с подвижной деформирующейся межфазной границей, в том числе типа газ-жидкость [10-22]. Применяется подход без явного выделения межфазной границы. Для сквозного интегрирования уравнений газовой динамики применяется метод CIP-CUP (Yabe T., Wang P.Y. Unified numerical procedure for compressible and incompressible fluid // J. Phys. Soc. Japan. 1991. V. 60. N. 7. P. 2105–2108). Используется адаптивная сетка специального вида (soroban-сетка, Takizawa K., et. al. Computation of free–surface flows and fluid–object interactions with the CIP method based on adaptive meshless Soroban grids // Comput. Mech. 2007. V.40. P. 167–183), которая представляет собой адаптирующийся к особенностям решения набор узлов на параллельных направляющих линиях.

Рис.6. Высокоскоростной удар капли по жидкой прослойке на стенке: a) – результаты расчета по созданной методике, b) – из работы Чижов А. В., Шмидт А. А. Высокоскоростной удар капли о преграду // ЖТФ. 2000. Т. 170. Вып. 12. С.18–27.

6. Исследована начальная стадия ударного воздействия высокоскоростной струи на жесткую стенку с прослойкой жидкости [6, 10]. Струя образуется при сжатии слегка сплюснутого микронного кавитационного пузырька у стенки при комнатных условиях, скорость струи 250 м/c. Установлено, что
- если пузырек касается стенки, то давление на стенке неоднородно с минимумом в центре удара и максимумом на периферии (более чем в 3 раза превышающим давление гидроудара).
- при наличии между пузырьком и стенкой жидкого слоя толщиной меньше 0.03 радиуса струи распределение давления на стенке также имеет вышеописанный характер с максимумом, также превышающим давление гидроудара.
- Для прослоек большей толщины давление на стенке распределено почти однородно и не превышает давления гидроудара.

Рис.7. Ударное воздействие высокоскоростной струи на жесткую стенку с прослойкой жидкости толщиной d=0.0063R (R – радиус струи): a) – изолинии давления p /pwh (pwh = 5 кбар) и векторы скорости в малой окрестности области контакта струи и прослойки (черные жирные линии – межфазная поверхность, белые линии – ее начальное положение), b) – радиальные профили давления на поверхности стенки в 5 моментов времени и огибающая их максимумов (красная кривая). a) соответствует моменту 4 на b).

Рис.8. Изменение максимума даления на стенке и давления в центре области удара (черная кривая) при воздействии на стенку высокоскоростной струи при наличии на стенке прослойки жидкости разной толщины, τ = R/Cs (Cs – скорость звука в жидкости).

7. Разработана методика расчета динамики приповерхностной области металлического тела при заданном законе нагружения. Тело - изотропное идеальное упруго-пластическое, деформации и перемещения малы. Динамика тела описывается линейными уравнениями, пластичность учитывается с использованием условия Мизеса.

Рис. 9. Распространение в теле импульсов продольных напряжений σz в два момента времени (кривые 1 и 2): (a) – без превышения предела текучести материала, (b) – с превышением предела текучести в переднем фронте импульса, (c) – с превышением предела текучести в переднем и заднем фронтах импульса.

8. Исследована динамика приповерхностного слоя металлического тела при приложении нему нагрузки, характерной для удара кумулятивной струи, образующейся при схлопывании пузырька у поверхности тела, а также при различных аппроксимациях и упрощениях такой нагрузки [11, 24-28]. Показано, что
- наибольшие напряжения возникают как на поверхности тела, так и на удалении от нее порядка радиуса круговой области нагружения [24, 25].
- предельная интенсивность напряжений в теле (на уровне предела текучести) достигается при значительно (примерно в два раза) меньших скоростях струи, чем в соответствующем одномерном приближении [27].
- наибольшее влияние на величину напряжений и остаточных деформаций из характеристик материала тела (модуль Юнга, коэффициент Пуассона, плотность, предел текучести) оказывает модуль Юнга (в частности, чем больше значение модуля Юнга, тем меньше глубина вмятины), а наименьшее – его плотность [28].
Для тел из алюминиевого и никелевого сплавов при различных скоростях струи определены размеры и конфигурация зон, где возникали напряжения с предельной интенсивностью [26, 27].

Рис.10. Поле интенсивности напряжений в теле σz (a) в один из моментов его нагружения по закону, характерному для удара кумулятивной струи, образующейся при схлопывании пузырька у поверхности тела, и, более подробно, соответствующая деформация его поверхности (b, масштаб тела по вертикальной оси увеличен в 104 раз). На рис.10a черная линия – граница области со значениями на уровне предела текучести.